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SECTION  1 


INTRODUCTION 


This  document  describes,  from  a  programmer's  view,  the  coding, 
function,  and  logic  of  the  Spiral  Decay  computer  program.  The  Spiral  Decay 
program  operates  in  the  SPACETRACK  B-3  Semi-Automatic  Programming  System;  and 
will  accept  sensor  observations,  assign  weights  and  biases  to  the  data,  perform 
a  differential  correction  and  predict  past  or  future  position  and  velocity. 

Section  2  of  this  document  contains  a  functional  description  and 
basic  flow  charts  of  the  program.  Section  3  defines  symbolic  names,  data  and 
tape  formats,  and  memory  assignments.  A  detailed  description  of  each  sub¬ 
routine  used  by  Spiral  Decay  is  given  in  Section  4.  The  Appendix  describes 
the  Adams- Bashf or th  integration  routine  which  is  used  by  Spiral  Decay  for 
ephemeris  integration. 
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SECTION  2 

PROGRAM  DESCRIPTION 


2.1  FUNCTIONAL  DESCRIPTION 

The  Spiral  Decay  Program  contains  three  main  functions.  These  are 
(1)  OBSERVATION  WEIGHTING,  (2)  DIFFERENTIAL  CORRECTION,  and  (3)  PREDICTION. 

This  program  is  designed  to  operate  in  the  Schedule  Tape  Mode  under  control 
of  the  B-3  Executive  Program, 

a.  Observation  Weighting 

The  Observation  Weighting  function  of  the  program  is  to  apply 
corrections  (biases)  to  the  observational  data  used  for  differential  correction, 
and  to  apply  statistical  weights  to  the  data  used  in  the  correction  procedure. 

The  Spiral  Decay  Program  must  be  supplied  with  weight  and  bias 
information  for  each  observation.  If  none  is  available,  that  observation  will 
be  omitted  from  the  differential  correction. 

Within  Spiral  Decay,  there  is  room  for  weight  and  bias  data  for 
thirty  sensors.  There  are  approximately  fifteen  sets  of  sensor  weighting  data 
assembled  within  Spiral  Decay.  The  weighting  data  contained  in  Spiral  Decay 
may  be  changed  and/or  increased  by  introducing  a  "weight  tape"  on  logical 
tape  7. 
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b.  Differential  Correction 


The  differential  correction  function  accepts  observational  data  for 
vhich  weight  and  bias  information  is  available,  and  performs  a  weighted 
differential  correction  on  the  crbit  element  set.  Ephemeris  computation  is 
carried  out  in  a  special  perturbations  variation-of-parameters  formulation. 

This  formulation  numerically  integrates  the  perturbative  accelerations 
influencing  the  parameters  of  an  instantaneous  two-body  reference  orbit. 

The  numerical  integration  is  performed  by  a  sixth-order  Adams- 
Bashforth  integration  routine.  This  routine  works  strictly  in  a  variable 
step  size  mode,  based  on  error  control  parameters. 

Perturbations  handled  in  the  program  include:  zona l,tes serai,  and 
sectorial  harmonics,  atmospheric  drag,  and  solar  radiation  pressure.  These 
effects  may  be  controlled  to  incorporate  any  combination  in  a  particular  case. 

A  correction  can  be  obtained  for  any  or  all  of  the  six  orbital 
elements,  and  may  also  include  the  satellite  mass.  This  latter  parameter 
actually  represents  a  correction  to  the  ballistic  parameter,  C^A/m,  where 
CQ  is  the  drag  coefficient,  A  is  the  satellite  cross-sectional  area,  (mZ) 
and  m  is  the  satellite  mass  (kg.). 

The  correction  process  will  repeat  up  to  twice  the  maximum  number  of 
iterations,  as  specified  by  the  user,  or  until  the  root-mean-square  (RMS)  of  the 
accepted  residuals  has  converged.  Accepted  residuals  are  those  that  are  not 
rejected  on  the  basis  of  an  absolute  maximum  or  a  relative  check  based  on  a 
multiple  of  the  residual  RMS.  When  the  correction  process  has  converged,  the 
epoch  may  be  updated  to  any  revolution  number  or  time. 

c .  Prediction 


The  prediction  function  of  the  program  is  used  to  obtain  future  or 
past  position  and  velocity  data  from  an  input  element  set  or  a  corrected 
element  from  the  differential  correction  section. 

Ephemeris  computation  in  the  prediction  function  is  carried  out  in 
the  same  manner  as  in  the  differential  correction,  with  one  addition: 
ephemeris  integration,  while  a  satellite  is  in  the  final  decay  stage,  is 
performed  using  the  Cowell  equations.  The  Cowell  integration  uses  a  fourth- 
order  Runge-Kutta  integration  scheme  instead  of  Adams-Bashf orth  because  of  the 
need  for  frequent  interval  change. 

In  addition  to  a  prediction  ephemeris  listed  on  hard  copy,  Spiral 
Decay  can:  (1)  produce  a  binary  ephemeris  tape  to  be  used  with  the  XYZLA 
program;  (2)  leave  in  core  five  time  points  for  use  in  a  subsequent  GIPAR  run. 
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2.2  FLOW  DIAGRAMS 


The  following  two  pages  display  the  functional  and  subroutine  level 
flew  diagrams  for  Spiral  Decay.  Additional  flow  diagrams,  written  in  further 
detail,  are  included  in  the  individual  subroutine  descriptions  in  Section  4 
and  the  Appendix  of  this  manual. 

The  following  symbol  conventions  have  been  adopted: 


contains  processing  function 


indicates  a  subroutine  which  is 
documented  within  the  section  or 
in  a  referenced  system  or  publi¬ 
cation 


indicates  a  decision  or  switch, 
including  title  on  the  appropriate 
path  (may  be  more  than  two  exits 
if  a  multiple  switch) 


program  or  flow-diagram  connector 


program  terminal  (start  or  exit' 
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SPIRDEC 


FIGURE  I.  SPIRAL  DECAY  BASIC  FUNCTIONAL  DIAGRAM 


10  U  OM 


SECTION  3 

MEMORY  ASSIGNMENTS 


3.1  SYMBOLIC  NAMES  AND  DEFINITIONS 

This  section  provides  a  listing  and  definitions  of  the  symbolic 
names  which  have  been  created  for  the  Spiral  Decay  Program,  and  are  relative 
only  to  the  coding  of  the  program.  The  symbolic  names  which  refer  to  address 
constants,  input/output  parameters,  command  constants,  and  temporary  storage 
will  not  be  listed  in  this  section,  but  may  be  located  directly  in  the  code 
edit  of  Spiral  Decay. 

Definitions  of  some  of  the  symbolic  names  contain  mathematical 
symbols  which  are  relative  to  the  formulation  of  Spiral  Decay.  Definitions 
of  these  mathematical  symbols  may  be  found  in  Appendix  VII,  Volume  1. 

The  symbolic  names  which  appear  in  the  SPS  Master  Assign  Deck  may 
be  referenced  in  either  Aeronutronic  report  U-1691,  Section  5-69,  or  System 
Development  Corporation  report  TM-LX-38/000/00 ,  Appendix  E. 
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NAME 

CELLS 

DEFINITIONS 

All 

49 

Least  Squares  Matrix  (N  x  N) 

ABRHO 

1 

(  a  •  B  •  p  )/2 

TT 

ACCCNT 

1 

Accepted  Residual  Counter 

ADDT02 

1 

Integration  Step  Size 

ADDT03 

1 

“  Integration  Step  Size 

ADLLEPS 

1 

Criteria  for  Increasing  Integration  Step 
Size 

ADPERRK 

1 

Ratio  of  Runge  Kutta  to  Adams  Bashforth 
Integration  Steps 

ADP 

1 

Control  Flag  for  Integration  Routine 

ADRKSTP 

1 

Indicator  for  Runge  Kutta  Integration 

ADTEST 

1 

Integration  Error 

ADXR1 

\ 

1 

j 

ADXR3 

1 

^  Index  Register  Storage 

1 

ADXR5 

1 

i 

AE 

1 

a  ,  Earth  Radius  =  1 
e 

ALFBUF 

7 

2 

Buffer  of  Powers  of  O'  ,  a  =  -f(sin  i)cos: 

2h 

ALFLG 

1 

a  Residual  Rejection  Indicator 

ALSUN 

1 

a  =  Right  Ascension  of  Sun 

o 

AOVH 

1 

a/H 

AP 

1 

Planetary  Magnetic  Index 
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NAME 


CELLS 


DEFINITIONS 


ASCON2 

1 

6378165  =  meters/e. r. 

ASTK 

1 

One  V,ord  of  -• 

AXGR 

1 

\ 

a 

X 

AXO 

1 

a 

X 

o 

AXPRINT 

1 

a  in  Output  Format 

X 

AX 

1 

a 

X 

AYGR 

1 

a' 

y 

AYO 

1 

a 

yo 

AYPRINT 

1 

a  in  Output  Format 

y 

AY 

1 

a  y 

AZGR 

1 

a' 

z 

AZO 

1 

a 

zo 

AZ 

1 

a 

z 

Bll 

7 

Least  Squares  Matrix  Cl  x  N) 

BFLAG 

1 

Indicator  for  Inclusion  of  Bulge  Perturbatioi 

BIAS AD 

1 

Address  of  Bias  Buffer 

BIBUF 

181 

Sensor  Bias  Buffer 

BNADR 

1 

Address  of  BNBUF  and  BNNBUF 

BNBUF 

7 

Buffer  of  (2n  -  3)  B  ,  B  =  Bessel 

n  n 

n  =  0 . 

functions 

1  j  • » •  6 

BNNBUF 

7 

Buffer  of  (2n  -  1)  B  ,  B  =  Bessel 

n  n 

n  =  0, 

functions 

1  j  *  •  •  6 
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NAME 

BOO 

BPRINT 

B 

BTFLAG 

BT 

BUF 


CELLS 

56 

1 

1 

1 

i 
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DEFINITIONS 

Buffer  Saving  New  Epoch  Elements  and  ICK 
Buffer 

B  in  Output  Format 

Ballistic  Parameter 

Indicator  for  Binary  Tape  Output 

Output  Time  for  Binary  Tape 

Temporary  Buffer  for  Computing  Correlation 
Matrix 


CAPD  1 
CAPR  1 
CAPX  1 
CAPY  1 
CAPZ  1 
CARDS  1 
CDLTAXN  1 
CDLTAYN  1 
CDLTN  1 
CNFLAG  1 

COMPARE  1 
CONI  1 
CONBUF  128 
CONTEST  1 


D 

R  =  Magnitude  of  Vector  to  Sensor 
X 

Y  ^  Location  of  Sensor 
Z 

4 

Indicator  for  "P  Card"  Input 
CAa 

xn 

CAa 

yn 

CAn/n 

Indicator  for  n  Only  Correction  on  First 
Interation 

Counter  for  Subroutine  MARTINI 
-(sin^  i)f/2h 
Storage  for  "P  Cards" 

Criteria  for  D.  C.  Convergence 
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NAME 

CORRDY 

CORRDZ 

COS20M 

COSEO 

COSI 

COSOM 

COSO 

COSPH 

COSPSI 

COSTH 

COSU 

COUNTLR 


CELLS 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


COUNTL  1 
COUNTLX  1 
COUNTR  1 
CRAIG  1 
CSALS  1 
CSNK  16 

DC FLAG  1 
DDGR  1 
DELTAQ  1 
DENOM  1 


DEFINITIONS 


Temporaries  used  by  CORRD  subr. 

Cosine  2 
Cosine  e  +  X 
Cosine  i 
Cosine  x 
Cosine  Q 
Cos  0 
Cos  i 
Cosine  0 
Cosine  u 

Number  of  Elements  to  Correct  -  Left  and 
Right  Address 

Number  of  Elements  to  Correct  -  Left  Address 

Counter  for  MATRIX  Subroutine 

Number  of  Elements  to  Correct  -  Right  Address 

Indicator  for  Printed  Observations 

Cos  ctQ 

Tesseral  Coefficient  Buffer 

Indicator  for  a  Differential  Correction 

•\ 

D 
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NAME 

DFLAG 

DGR 

DLFLG 

DLSUN 

DLTB 

DOTPRX 

DOTPRY 

DOTPR2 

DQFLAG 

DQN 

DR  SDL 

D'JBMAT 

E2VGR 

ECOSV 

ENDT 

E01 

EO 

EPS1 

EPSILON 

F10AV 

FIO 


CELLS 

1 

1 

1 

i 

1 

1 

1  > 
1 

i 

1 

1 

1 

36 

1 

1 

I 

1 

1 

1 

1 

1 

1 


DEFINITIONS 

Indicator  for  Inclusion  of  Drag  Perturbations 

d' 

^  Residual  Rejection  Indicator 
60  =  Declination  of  Sun 

A  b/b 


Temporaries  for  DOTPR  Subroutine 

Indicator  for  Delta  Q  Check 
Maximum  Delta  Q 
8  Residual 

Matrix  for  Increasing  Integration  Step  Size 

2  \ 

-e  v 

p/r  -  1 

End  Time  for  Printed  Prediction  Interval 
Output 

U  =  Used  in  Keplers  Equation 
e 

o 

Convergence  Criteria  for  Iteration  on  Keplet:  ' 
Equation 

e  =  0.00001 

Average  Solar  Radiation  Constant 
Solar  Radiation  Constant 
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definitions 


NAME 

FIRST 

FOBHR 

FOBMIN 

FOBSEC 

FTFLAG 

FX 

GAMMA 

GE 

GEZ 

GITADR 

GIT 

GS 

GSZ 

GU 

GUZ 

H10VZ 

HI 

HAFMAT 

HEAD5 

HEDLIN1 

HEDLIN 

HSUBQP 


Forward  or  Reverse  Integration  Indicator 
Floating  Point  Hour,  Min.,  Sec.  in  Output 


Indicator  for  Prediction  by  Revolution  or 
Temporary  Used  by  DIVFIF  Subroutine 

Reflectivity  Constant 

8e  “  With  Zonal  and  Tesseral  Effects 

g£  -  With  Only  Zonal  Effects 

Address  of  GIT  Buffer 

GIPAR  Output  Time  Buffer 

gg  -  With  Zonal  and  Tesseral  Effects 

8g  "  with  Only  Z^nal  Effects 

8u  "  With  Zonal  and  Tesseral  Effects 
gu  -  With  Only  Zonal  Effects 

250,000  Ft.  in  Earth  Radii 

500,000  Ft.  in  Earth  Radii 

Matrix  for  Decreasing  Integration  Step  Size 

Storage  for  BCD  Page  Heading 

Parameters  for  Prediction  Output  Control 

Perigee  Altitude  (km)  for  Output 


13 


NAME 

CELLS 

DEFINITIONS 

HSUBS 

1 

1/Hp 

HX PRINT 

1 

in  Output  Format 

HYPRINT 

1 

hy  in  Output  Format 

HZPRINT 

1 

in  Output  Format 

ICK 

47 

Interpolation  Buffer 

INELT 

9 

Initial  Element  Buffer 

JBSCHNO 

1 

Input  Option  -  From  SPSJOB  Card 

JMPMOD 

1 

Jump  Table  Modifier  in  CDLTB  Subroutine 

JNSAVE 

1 

Save  Index  Registers 

KAPPB 

1 

B  Modifier  (Altitude  Dependent) 

KAPPA 

1 

Upper  Bound  for  B  Modifier 

KEORTM 

1 

.07436662 

KNTRL 

1 

Indicator  for  Elements  Being  Corrected 

LCOUNT 

1 

Counter  for  MARTINI  Subroutine 

LFLAG 

1 

Output  Indicator  for  East  or  West  Longitude 

LINECNT 

1 

Current  Line  Position  on  Page 

LINECT 

1 

Line  Counter  for  RESOUT  Subroutine 

LJBUF 

7 

Temporary  Storage  for  JNDRAG  Subroutine 

LJDC 

1 

(v/2)  ■  B  •  P  •  .  •  6378165 

TT  D 

LLO 

1 

*o 

LNIOBE 

1 

77:  LOG  10 

10  e 

LOGRHOH 

1 

Log  P  (h,T) 

LOLIMTT 

1 

Minimum  Altitude  Allowed  for  Drag  Perturbs ci 
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NAME 

CELLS 

DEFINITIONS 

LOLMT 

1 

Flag  for  Decay  Corridor 

LPREV 

1 

Temporary  in  LPMCD  Subroutine 

LPRINT 

1 

L  in  Output  Format  (deg.) 

LPSUB 

1 

+  OR  -  2TT  Used  in  LPMGD  Subroutine 

LPTEMP 

1 

Temporary  in  LPMOD  Subroutine 

LSQBUF 

14 

Buffer  for  Matrix  Inversion 

LSQX 

7 

Buffer  of  Delta  Elements 

LTBUF 

6 

"P  Card  4"  Parameter  Storage 

MARTSAV 

1 

Index  Register  Storage 

MATRIXB 

28 

Buffer  for  Correlation  Matrix 

MAX 

1 

Maximum  Acceptable  Residual 

MCOUNT 

1 

Counter  for  Matrix  Subroutine 

M 

1 

Counter  in  MARTINI  Subroutine 

MS2ERK 

1 

.00012649618  -  Conversion  from  m/s  to  er/k 

MU 

1 

M-  =  1 

NEWPAGE 

1 

Line  Counter  for  Page  Control 

NICOLET 

200 

Atmospheric  Density  Table 

NPRFAC 

12 

Buffer  of  1/ (n  +  r) \ ;  0  £  (n  +  r)  £  10 

N 

1 

Counter  in  MARTINI  Subroutine 

NU 

1 

V 

NUX 

1 

V 

X 

NUY 

1 

V 

y 

15 


NAME  CELLS 

NUZ  1 


OANDE 

OBDAY 

OBFLG 

OBMO 

OBSLEFT 

OBSPROC 

OBSREJ 

OB YEAR 

OCOUNT 

OLDRMS 

OLDUZ 

OLINCNT 

OPRTESQ 

ORGDAY 

ORMS 

OVALA 


1 

1 

1 

1 


1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


P10LADR  i 
P1LADR  1 
P2LADR  1 
P3LADR  1 
P4LADR  1 


DEFINITIONS 

v 

z 

Address  of  OBLOC  and  EBLOC 
Day  Number  in  BCD  Output  Format 
Indicator  of  Observed  Quantities 
Month  Number  in  BCD  Output  Format 

Counters  for  PROOBS  Subroutine 

Year  Number  in  BCD  Output  Format 
Temporary  Counter 
RMS  of  Previous  D.C. 

Previous  Value  of  u 

z 

Line  Count  of  Printed  Observations 


Days  from  Beginning  of  Year  to  Epoch 

RMS  of  1st  Pass  for  Output 

Output  Control  for  WILBUR  Subroutine 

10  Scaled  T15 

1  Scaled  T15 

2  Scaled  T15 

3  Scaled  T15 

4  Scaled  T15 
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NAME 

CELLS 

DEFINITIONS 

P55LADR 

1 

53  Scaled  T15 

P5LADR 

1 

5  Scaled  T15 

P6LADR 

1 

6  Scaled  T15 

P8LADR 

1 

8  Scaled  T15 

PAGENO 

1 

Output  Page  Counter 

PAPRINT 

1 

Period  -  for  Output 

PAP 

1 

Planetary  Magnetic  Index  -  for  Prediction 

PARTA2 

1  \ 

PARTA3 

1 

1 

Temporaries  for  CDLTB  Subroutine  -  Contain 

PARTA 

1 

[  Segments  of  fa/  ft  and  6e/  ft 

r 

PARTE 2 

1 

PARTE3 

1 

PARTE 

1  , 

PB 

1 

Ballistic  Parameter  -  for  Prediction 

PBUF 

6 

Buffer  of  P  ,  n  =  0,  1,...,  5 
n 

PDAY 

1 

BCD  Day 

PFIOAV 

1 

Average  Solar  Radiation  Constant  for 
Prediction 

PF10 

1 

Solar  Radiation  Constant  for  Prediction 

PFLAG 

1 

Indicator  for  Prediction  Output 

P  GAMMA 

1 

Reflectivity  Constant  for  Prediction 

PHI 

* 

X 

0 

PKAPPA 

1 

Lower  Bound  for  B  Modifier  -  for  Prediction 
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NAME 

CELLS 

DEFINITIONS 

POBUF 

27 

Buffer  for  Output  of  Observations 

POVH 

1 

p/h 

PPBUF 

6 

Buffer  of  P  n  =  1,...,  5 

n 

PPPFLAG 

1 

Indicator  for  a  Prediction 

PRDSAV 

1 

Indicator  for  New  Epoch  Elements 

PREDBF 

11 

Temporary  Buffer 

PREDFLG 

1 

Indicator  for  Prediction  by  Time  or 
Revolution  Number 

PREV 

1 

Revolution  Number  for  Prediction  Output 

PRTADR 

1 

Address  of  PRT  Buffer 

PRTIML 

2 

Revolution  Number  or  BCD  Date  of  New  Epoch 
(D.C) 

PRTIM 

2 

Revolution  Number  or  BCD  Time  of  Final  Epoch 
(Prediction) 

PRT 

19 

Buffer  of  Prediction  Intervals  for  Output 

PTCOUNT 

1 

Counter  of  Interpolation  Points 

PW 

1 

Indicator  to  Print  Weights  and  Biases 

PYEAR 

1 

BCD  Year 

QQ 

8 

Temporary  Storage 

RANGE 

1 

P 

RBUF 

9 

Buffer  of  R 

n  ,m 

RCNT 

1 

\ccepted  Residual  Counter 

RDOTOR 

1 

f/r 

REJCNT 

1 

Rejected  Residual  Counter 
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NAME 
REJFLG 
RESCNTL 
RES OPT 

REV 

RGFLG 

RGSDL 

RHO 

RKINDT 

RMS 

ROFLAG1 

ROFLAG 

ROVA 

RPBUF 

RPCON3 

RPFLAG 

RPT 

RRFLG 

RRSDL 

RSKNTRL 

RTIMU22 


CELLS 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

8 

1 


DEFINITIONS 

Rejection  Indicator  for  any  Residual 

Control  for  Residual  Output 

Indicator  for  Type  of  Units  of  Printed  Angle 
Residuals 

Revolution  Number 

Rejection  Indicator  for  Range  Residuals 
Range  Residual 
Atmospheric  Density  =  p 
Initial  Integration  Step  Size 
Root  Mean  Square  of  Residuals 

i 

Residual  Output  Flags 
r/a 

Buffer  of  R 

n  m 

Radiation  Pressure  Constant 


1 

1 

1 

1 

1 

1 


Indicator  for  Inclusion  of  Radiation 
Pressure  Perturbations 

Maximum  Number  of  Iterations  for  a 
Differential  Correction 

Rejection  Indicator  for  Range  Rate  Residuals 
Range  Rate  Residual 

Indicator  for  Elements  to  be  Corrected 

Vi  -  u2  Used  in  MARTINI  Subroutine 
z 
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NAME 

CELLS 

DEFINITIONS 

SATEL 

1 

■Satellite  Number  (Left  Justified) 

SAVEO 

1 

Cell  0  Storage 

SAVE3 

1 

Cell  3  Storage 

SAVE  LEM 

1 

Indicator  for  Integrating  to  New  Epoch 

SAVEM 

I 

Temporary 

SAVOBS 

I 

Temporary 

SAVT 

1 

Temporary 

SBUF 

180 

Buffer  of  Sensor  Information 

SCONBUF 

70 

Storage  for  "P  Card"  Images 

SCOUNTR 

1 

Counter  for  MATRIX  Subroutine 

SIGMA 1 

1 

l/ffp  (e.r.)"1 

SIGMA2 

I 

l/9p  (e.r./ke)  1 

SIGMA 3 

l 

1/<?A  (rad)’1 

SIGMA4 

I 

l/o.  (rad)’1 

n 

SIGMAI 

1 

I/O.,  o  =  Weight  for  P,  A,  h,  or  p 

SIGN 

7 

Output  Buffer  for  Standard  Deviat’^ns 

SIN20M 

I 

Sine  2u) 

SINBOTH 

I 

Sin  (Y  +  ?) 

SINCOS 

SKNTRL 

9 

Buffer  of  sin  m  >  and  cos  m  \ ;  m  =  1,2 

Indicator  for  Elements  to  be  Corrected 

SMLGR 

I 

1' 

SNALS 

1 

Sin  o' 

9 
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NAME 

CELLS 

DEFINITIONS 

SNDLS 

1 

Sin  *© 

SOBFLG 

1 

OBLFG  Storage 

SOLDUZ 

1 

OLDUZ  Storage 

SREVF 

1 

Final  Revolution  Number  to  End  Prediction 

SREV 

1 

REV  Storage 

SUNLX 

1 

L 

X 

0 

SUNLY 

1 

L 

ye 

SUNLZ 

1 

L 

TAPBUF 

128 

0 

Buffer  for  Ephemeria  Tape  Output 

TAPCNT 

1 

Current  Address  in  TAPBUF 

TEMPO 

1 

Atmospheric  Temperature  at  Epoch 

TEMPI 

1 

Temporary 

TEMP2 

1 

Temporary 

TEMP3 

1 

Temporary 

TEMPT 

1 

Atmospheric  Temperature  at  Any  Time 

TERMS 

8 

Buffer  of  7  Coefficients  and  1  Residual 

TESRAL 

1 

Indicator  for  Inclusion  of  Tesseral  Harmonic: 

THETA 

1 

0 

VALLINI 

'} 

Control  for  Prediction  Output 

VALLIN 

1  J 
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NAME 

CELLS  DEFINITIONS 

WAND B I 

I 

Address  of  WBUF  and  BIBUF 

WBUF 

15 1 

Buffer  of  Weights 

WOBMARK 

I 

Temporary 

XBDGR 

I 

% 

XDD 

I 

•  • 

X 

XDGR 

I 

x' 

XDTGR 

I 

XKE 

1 

k  =  .07436662 
e 

XLAMD 

I 

X  for  Output  (deg) 

XLSUNT 

1 

L  at  Time  t 

XMPER 

1 

m/e.  r.  =  6378165 

XRDGR 

1 

•x 

*R. 

XRMODA 

t  1 

1 

XRMOD 

i  j 

/  Counters  for  MATRIX  Subroutine 

YBDGR 

1 

YDD 

1 

y 

YDGR 

I 

•\ 

y 

YDTGR 

i 

•  \ 
yD 

YRDGR 

I 

•  N 

yR 

YY 

2 

BCD  Epoch  Date 

Z2BUF 

7 

Buffer  of  Powers  of  z/2,  where 

=  ae 
H 
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i 


] 


NAME 

CELLS 

DEFINITIONS 

ZBDGR 

1 

ZDD 

1 

•  • 

Z 

ZDGR 

1 

z' 

ZDTGR 

1 

ZRDGR 

1 

ZrV 

ZZ 

1 

z  =  ae/Hp 
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3.2  DATA  BLOCK  AND  TAPE  FORMATS 


This  section  contains  a  definition  of  each  of  the  buffers  and 
tapes  used  by  the  Spiral  Decay  Program  which  are  not  common  to  the  B-3 
operating  system. 

Mathematical  symbols  used  in  the  buffer  descriptions  are  defined 
in  Appendix  VII,  Volume  I. 


Buffer  Name 

Related  Subroutines 

Page 

All 

LSQ,  LSQS ,  LSQR 

26 

ALFBUF 

CDLTB 

27 

Bll 

LSQ,  LSQS,  LSQR 

28 

BIFUF 

WEIGHT,  BIAS 

29 

BNBUF 

CDLTBIN,  CDLTB 

30 

BNNBUF 

CDLTBIN ,  CDLTB 

31 

BOO 

SAVICK,  CNTRL,  PCONTRL 

32 

BUF 

MATRIX 

33 

CSNM 

MARTINI 

34 

DUBMAT 

AD BASH 

35 

EBLOC 

PROOBS,  NXTOB 

36 

GIT 

PCONTRL 

37 

HAFMAT 

AD BASH 

38 

ICK 

DIVDIF ,  CNTRL,  PCONTRL 

39 

JBUF 

MARTINI 

40 

LJBUF 

JNDRAG 

41 

LOGRHO 

JNDRAG 

42 

LSQBUF 

LSQS 

43 

LTBUF 

PCONTRL 

44 

MATRIXB 

MATRIX 

45 

NICOLET 

JNDRAG 

46 

PBUF 

MARTINI 

47 

POBUF 

WILBUR 

48 

PPBUF 

MARTINI 

49 

PREDBF 

PCONTRL,  CNTRL 

50 
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Buffer  Name 

Related  Subroutines 

Page 

RBUF 

MARTINI 

51 

RPBUF 

MARTINI 

52 

SBUF 

SETSBUF,  NXTOB 

53 

SIGN 

MATRIX 

54 

SINCOS 

MARTINI 

55 

TAPBUF 

TAPEW 

56 

TBUF 

DIVDIF 

57 

TERMS 

CDLTB,  LSQR,  CMPCF 

58 

W 

AD BASH 

59 

WBUF 

WEIGHT,  GETWGT 

60 

Z2BUF 

CDLTBIN 

61 

Tape  Number 

Contents 

Page 

7 

Sensor  Weights  and  Biases 

62 

12 

Binary  Ephemeris  for  XYZLA 

63 

25 


All  49  Cells 


26 


f 


I  ALFBUF  7  Cells 

\  ALFBUF  is  used  in  the  CDLTB  subroutine  and  contains  powers  of 

|  (or’Vn!  for  n  =  0,  1,  2,  3,  4,  5,  6. 

ALFBUF  +  0  O'  °/0!  =  1 

+  1  a  1/1! 

+  2  a  2/2! 

,  +  3  a  3/3! 

+  4  a  4/4i 

+  5  «  5/5‘. 

+  6  a  6/b! 
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Bll 


7  Cells 


Bll  contains  th^  least  squares  vector  used  to  solve  for  the 
corrections  to  the  elements.  For  each  accepted  observation,  a  set  of 
scalar  differential  coefficients  (C)  and  a  residual  (R)  are  computed. 
The  products,  C*R,  are  summed  for  each  observation,  and  the  result  is 
the  BIT.  oatrix. 


Bll  +  0  E  CAn/n  •  R 


+  1  E  Ca  •  R 

xn 

+  2  E  Ca  •  R 

yn 

+  3  E  CU  •  R 


+  4  E  cn  •  R 


+  5  E  ct  •  R 

+6  E  CAB/b  •  R 


Bll  would  appear  as  above  if  all  seven  elements  were  being 
corrected.  Bll  is  actually  compacted  to  represent  only  those  elements 
being  corrected.  If  only  UQ  and  i  were  being  corrected,  Bll  would  be: 


Bll  +  0  E  CU  •  R 

o 

+1  E  •  R 
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BIBUF  181  Cells 


BIBUF  contains  sensor  bias  information  for  range,  azimuth, 
elevation,  range  rate  and  time.  There  are  about  16  sets  of  sensor’ 
biases  assembled  in  the  program;  these  may  be  changed  or  extended  by 
introducing  weight  and  bias  cards  on  logical  tape  7. 


BIBUF  +  0 

+  1 
+  2 
+  3 

+  4 
+  5 


00000SSS  (BCD  Sensor  Number) 


-B 

(e.r.) 

r 

-B 

(e.r.) 

a 

-B 

(e  .r .) 

n 

-B. 

r 

(e ,r ./ke) 

-bt 

(min.) 

Up  to  29  sets  of  sensor  number, 

Br>  Ba»  Bh>  Bf»  Bx»  followed  by 
l  BCD  word  of  00000ZZZ 


+  180 


Bit  position  0  of  each  word 
VERLORT/PRELORT  sensor. 


containing  sensor  number  is  set  to  1  for  each 
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BNBiT 


7  Cells 


BNBUF  contains  factors  of  B  (z)  which  are  computed  in  subroutine 
CDLTB1N  and  used  in  subroutine  CDLTB . 

BNBUF  +0  Bq  (z) 

+  1  Bl  (z) 

+2  B2  (z) 

+3  3  •  B3  (z) 

+  A  3  •  5  •  B^  (z) 

+  5  3  •  5  •  7  •  B5  (z) 

+  6  3  *  5  •  7  •  9  •  Bfi  (z) 
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BNNBUF  7  Cells 

BNNBUF  contains  factors  of  B  (z)  which  are  computed  in  subroutine 
CDLTBIN  and  used  in  subroutine  CDLTB .  n 

BNNBUF  +0  Bq  (z) 

+1  Bl  (z) 

+  2  3  •  B2  (z) 

+3  3  •  5  •  B3  (z) 

+4  3  *  5  •  7  •  B.  (z) 

4 

+  5  3  •  5  •  7  •  9  •  B5  (z) 

+  6  3  •  5  •  7  •  9  •  11  •  B  (z) 

o 
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BOO 


56  Cells 


BOO  is  used  to  store  the  interpolation  buffer  (ICK)  and  new  epoch 
elements  as  found  during  a  D.C.  This  is  done  to  save  unnecessary  integration 
when  starting  a  prediction. 


BOO  +  0  a 

xn 


+  1 


a 

yn 


+  2  B 


+  3  a 


+  4 


2 

e 


+  5  sfn  i 


+  6 
+  7 

+  8 


P 

n 

time  (min) 


new  epoch 
elements 


+  9  L 

+  10  a 

x 

+  11  a 

y 

+  12  a 

z 

+  13  h 

x 

+  14  h 

y 

+  15  h 

z 


+  16  to  +  55  5  sets  of  t,  L,  a,  h,  taken 

directly  from  ICK+0  to  ICK+39 
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BUF  7  Cells 

BUF  contains  the  standard  deviations  of  the  delta  elements  for 
those  elements  being  corrected. 

If  seven  elements  were  being  corrected: 

BUF  +  0  0  n 

+  1  °  a 

xn 

+  2  a  a 

yn 

+  3  a  U 

o 

+  4  a  n 

+  5  a  i 

+  6  a  B 

If  only  a  ,  fi,  B  were  being  corrected: 

xn’ 

BUF  +  0  a  a 

xn 

+  1  cn 

+  2  a  b 
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CSNM 

16  Cells 

CSNM  contains  the  tesseral  coefficients  as  input  from  P  cards 

7  and 

8. 

CSNM  +  0 

C22 

+  1 

S22 

+  2 

C31 

+  3 

S31 

+  4 

C32 

+  5 

S32 

+  6 

C33 

+  7 

S33 

+  8 

C4X 

+  9 

S41 

+  10 

C42 

+  11 

S42 

+  12 

°43 

+  13 

S43 

+  14 

C44 

+  15 

S44 

34 

DUBMAT 


36  Cells 


DUBMAT  contains  a  6  X  6  matrix  used  to  change  the  difference  table 
in  subroutine  AD BASH  to  represent  a  step  size  twice  as  large  as  the  previous 
step  size. 

This  matrix  is  stored  by  rows  starting  with  the  first  element. 


1 

0 

0 

0 

0 

4 

-4 

1 

0 

0 

0 

8 

-12 

6 

-1 

0 

0 

16 

-32 

24 

0 

0 

0 

32 

i 

00 

o 

0 

0 

0 

0 

1 

64 

33 


I 


EBLOC  11834  Cells 

EBLOC  is  used  for  observation  storage.  The  original  10  words 
per  observation  is  cut  to  6  words  per  observation  by  subroutine  PR00BS. 
Up  to  984  observations  are  stored  chronologically  in  EBLOC  as  follows: 

EBLOC  +  0  Station  I.D.  (BCD,  right  adjusted) 

+  1  time  (minutes  since  epoch) 

+2  p  (range  in  e.r.) 

+3  P  (range  rate  in  e.r./ke) 

+  4  Of  azimuth  or  rt.  ascen.  (rad.) 

+  5  ^  elevation  or  declin.  (rad.) 


Up  to  983  sets  of  Sta.  I.D.,  t, 
P,  P,  Of,  6  followed  by  one  word 
/  of  00000ZZZ 
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t 


GIT  11  Cells 


GIT  contains  the  output  times  requested  on  "P  card  10"  to  L>e 
left  in  core  for  a  subsequent  GIPAR  run. 


GIT  +  0  One  bit  for  each  GIPAR  time  requested, 

right  adjusted  to  T47;  i.e., 

0  0  1111  for  four  output  times. 


GIT  +  1 

+  2 
+  3 
+  4 
+  5 


Up  to  five  points  (minutes  from 
prediction  epoch);  terminated  by 
a  word  of  0  0 


During  initialization,  all  11  cells  of  GIT  are  used  for  input  and 
conversion  of  time  from  date  to  minutes  from  epoch. 
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HAFMAT 


36  Cells 


HAFMAT  contains  a  6  X  6  matrix  used  to  change  the  difference  table 
in  subroutine  ADBASH  to  represent  a  step  size  half  as  large  as  the  previous 
step  size. 

This  matrix  is  stored  by  rows  starting  with  the  first  element. 


HAFMAT  « 

.5 

.125 

.0625 

.0390625 

.02734375 

.0205078125 

0 

.25 

.125 

.078125 

.0546875 

.041015625 

0 

0 

.125 

.09375 

.0703125 

.0546875 

0 

0 

0 

.0625 

.0625 

.0546875 

0 

0 

0 

0 

.03125 

.0390625 

0 

l _ 

0 

0 

0 

0 

.015625 
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ICK  47  Cells 


ICK  is  used  as  a  buffer  containing  the  last  five  sets  of  elements 
produced  by  the  integration  routine.  When  an  output  time  falls  within  the 
time  span  covered  by  the  five  points  in  the  ICK  buffer,  a  set  of  elements 
for  the  output  time  are  obtained  by  interpolation  and  stored  in  the  last  7 
cells  of  ICK. 


ICK  +  0  Time  \ 

+  1  L 


+  2  a 

x 

+  3  a 

y 

+  4  a 

z 

+  5  h 

x 


most  recent  point  obtained 
by  the  integration  routine 


+  6  h 

y 

+  7  h 

z 

+  8 


+  39 


/ 

\ 


) 

t 


4  more  sets  of  T,  L,  a,  h 
in  the  order  they  were 
obtained 


+  40 
+  41 
+  42 
+  43 
+  44 
+  45 
+  4b 


L 


a 

x 


a 

y 


a 

z 


h 

x 


h 

y 

h 

z 


obtained  by  interpolation 
for  any  output  time  T. 
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JBUF 


4  Cells 


JBUF  contains  zonal  coefficients  either  as  assembled  in  the 
program  or  as  input  from  P  card  6. 

JBUF  +  0  J2 

+  1  J3 

+  2  J, 

4 

+  3  Jc 
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U3UF 


7  Cells 


UBUF 


is  used  In  the  JNDRAG  subroutine  as  a  set  of  temporaries. 


used  for  temperature 
computations 


+  4 
+  5 


LJBUF  + 

0 

sin 

2.5 

%/2 

+ 

1 

cos 

2.5 

Tl/2 

+ 

2 

cos 

2.5 

T'2  , 

+ 

3 

'I 

used  as  temporaries 
for  interpolation  to 
find  density  coefficients 


+  6 
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LOGRHO 


16  Cells 


LOGRHO  contains  four  sets  of  atmospheric  density  coefficients 
extracted  directly  from  the  NICOLET  table.  The  coefficients  inserted  in 
LOGRHO  by  subroutine  JNDRAG  are  those  which  correspond  to  altitudes 
surrounding  one  point  for  which  the  density  coefficients  will  be  inti rpolated . 


LOGRHO  + 
+ 
+ 
+ 


0  log  P 

1  A 

2  B 

3  C 


/ 


+  4  log  P  \ 


+  5  A 

+  6  B 


0 


+  7  C 


+  8  log  P 

+  9  A 

+  10  B 


+1 


+  11  C 


point  to  be 
interpolated 


+  12  log  P 

+  13  A 

+  14  B 

+  15  C 


/  +2 


increasing 

altitude 


LSQBUF  14  Cells 


LSQBUF  is  a  buffer 
subroutine  FMAIN.  LSQBUF  is 


for  temporary  storage  as  required  by  Phiico 
used  in  subroutine  LSQS. 
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LTBUF 


6  Cells 


I 
I 
I 

LTBUF  contains  prediction  controls  as  extracted  from  P  card  4.  g 

LTBUF  +0  not  used 

+  1  prediction  output  interval  (floating  point)  | 

+  2  bulge  perturbation  flag  (T47) 

+  3  drag  perturbation  flag  (T47)  ^ 

+  4  radiation  pressure  flag  (T47)  j 

+  5  prediction  output  interval  for  a 

binary  ephemeris  tape  (floating  point) 
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MATRIXB 


28  Cells 


MATRIXB  contains  the  correlation  matrix  of  the  delta  elements 

is  a ° lower  half0^  MATRIXB  ls  computed  by  subroutine  MATRIX  and’ 

is  a  lower  half  matrix  compacted  to  represent  the  number  of  elements 

™  f  6  Correlation  half  ^trix  is  stored  in  MATRIXB  by 

ows  for  the  elements  being  corrected  in  the  following  order  n,  a  , 

yn’  o’  xn 

If  n,  Uo,  fi,  and  i  were  being  corrected  MATRIXB  would  be: 


MATRIXB  +  0 

+  1 
+  2 
+  3 
+  4 
+  5 
+  6 
+  7 

+  8 
+  9 


r  =  1 

n  ,n 


1  U 

o,n 

ru„,  U_  =  1 


n 

-i 

n,  u 

c 

n,  n 


l  ,n 


U 


r . 


i,  n 

‘  i,  t =  1 


In  printed  form  this  would  be: 


n,n 

u0,n 

ru  ,  u 

o’  o 

0,n 

rQ,  U 

O 

rn,  n 

i,n 

I\  n 
l’Uo 

ri,n 

ri,i 
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NICOLET 


202  Cells 


+  200  word  1  |  , 

)  at  h  =  1000  km 

+  201  word  2  J 
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PBUF 


6  Cells 


PBUF  is  used  in  subroutine  MARTINI  for  use  in 

of  zonal  harmonics.  PBUF  contains  values  of  P  for  n  = 

n 


the  computation 
0 ,  1 ,  2  ,  3 ,  4  and 


5. 


PBUF  +  0  PQ 


+  1 


+  2 
+  3 
+  4 


+  5 


1 

U 
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POBUF 


27  Cells 


POBUF  is  used  by  subroutine  WILBUR,  and  contains  the  data  from 
each  observation  card  (one  at  a  time)  so  that  a  hard  copy  listing  of  the 
Input  observations  may  be  obtained. 


POBUF 


+  0 

Satellite  No. 

BCD 

1 

Equipment  type 

BCD 

2 

Station  ID 

BCD 

3 

Accuracy 

BCD 

4 

YYMMDD  (year ,mo . ,day) 

BCD 

5 

Hour 

FI.  Pt, 

6 

Minutes 

FI.  Pt 

7 

Seconds 

FI.  Pt 

8 

Elevation  or  declination 

FI.  Pt 

9 

Range 

FI.  Pt 

10 

Range  rate 

FI.  Pt 

11 

Max.  frequency  shift 

BCD 

12 

Brightness 

BCD 

13 

Maximum 

BCD 

14 

Minimum 

BCD 

15 

Time  interval 

BCD 

16 

Not  used 

17 

Message  No. 

BCD 

18 

Equinox 

BCD 

19 

Not  used 

20 

Observation  number 

BCD 

21 

Not  used 

22 

Right  ascension  indicator 

BCD 

23 

Azimuth  or  Rt .ascen. seconds 

FI.  Pt 

24 

Rt.  Ascen.  minutes 

FI.  Pt 

25 

Rt.  Ascen.  hours 

FI.  Pt 

26 

Observation  type 
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PPBUF 


5  Cells 


PPBUF  is  used  in  subroutine  MARTINI  for  use  in 

of  zonal  harmonics.  PPBUF  contains  values  of  P*  for  n  = 

n 


the  computation 
1,  2,  3,  4,  5 


PPBUF  +  0  P'  =  1 


+  1  P'2 

+  2  P' 

+  3  P' 

4 

+  4  P' 
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PREDBF 


11  Cells 


PREDBF  has  two  uses: 

(1)  Whan  the  weight  and  bias  information  is  printed,  PREDBF 
is  used  to  hold  information  in  output  format. 


PREDBF  + 

0 

o  range 

(km.) 

+ 

1 

azimuth 

(deg.) 

+ 

2 

°  elevation 

(deg.) 

+ 

3 

a  range  rate 

(km/ sec) 

+ 

4 

range  bias 

(km.) 

+ 

5 

azimuth  bias 

(deg.) 

+ 

6 

elevation  bias 

(deg.) 

+ 

7 

range  rate  bias 

(km/sec) 

+ 

8 

time  bias 

(sec .) 

(2)  PREDBF  is  used  to  hold  initial  or  final  conditions  as  used 
in  the  integration  buffer.  This  allows  the  integration 
scheme  to  be  restarted  with  a  minimum  of  computation. 


PREDBF  + 

0 

time 

*  W 

+ 

1 

+ 

1 

A  time 

=  W 

+ 

2 

+ 

2 

last  time 

=  w 

+ 

3 

+ 

3 

direction 

=  w 

+ 

4 

+ 

4 

'  L  ) 

- 

=  w 

+ 

5 

+ 

5 

a 

X 

=  w 

+ 

6 

X 

+ 

6 

a 

y 

=  w 

+ 

7 

y 

+ 

7 

a 

)  or 

z 

=  w 

+ 

8 

z 

+ 

8 

h 

X 

=  w 

+ 

9 

X 

+ 

9 

h 

• 

y 

=  w 

+ 

10 

+ 

10 

z 

• 

z 

=  w 

+ 

11 
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RBUF 


9  Cells 


harmonics 


RBUF  is  used  in  subroutine  MARTINI  for  the  computation  of  tesseral 

.  RBUF  contains  entries  of  R  (U  )  for  n  and  in  up  to  4  . 

n,m  z  r 

R  (U  )  =  Pn,m  (V 

n,m  z  — 2 - - 


VTT?- 


RBUF  +  0  R 

+  1 

+  2  R 

+  3  R 

+  4 

+  5 

+  6  R 

+  7  R 

+  8  R 


1,1 


2,2 


3,1 


3,2 


3,3 


4,1 


4,2 


4,3 


=  1 


4,4 
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RPBUF 


8  Cells 


harmonics 


RPBUF  is  used  in  subroutine  MARTINI  for  the  computation  of  tesseral 

.  RPBUF  contains  entries  of  R’  for  n  and  m  up  to  A. 

n,m 


RPBUF  +  0 

+  1 
+  2 
+  3 
+  A 
+  5 

+  6 

+  7 


R,2,2 


R' 


3,1 


3,2 


3,3 


A,1 


A, 2 


R* 


A, 3 


R' 


A, A 
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SBUF 


151  Cells 


SBUF  contains  sensor  information  which  was  originally  stored  in 
SBLOC .  SBUF  and  BIBUF  occupy  the  same  locations  in  core;  however,  SBUF  is 
not  set  until  BIBUF  is  no  longer  needed.  After  SBUF  is  set  up,  the  SBLOC 
area  is  used  for  temporaries  to  save  space  in  core. 


SBUF  +  0  00000SSS,  BCD  Sensor  Number 

+  1  0 

+  2  X 


+  3  X/cos  0 


+  4  Z 

Up  to  29  more  sets  of 
sensor  no.,  0,  X,  X/cos0, 
^  Z,  followed  by  one  BCD 

word  of  00000ZZZ 

+  150  ' 
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SIGN 


7  Cells 


SIGN  contains  the  standard  deviations  of  the  corrections  to  the 
elements  after  each  iteration.  SIGN  is  used  only  for  output  purposes.  If 
an  element  is  not  being  corrected,  the  entry  in  the  buffer  is  set  to  octal  0 
to  suppress  printing. 


SIGN 


+ 

0 

a 

n 

+ 

1 

a 

a 

xn 

+ 

2 

c 

a 

yn 

+ 

3 

c 

U 

0 

+ 

4 

a 

n 

+ 

5 

a 

i 

+ 

6 

a 

B 
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S INC OS 


8  Cells 


SINCOS  contains  the  sines  and  cosines  of  X£,  2X£,  3X£,  4X£  which 
are  necessary  for  the  computation  of  tesseral  harmonics  in  the  MARTINI 
subroutine . 

SINCOS  +  0  cos  K 

Cj 

+  1  sin  X£ 

+  2  cos  2  X£ 

+3  sin  2  X„ 

Li 

+  4  cos  3  X 

L 

+  5  sin  3  Xr 

E 

+  6  cos  4  X£ 

+  7  sin  4  X 

E 
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TAPBUF 


128  Cells 


TAPBUF  is  used  as  an  output  buffer  to  accumulate  one  block  of 
prediction  points  to  be  written  on  the  binary  ephemeris  tape  (logical  12). 


TAPBUF  +  0 

+  1 
+  2 
+  3 

+  4 

+  5 

+  6 


t  (min) 
x  (e.r.) 
y  (e.r.) 
z  (e.r.) 
x  (e.r./kc.) 
y  (e.r./ke.) 
z  (e.r./ke.) 

Up  to  17  sets  of 

t,  x,  y,  z,  x,  y,  z 


+  126 


0“ 


D 


+  127 


* 
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TBUF 


15  Cells 


TBUF  is  used  in  subroutine  DIVDIF  for  temporary  storage  while 
interpolating  with  a  fourth  order  divided  difference  method.  The  contents 
of  TBUF  are  quite  variable  and  cannot  be  described  as  single  quantities. 

TBUF  +  0  thru  TBUF  +  14  Temporaries 
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TERMS 


8  Cells 


TERMS  contains  the  coefficients  to  he  included  in  the  least  squares 
solution  matrix.  These  coefficients  are  computed  by  CMPCF  or  CORRD  depending 
on  the  type  of  observations.  Only  the  coefficients  related  to  elements  being 
corrected  are  computed;  these  are  packed  into  the  TERMS  buffer  in  the 
following  order:  n,  a 


n,  a  ,  a 

i  , 

u  ,  n,  i, 

*  xn’ 

yn’ 

0 

TERMS  + 

0 

^An/n 

+ 

1 

C. 

Aa 

xn 

+ 

2 

CAa 

yn 

+ 

3 

n 

cf 

o 

+ 

4 

ca  n 

+ 

5 

CAi 

+ 

6 

cab/b 

+  7 


R  (weighted  residual) 


TERMS  is  set  up  as  above  when  all  seven  elements  are  being  corrected 

If  a  and  i  were  not  being  corrected  TERMS  would  be: 
yn 


TERMS  +  0 

+  1 
+  2 
+  3 

+  4 
+  5 

+  6 
+  7 


An/n 

CA 

Aa 

xn 

CAU 

0 

CM) 

CAb/B 

R  (weighted  residual) 
0 
0 
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W  203  Cells 


See  Appendix  for  description. 
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WBUF 


151  Cells 


WBUF  contains  sensor  weighting  information  for  range,  azimuth, 
elevation,  and  range  rate.  There  are  about  16  sets  of  sensor  weights 
assembled  in  the  program;  these  may  be  changed  or  extended  by  introducing 
weight  and  bias  cards  on  logical  tape  7. 


WBUF 


+ 

0 

00000SSS,  BCD  sensor  number 

+ 

1 

1/0p 

(e.r  .-1) 

+ 

2 

l'°A 

(rad.  S 

+ 

3 

l/°h 

(rad.  S 

4- 

4 

l/°p 

(ke/e.r .) 

.  i 

• 

Up  to 

29  sets  of  sensor  number, 

a 

1 fo  l/o  1 /o,,  followed 

. 

by  1  BCD  word  of  00000ZZZ . 

+ 

. 

150 
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Z2BUF 


7  Cells 


Z2BUF  is  used  in  the  CDLTBIN  subroutine  for  temporary  storage 
while  computing  Bessel  Functions  for  AB/B  correction.  Z2BUF  is  used  for 
two  different  purposes,  depending  on  the  magnitude  of  Z/2. 


For  Z/2  <  1 
Z2BUF  +0  1 

+  1  (Z/2)2/l! 

+  2  (Z/2)A/2 1 

+  3  (Z/2) 6/3 ! 

+  4  (Z/2) 8/4 1 

+  5  (Z/2) 10/5l 

+6  b lank 


For  Z/2  ^  1 
1 

2/Z 

(2/ Z) 2 

(2/Z)3 

(2/Z)A 

(2/Z)5 

(2/Z)6 
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Logical  Tape  7 


Position 
Block  1 

Block  2 
to  Block 


Contents 

Identification 


1  to  30  weight  and  bias 
cards  followed  by  an 

ENDSIGMA  card:  stored 
in  code  mode,  12  cards 
per  block. 


FIGURE  3.  SENSOR  WEIGHT  AND  BIAS  TAPE  FORMAT 


62 


Position 


Contents 


t  Time  since  epoch  (minutes) 
x  (earth  radii) 

y  (earth  radii) 

z  (earth  radii) 

x  (e.r./ke) 

y  (e.r./ke) 

z  (e.r./ke) 

17  more  sets  of  t,  r,  t 
followed  by  2  words  of 
fixed  point  0 


Up  to  18  sets  of  t,  jr, 
followed  by  1  word  of 
alpha-numeric  Z's 


r 


FIGURE  4.  BINARY  EPHEMERIS  TAPE  FORMAT 


63 


3.3  MEMORY  STORAGE  DIAGRAM 


EXECUTIVE 


SYSTEM 

FILES 


PROGRAM 


BMEWS 

AREA 


NO.  OF 

ASSIGNMENT  WORDS 


OCTAL 

LOCATION 


r 

< 

k. 

r 


< 


r 


< 


EXECMOD1 

EXECM0D2 

MASTER  ASSIGN 

6707io 

AREA 

00000 

EBLOC 

o400io 

15064 

OBLOC 

50°°io 

31464 

SBLOC 

l388w 

43274 

SPIRDECA 

6708io* 

46050 

BMEWS 

650010 

63234 

"  6708  does  not  include  785  locations  assigned  within  SBLOC  that  are 
used  by  SPIRDECA  after  initialization. 


FIGURE  5.  "EMi'i  i'  ;rni!\Gl,  AG  RAM 
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SECTION  4 


SUBROUTINE  DESCRIPTIONS 


This  section  contains  a  description  of  each  subroutine  used  by  the 
Spiral  Decay  Program.  Each  subroutine  description  lists  the  purpose,  usage, 
subroutines  used,  size,  description;  flow  charts  are  included  when  necessary. 

4.1  ALPHABETICAL  INDEX  TO  SUBROUTINE  DESCRIPTIONS 


Name 

Func  tion 

Page 

AD BASH 

Numerical  Integration  Routine 

70 

ALREC 

Compute  A,  Ds  L  from  * 

71 

ANGSUN 

Compute  cy,  f  of  the  sun 

72 

AZREC 

Compute  A,  D,  L  from  a,  h 

74 

BCD TIM 

Convert  Time  to  BCD  Date 

76 

BEGIN 

Compute  classical  elements  at  Epoch 

77 

BIAS 

Apply  Biases  to  the  Observations 

80 

CALH 

Compute  Height 

81 

CALL! 

Compute  U 

82 

CARDER 

Output  Card  Error  Comments 

83 

CDERIV 

Derivative  Routine  for  Cowell  Integration 

84 

C.DLTB 

Compute  CAB/g 

86 

ODLTBTN 

Initialize  Subroutine  CDLTR 

89 
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Name  _ Function _  Page 


CHECK  Test  for  Convergence  91 
CHECKI  Initialize  Subroutine  CHECK  93 
CHGCWD  Modify  Output  Routine  PRNTMAT  94 
CHGNXN  Setup  Least  Squares  Matrix  95 
CMPCF  Compute  Coefficients  for  p,  At  h,  p  95 
CNTRL  Differential  Correction  Control  9g 
C0EFF2  Control  Coefficient  Computation  105 
COELTS  Convert  Elements  for  Output  107 
COMDEL  Compute  A  L  108 
COMDQ  Compute  Aq  109 
COMRMS  Compute  Root -mean -square  110 
COMXLX  Compute  L  HI 
CORRD  Compute  Range  Rate  Coefficients  112 
CREKT  Apply  Corrections  to  the  Elements  H4 
DELOUT1  Print  Delta  Elements  116 
DELOUT2  Print  Final  Results  H7 
DELTAU  Compute  A  u  Hg 
DERIV  Compute  Derivatives  H9 
DIVDIF  Divided  Difference  Interpolation  120 
DOTPR  Compute  Dot  Products  122 
ELMOUT  Print  Corrected  Elements  after  Convergence  123 
ELM0UT1  Print  Corrected  Element  after  Each  Pass  124 
FTAPEW  Wrap-Up  Ephemeris  Tape  126 
GETSEN  Retrieve  Sensor  Information  127 
GETWGT  Retrieve  Weighting  Information  128 
GIPAR  Store  Information  for  Program  GIPAR  129 
HEAD  Output  Page  Headings  130 
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Name 

F  unction 

Page 

IHEAD1 

Initialize  for  Page  Headings 

131 

INITIAL 

Initialize  for  Integration 

132 

INITL 

Print  and  Process  Cards 

133 

INPUT 

Unpack  and  Validate  P  Cards 

134 

ITAPEW 

Initialize  Ephemeris  Tape 

144 

JNDRAG 

Compute  Drag  Perturbations 

145 

LPART 

Compute  L  Modulo  2 tt 

149 

LSQ 

Clear  Least  Squares  Matrix 

150 

LSQR 

Build  Least  Squares  Matrix 

151 

LSQS 

Solve  Least  Squares  Matrix 

152 

MARTINI 

Compute  Bulge  Perturbations 

153 

MATRIX 

Compute  Correlation  Matrix 

158 

MOD 2 PI 

Modulo  2™  Subroutine 

160 

MOVBUF 

Move  Observation  Buffer 

161 

MOVDAT 

Retrieve  Processed  Observations 

163 

NXTOB 

Retrieve  Observations  in  OBLOC  Format 

165 

OBVEC 

Reformat  Observations 

167 

PAGECON 

Page  Control  for  Output 

169 

PCONTRL 

Prediction  Control 

170 

PHLAH 

Compute  0,  h 

182 

PRERES 

Compute  New  Epoch  Elements 

183 

PRINTW 

Print  Weights  and  Biases 

186 

PROOBS 

Reformat  Observation  Buffer 

187 

PRNTMAT 

Print  Correlation  Matrix 

189 

RDPRES 

Compute  Radiation  Pressure 

190 

RDTSB 

Compute  Range  Rate 

191 

READOBS 

Read  Observations  from  Tape 

192 
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Name 

Function 

Pag 

REJECT 1 

First  Pass  Rejection 

193 

REJECT2 

Range  Rate  Rejection 

194 

RESICK 

Restore  Integration  Buffer 

195 

RES0UT1 

Output  Residuals 

196 

RES OUT 2 

Convert  Time  to  BCD  for  Residuals 

198 

RESOUT4 

Print  Page  Headings  for  Residuals 

198 

RESREJ1 

Residual  Rejection 

199 

RESW 

Restore  Integration  Buffer 

200 

RESWBF 

Restore  Integration  Buffer 

201 

REVSUB 

Update  Revolution  Number 

202 

RHOSB 

Compute  P 

203 

RINEL 

Restore  Initial  Elements 

204 

RR2AHL 

Convert  Elements 

205 

SAVCON 

Save  Card  Images 

207 

SAVELM 

Save  Elements 

208 

SAVICK 

Save  ICK  Buffer 

209 

SAVW 

Save  Integration  Buffer 

210 

SAVWBF 

Save  Integration  Buffer 

211 

SAV5PTS 

Save  Elements  for  Interpolation 

212 

SENLOC 

Compute  Sensor  Location 

213 

SETSBUF 

Store  Sensor  Information 

214 

SETW 

Set  Buffer  for  SUBXYZ 

216 

SINEL 

Save  Initial  Elements 

217 

SORTOB 

Sort  Observations  by  Time 

218 

SUB  OUT 

Output  Prediction  Ephemeris 

219 

SUBOUTI 

Initialize  for  Output 

220 

SUBXYZ 

Compute  Position  and  Velocity 

221 
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Name 

Function 

Page 

TAPE* 

Write  Ephemeris  Tape 

224 

TEMP 

Compute  Epoch  Temperature 

225 

THGRC 

Compute  8 

GR 

226 

WEIGHT 

Read  Weights 

227 

WILBUR 

T. int  Observations 

229 

WSETU? 

Set  Integration  Buffer 

230 

SUBROUTINE 

DESCRIPTIONS 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 

SUBROUTINES: 

STORAGE 
r  'QUIREMENTS : 

METHOD : 


REFERENCE: 


AD BASH 
SPIRDEC 


To  integrate  a  system  of  seven  first  order 

differential  equations  from  X  =  X  to  X  =  X  , . 

n  n+1: 


Yi  "  fi  (X-  V  Y2 1  • 


Y?);  i  -  1,  2, 


After  initialization  (see  "How  to  use  ADBASH"), 
co  start  or  continue  integration  simply, 

O'  JMP  AD  BASH 

O'  +  1H  Error  Return 
O'  +  2H  Normal  Return 


Yl3  Y2  5  Y3 '  *  *  * f  Y7  and  X> 

Y. ,  Y  ,  Y0 ,  . .  . ,  Y_  and  X, 
l  Z  u  / 


for  X  =  X 

n 

for  X  =  X 

n+l. 


Program:  LPART ,  DERIV ,  CDERIV 


520  decimal  or  1010  octal  locations. 

The  Adams -Bashf or th  (A-B)  method,  together  with  the 
Runge-Kucta  (R-K)  method  of  integration  are  employed 
in  this  subroutine.  R-K  is  used  as  a  starting  proce¬ 
dure  so  that  an  adequate  number  of  step-wise  solutions 
may  be  obtained  to  build  a  difference  table  needed  by 
the  A-B  method.  The  interval  of  integration  is 
automatically  varied  to  keep  the  discrepancy  between 
the  integrated  values  and  an  absolute  error  check  within 
prescribed  limits. 

A  detailed  description  of  this  subroutine  is  included  in 
the  Appendix  section  of  this  manual. 


70 


PURPOSE : 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 


SUBROUTINES: 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


ALREC 

SPIRDEC 


To  compute  A,  D,  L  if  a,  are  observed. 
JMP  ALREC 


ALPHA  =  a  (rad) 
DELTA  =  A  (rad) 


ASUBX 

= 

-sin 

(ALPHA) 

ASUBY 

= 

cos 

(ALPHA) 

ASUBZ 

- 

F/O 

DSUBX 

= 

C-sin 

(DELTA)] 

X 

[ASUBY] 

DSUBY 

= 

[  sin 

(DELTA)] 

X 

[asubx] 

DSUBZ 

= 

cos 

(DELTA) 

XLSUBX 

= 

[  cos 

(DELTA)] 

X 

[ASUBY] 

XLSUBY 

= 

i — i 

■ 

o 

o 

cn 

(DELTA)] 

X 

[ASUBX] 

XLSUBZ 

= 

sin 

(DELTA) 

PHILCO 

- 

FSIN , 

FCOS 

15  Cells 

Computes  the  values  listed  under  Output. 
Called  by  subroutine  MOVDAT. 
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PURPOSE: 

To  compute 

the  sun's  position  at  time  t. 

CALL  SEQUENCE: 

JMP  ANGSUN 

INPUT: 

W  +  I 

t  (minutes  since  epoch) 

EPOCH 

days  and  frac  from  1950 

to  epoch 

XLSUNO  = 

Lq  at  beginning  of  year 

c 

Cl 

.985647346  deg/day 

C2 

1.91633  deg 

C3 

C^  -  from  TLC  subroutine 

C4 

2.578909 

C5  = 

.4336428 

OUTPUT: 

XLSUNT  = 

Lq  at  time  t 

SUNLX 

L  CSALS  = 

X 

cos  O'  q 

SUNLY 

SUNLZ 

L  SNALS  = 

y 

L  CSDLS  = 

2 

sin  at 

CDS  F  8 

ALSUN 

O'  SNDLS  = 

sir,  «  e 

DLSUN 

0 

SUBROUTINES: 

PH  IT, CO  -  FSIN  j  FCOS,  FATaN 

STORAGE 
REQUIREMENTS : 

23  Ce I id 

DESCRIPTION: 

Computes  th 

e  following: 

\  "  °1 

(t)  +  L  +  C-„  sin 

c  i. 

o 

57.295. . 

i 

where  t  -  days  since  Jar..  0.0 


ANCSUN 
SPIRDEC 
1  of  2 
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ANGSUN 
SPIRDEC 
2  of  2 


DESCRIPTION: 

(continued) 


1 

r 

2Le  1 

t 

1 

=  -  < 

C  sin 

"M 

57.295... 

4 

57 .295 j 

SUNLX  = 

cos 

“e 

COS 

6e 

SUNLY  = 

cos 

6e 

sin 

“e 

SUNLZ  = 

sin 

s0 
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AZPEC 
5FIRDEG 
1  of  2 


PURPOSE : 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 


SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


To  comp, te  A,  D.  L  if  azimuth  and  elevation 
are  observed. 

.IMP  AZREG 

ALPHA  -  A  (rad) 

DLLrA  =  :i  (rad)  FKIRD  =  0 
COS IH  *  9 

S I  H  -  s  I  n  9 

AS  VEX 
AS  l  'By 
ASVLZ 
us  ex 

DS'JB'f 

DSLSZ 

XLSUEX 

XL-  jp; 

XLS f.Z 

PPTLT  -  FSIN,  EGOS 


o  Ceils 

v  -npures  A,  D,  L  from  the  following  formulation. 
Used  by  sTrr.utine  MOVDAT . 

A^  "s'?  =  sir  (ALPHA) 

AS"-  =  :o?  (ALPHA  > 

AS  2.T  -  I/O 

IT  *.\R  -  -c  j s  (DELTA)  cos  (ALPHA) 

XLS  \H  -  c  s  (DELTA)  sin  (ALPHA) 

XLS  ZH  =  sin  (DELTA) 

DS  XL-  - i n  (DEL FA)  cos  (ALPHA) 

DS.  ■  =  :<r  (DELTA)  sin  (ALPHA) 

PS  ZT  -  r.od  (DEL jA) 

SV\PH  -  sin  0 
COSPH  --  co f  0 


'4 


AZPEC 
SPIRDEC 
2  of  2 


DESCRIPTION: 
(continued ) 


ZSU&X  =  cos  0  cos  9 

ZSUBY  =  cos  0  sin  0  Z 

ZbluZ  =  sin  0  j 


SSUEX  =  sin  0  cos  9 
SSUBY  =  sin  0  sin  9 
SSUBZ  =  -cos  0 

ESUBX  =  -sin  9 
ESl'BY  =  cos  9 
ESUBZ  =  F/0 


XLSCBX 

ss 

L  Vs 

xh  x 

+ 

L 

yh  x 

+ 

L  ,  Z 
zh  x 

XLSl BY 

= 

L 

xh  y 

+ 

L  uE 

yh  y 

+ 

L  .  Z 
zh  y 

XLSUBZ 

~~ 

L 

xh  z 

+ 

L  UE 

yh  z 

L  uZ 
zh  z 

ASUBX 

= 

A  S 

Xt  X 

+ 

A  E 
yt  x 

+ 

A  Z 

Zt  X 

ASUBY 

= 

A  S 
xt  y 

+ 

A  E 
yt  y 

+ 

A  Z 
zt  y 

AST'BZ 

A  S 
xt  z 

+ 

A  E 
yt  z 

+ 

A  Z 
zt  z 

DSYBX 

= 

D  S 
xt  X 

+ 

D  E 
yt  x 

+ 

D  Z 

Zt  X 

DSl'BY 

= 

D  S 
xt  y 

+ 

D  E 
yt  y 

+ 

D  Z 
zt  y 

DSYBZ 

D  S 
xt  z 

+ 

D  E 
y  t  z 

+ 

D  Z 
zt  z 

S 


E 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 


SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


BCDTIM 

SP1RDFC 


Tc  convert  time  (min  since  epoch)  to  BCD  date. 

TMA  (Time  in  min.) 
jMP  BCDTa'M 


EPOCH 

OBYEAR 

OBMO 

OBDAY 

FOBHR 

FOBMIN 

FOBSEC 


days  and  fractions  since  1950 


YY 

M>1 

DD 

F /hours 
F /min 
F/ sec 


FCD 

left  justified 


System  -  DKLOK,  SEPSUB 


17  Cells 

Time  is  converted  to  days  and  added  to  EPOCH.  This 
time  is  converted  to  BCD  (YYMMDD) ,  by  DKLOK  subroutine. 
The  fraction  of  a  day  remaining  is  then  converted  to 
flcat'ngpt.  hours,  min.,  sec..  If  the  seconds  >  59.9 
the  time  is  rounded  to  the  nearest  minute. 
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BEGIN 
SPIRLLC 
1  of  3 


PURPOSE: 

To  compute  initia 
elements . 

CALL  SEQUENCE : 

JMP  BEGIN 

INPUT: 

HXO 

=  h 

xo 

HYO 

=  h 

yo 

HZO 

=  h 

zo 

AX  NO 

=  a 

xno 

AYNO 

=  a 

yno 

XLO 

=  L 

o 

XKERTM 

QJ 

II 

OUTPUT: 

P 

P 

RTP  = 

Y7 

WX 

W 

X 

WY 

W 

y 

wz 

w 

z 

COST  = 

cos  i 

SIM  - 

sin  i 

XINCL  = 

=  i 

SING  - 

sin  n 

XNODEO 

-  n 

ESQ  - 

9 

e  “ 
o 

EO 

e 

o 

> 

o 

ii 

a 

o 

R'A  = 

Vao 

SUBROUTINES : 

Program  -  ARC TAN 

Phi  lco 

-  FSQRT 

STORAGE 

REQUIREMENTS: 

43  Cells 

orbit  conditions  from  the  input 


UO 

U 

o 

XNO 

= 

n 

o 

AXO 

= 

a 

xo 

AYO 

— 

a 

yo 

AZO 

— 

a 

zo 

QO 

= 

qo 
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BEGIN 
SPIRDEC 
2  of  3 


DESCRIPTION: 


-V 


2  2  2 
h  +  h  +  h 
xo  yo  zo 


P 


Yp~ — >  RTP 


W 


x  = 


xo  — >  WX  ,  x  ->  y ,  z 

vr 


cos  i  =  W  — >  COSI 
z 


sin 


v;  ~ 

1  =  *  1-COS  1 


SINI 


i  =  tan 


-I 


sin  i 
cos  i 


•XINCL 


sin 

n  = 

W  / 

sin 

i  — 

->  SINO 

X 

cos 

n  = 

W  / 

sin 

i  COSO 

y 

n  = 

tan  ^ 

["  sin 

JL 

- XNODEO 

cos 

n 

e  =  a 
o 


2  ,  2 
+  a 


xno  yno 


ESQ 


EO 


a  =  p/ -:l-e  — >  AO 
o 


yiT"  - >  RTA 

U  =  L  -Q  >  UO 

o  o 

n  =  k  /a  3/2 — 3*-XN0 
o  o 
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BEGIN 
SPIRDEC 
3  of  3 


DESCRIPTION: 

(continued) 


xo 


yo 


20 


-sin  G  •  cos  i 


+cos  i  *  a 

yno 

sin  i  •  a 

yno 


a  (1  -e) 
o  o 


•  a  +  cos  G  a  - AXO 

yno  xno 

cos  G  +  sin  G  a  - >  AYO 

xno 

- - >  AZO 

- >  QO 
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BIAS 

SPIRDEC 


PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


To  subtract  biases  from  the  observations. 

JMP  BIAS 
(Return) 

BIASAD  =  C/HLT,  EBLOC ;  C/HLT,  BIBUF 
BIBUF  =  buffer  of  biases  stored  by  sensor 
Observations  starting  in  location  EBLOC 

See  description 

None 


22  Cells 

The  sensor  number  is  extracted  from  an  entry  in  the 
observation  buffer  and  BIBUF  is  searched  for  a  match. 

If  a  match  is  found,  the  range,  angle,  and  range-rate 
biases  are  subtracted  from  the  observed  quantities. 

Convert  the  time  of  the  observation  to  minutes  since 
epoch,  subtract  the  time  bias,  and  store  in  same  cell  (2,1). 
If  bit  0  of  the  word  containing  the  station  number  is  1, 
then  the  VERLORT/PRELORT  corrections  are  applied.  Keep  a 
count  of  negative  observation  times  in  the  right  address  of 
MCOUNT  for  the  SORTOB  subroutine.  If  no  match  is  found, 
proceed  to  the  next  observation  without  any  error  indication. 

This  procedure  is  continued  until  the  sentinel  of  Z's  is 
found  at  the  end  of  the  observation  buffer. 
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CALH 

SPIRDEC 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 

SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


To  compute  height  above  sea  level. 
JMP  CALH 
UN  =  U 

z 

R  =  r 

F  =  f  =  1/298.3 
H  =  h  (e.r.) 

(None) 


8  Cells 

H  =  h  =  r-1  +  U  2  F  +  3/2  F2  U  2  (1-U  2) 
z  z  z 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 


SUBROUTINES : 

STORAGE 

REQUIREMENTS 

DESCRIPTION: 


CALU 

SPIRDEC 


To  compute  U. 

TMD  C/HLI,  W  +  5;  C/HLT,  X 


JMP  CALU 

W  +  5 

= 

X 

+  6 

= 

y 

+  7 

= 

z 

+  8 

= 

X 

+  9 

= 

y 

+  10 

= 

z 

R 

r 

X 

=  X 

XDOT 

=  X 

UX  = 

U 

X 

Y 

=  y 

YDOT 

=  y 

UY  = 

UZ  = 

u 

y 

u 

z 

Z 

=  z 

ZDOT 

=  z 

PHILCO 

= 

FSQRT 

12  Cells 

Moves  values  from  W  +  5  thru  W  +  10  to 
X,  Y,  Z,  XDOT ,  YDOT ,  ZDOT . 

Computes  the  following: 

-v  /  2  ,  2  ,  2 

r  =  y  X  +  y  +  z 

U  =  x/r 

x 

TJ  =  y/r 

y 

U  =  i.i  r 

z 
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PURPOSE : 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 
SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


CARDER 

SPIRDEC 


To  output  card  error  comments. 
(4  entrances) 


(1) 

TMA 

JMP 

(Word) 

CARDER 

(3) 

JMP 

CARDERB 

(2) 

TMA 

JMP 

(Word) 

CARDERA 

(4) 

TMA 

JMP 

(Word) 

CARDERC 

Word  «  BCD  at  T47 

XR3  =  location  of  1st  word  of  the  erroneous  card 

Comment  on  hard  copy  -  (see  description) 

System  -  PANT 
Program  -  PAGECON 


44  Cells 

The  subroutine  will  insert  the  given  parameter  into 
one  of  the  following,  corresponding  to  the  entrance: 

(1)  "Error  in  field  ending  in  Col.  XX" 

(2)  "System  expects  card  with  X  in  Col.  SO" 

(3)  (No  Comment) 

(4)  "Card  repeated,  too  many  cards  of  type  X." 

The  program  will  then  output  the  comment  (or  no  comments) 
followed  by 

1  8  16  24  32  40  48  56  64  72  80 

(80  Column  Card  Image - >  ) 

1  Line  Space 

Ce '  1  CARDSW  is  set  4  0  to  indicate  this  subroutine  has  been 
used  . 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


SUBROUTINES : 

STORAGE 

REQUIREMENTS 

DESCRIPTION: 


CDERIV 
SPIRDEC 
1  of  2 


To  compute  the  derivatives  of  position  and 
velocity  elements  with  respect  to  tim;,  using 
the  effects  of  radiation  pressure,  drag  and/or 
bulge  perturbations  as  specified  on  P  Cards. 

JMP  CDERIV 


W  +  1  =  time  since  epoch  (min.) 

5  =  F/0 

6  =  x 

7  =  v 

8  =  z 

9  =  x 

10  =  y 

11  =  z 

And  input  needed  for  subroutines  listed. 


W  187  =  F/0 

188  =  dx/dt 

189  =  dy/dt 

190  =  dz/dt 

191  =  dx/dt 

192  =  dy/dt 

193  =  dz/dt 


XdD  =  x 
YDD  =  y 
ZDD  =  z 
XDGR  =  xx 
YDGR  =  yx 
ZDGR  =  z' 


Program  -  CALL',  MARTINI,  JNDRAG ,  RDPRES 


25  Cells 

This  subroutine  is  called  only  when  executing  the 
Cowell  integration.  The  perturbative  effects  are 
computed  by  tbe  subroutines,  MARTINI,  JNDRAG,  RDPRES. 
The  derivatives  of  the  elements  are  then  computed  and 
stored  in  the  Adams-Bashfor th  buffer  (W  +  187  to  W  + 


193)  . 
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FORMULATION: 


ODER IV 
SPIRDEC 
2  of  2 


XDD 

X 

=  rx/r3 

X 

1 

v: 

N 

XDGR  = 

.  s 
X 

=  x'B  +  X 

V  •  • 

D  '  X 

X 

1 

N 

W  +  187 

188 

= 

F/0 

dx/dt  =  k 

X 

189 

= 

e 

dy/dt  =  kg 

y 

190 

— 

dz/dt  =  k 

• 

z 

191 

_ 

e 

dx/dt  =  k 

#  \ 

X 

192 

= 

e 

dy/dt  =  k^ 

.  n 

y 

193 

= 

dz/dt  =  k 

e 

.  > 
z 
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CDLTB 
SPIRDEC 
1  of  3 


PURPOSE:  To  compute  scalar  differential  expression  for 

B  =  CL  A 


CALL  SEQUENCE: 

For 

range  and  angles’ 

For 

range  rate: 

JMP 

CDLTB 

JMP 

CDLTBA 

TAM 

(CAb/B} 

TAM 

^AB/B^ 

INPUT: 

Output  of  subroutines: 

C0EFF2 

,  CORRD,  CDLTBIN 

OUTPUT: 

SUBROUTINES: 


(A)  reg  =  and  output  listed  under  description. 

Program  -  CORRD 


PLilco  -  FEX 


STORAGE 

REQUIREMENTS:  106  cells 

DESCRIPTION:  This  subroutine  must  be  initialized  by  CDLTBIN  subroutine. 

There  are  two  entrances  to  the  program:  (1)  for  range  and 
angles  -  JMP  CDLTB,  (2)  for  range  rate  -  JMP  CDLTBA . 

Formu lation: 


For  range  and  angles,  if  any  of  n/n  or  a^  or  a^  is  not 
being  corrected,  its  coefficient  will  be  computed  as  follows 


CDLTAXN  =  C 


Aaxn 


=  U  (L.  .  V)  +  R  (L.  .  U) 
xn  l  xn  l 


CDLTAYN  -  C 


Aayn 


=  l  U 


V)  +  R 


w^ere  i  =  A  or  h 

and  L^  is  L,  A  or  D  corresponding  to  the  observed 
quantity . 
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DESCRIPTION: 

(continued) 


CD  LIB 
SPIRDEC 
2  of  3 


If  the  range  rate  entrance  is  taken  and  the  3  elements 
are  not  being  corrected,  the  coefficients  will  be  computed 
by  osing  part  of  the  CORED  subroutine.  The  values  will  be 
stored  in  CDLTN ,  CDLTAXN,  CDLTAYN . 

Hove-er.  if  the  elements  are  being  corrected  for  p,  P,  A, 
or  h.  the  values  will  be  unpacked  from  TERMS,  1ERMS+1,  and 
TERMS+2  respectively. 


Then  compute  the  following. 


C0S2CM  =  cos  2  x 

> 

j 

SIN20M  -  sin  2 


xno 


u 


1  cos‘  2 


if  e  ^  o 
o 


if  e  =  0 
o 


Set  up  buffer  ALFBUF  as  powers  of  0/  with  Ot  defined  as: 
'  =  ~(sin^  i)  (f/2h)  cos  2X 
where  f  flattening  =  ^/298.3 
h  =  perigee  height 

A  IT  B  hr’  -  0f''/0:  =  F/l 

n 

Ai  r  B  .T+N  =  O'  1  <  n  <  6 


rompnte  the  following,  which  will  be  used  to  solve  for 

PAKTA  ■-  B  +  ■+  u  (3-  B2)  +  U*  (3-5B,)  .  .  .+t,6  (3 -5  7-9-UB,' 

o  1  — r  try  3  7-7  0 

2  .  3  .  o  . 


PARTA2  =  <"  B ,  -  1  /  2  •  3  B.U2(3B  -  1/2  3  5B, 
1  2  2  3 


)  .  .  .  .+i^(3-5-7-9B  - 1/2 
4 ^ 


3  •  5  ■  7 ' 9 . 1  IB  ) 

(V 
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DESCRIPTION: 

(continued) 


CDLIB 
SPIRDEC 
3  of  3 


PARTA3  =  (3B  -3  *5B  -1/4*  3  *5 • 7B  )+ « . Ay  (3-5-7B,  3*5*7*9BC 

2  3  4  T;  4  5 

+  1/4-3-5-7-9-1IB  ) 

PARTE3  =  B,  +  *B_  +  <2  3-B.  +  . . .  +  a5  3*5 *7 -9B, 

2  IT  3  IT  6 


PARTE2  =  (B  -  1/2  3-B.)  +  <*(3*B  -1/2  3*5B.)  + _ +  y 

2  3  3  4  3 T 

(3-5-7B  -1/2  3*5*7* 9B  ) 

D  D 


PARTE  =  (3 B  -3-5B.  -  1/4*3*5*7BC)  +a  (3-5B,  -  3*5*7BC 
3  4  5  4  5 


+  1/4'3*5*7*9B.) 

6 


where  I)  powers  of  y  are  in  ALFBUF 

2)  1-3 .  (2N- 1)  B  are  in  BNNBUF 

n 

3)  1-3 _ _  (2N-3)  B  are  in  BNBUF 

n 


Solve  for  ^a/Rt  and  oe/ 6 


5e , 


-  exp  '-z) 
fit 


9  / 

(PARTA)+(c»tan2u)  (PARTA2)+(r'tan2JJ)  (PARTA3) 


=  (PARTE)+(t/tan2a))"(PARTE2)+(ytan2u)‘+(PARTE3) 


nr  (I 

Solve  for  C^B  =  (M-Mq)  (aBDn) 


3C 


An/n 


(fa/^t)  -  (C 


Aaxn  xn 


fie. 


+  CAayn  ayn}  ^ 
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A  'RLE' 


f  rfjse: 


it i--  for  .Dlix'b  s-;bi :  -tins 


CALL  SEQ  3VCE:  J><?  CDL7BIN 


I  NT.  I’: 


SKJS ERL  =  C  OXXXXXXX 
AS'  'NS  =  Tyn3781o5  (-T./er) 

"  r.irdzd  r.r  S  .da\Z  s obr entire 


fc-.t  dfes^rltti'n 


5  BR  ILVES:  fr^gcvn  -  ?  BXiZ.  JNDRAG .  'FITS LB 


S  FORAGE 
REQ..  JREMEN'I  5 : 

dzscrtftefn: 


96  Cells 

‘?£3t  if  '.'rrecting  B;  if  net  exit.  Jse  part  of  JNDRAG 
s  irr.-zeine  no  compete  S  and  log  P  (h .  ,  T)  ,  i  -  1,  2,  3S  4 
and  perigcie  density  (Prr)  . 

Chen  com-oite  density  scale  height: 


HS.Bo  - 


6376.165  log^iO 


,0  (3SZ-6S+Z)  log  P  (h.,T> 

6  L 


“*3Sz-4S-i)  log  P  (VJ> 
2 

+  (3Sz-2S-2)  l.g  P  (h  i) 


-  (3S2  -1)  leg  P  (»*  j  I  ■ 
6  “ 


:  "r  „  n--.ar-t  valves  needed  i.r  CDLTB: 

ABRHC  =  BP-)  /  2 
POvH  =  p/h 

s 

A  (AH  =  a'b 

S 

ZZ  -  z  =  ae/h 
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CDLTBIN 
SFIRDEC 
2  cf  2 


DESCRIPTION: 

(continued) 


Test  2/ 2: 

If  z/2  <  1,  set  Z2BUF  as  follows: 


Z2BuT'  =  F/I 


Z2BUF 


r  =  A 

r 


2r 


r  =  1,  2, 


5 


Then  solve  for  B  , 
n 


B  =  E  (z/2)2r 
n  r^°  (ri)  (n+r) I 


n  =  I,  2 ,  ...  6 


If  z/2  >  1 


Z2B'V  =  I 

1 2  L 

Z2BlT+  .  =  |- 

Z 


i2' 


,  x  =  1,  . . . ,  5 


Set  up  BN BIT  and  BNNBUF: 


e  B 


2\  U  1 


E  L 


TT  Z 


where  I.  =1  and  L  =  L 


n,m  n,m-l 


where  n  =  1,  . . . ,  6 
(m-^'/2-n  m^^+n) 


2mz 


Pie  series  should  be  truncated  when: 
L 

n,  m  >1 

u  . 

n,  m-1 


Whatever  the  value  of  Z/2,  set  switch  JMPTBL  +  3  =  C/JMP, 
PA1;  C/JMP,  SOLPRTA  in  CDLTB  subroutine. 
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P  JR  FOSE : 

CALL  SEQ.ENCE: 
1  NP'  7: 


Oi  TP  T  • 

SI  BRO'j  :TNES : 

STORAGE 

REQ  j  LREMENTS: 

DESCRIPTION: 


CHEER 
SP7  RDr  7 

i  i:  2 


To  check  for  correction  cycle  convergence. 


IMP  CHECK 

RPI 

0  0X3  X  =  No . 

of  iterations 

RCNr 

No.  of  residuals 

RMS 

Current  root  mean 

square  of  residuals 

wLDRMS 

Previous  RMS 

DO?  TAG 

Flag  for  Eq  check 

DON 

EH  max 

C  JNTEST  - 

=■  for  allowable  RMS  change  for  convergence 

See  descrip' 

t  ior  . 

Program  -  CHECK!.  COMDQ 

65  Cells 

r-'is  is  not  a  closed  subroutine  in  all  cases.  If  the 
differential  correction  should  continue.,  the  subroutine 
exit  employs  the  Jump  Register  (closed).  If  the  D.C. 
converges.,  the  exit  is  to  location  CONCOM;  if  the  D.C. 
diverges >  to  location  DIVCOM. 

■'I'1  Ar'  CHSWJ-LH,  tve  second  pass  RMS  must  be  saved  for 
at:  the  end  of  the  D.C.  if  the  first  pass  is 
an  nlv"  correction.  Otherwise,  the  first  nass 
RMS  Is  saved . 

ac  chq- 

7 i  a  and  a  are  not  being  corrected,  there  is  no 
xn  yn 

need  rot  a  EH  check.  If  they  are  being  corrected^ 
and  a  El  test  is  requested,  the  following  procrdure 
is  foii  wed 
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DESCRIPTION: 

(continued) 


CHECK 

spirdf: 

2  of  2 


(a)  A9  is  computed  by  COMDQ  s  Enroctine 

(b)  a  may  increase  oy  any  amount  but  it  may  only 

decrease  by  Aq  max  (from  P  Card  3)  .  If  the 

decrease.  >  Aq  max,  Aa  and  Aa  are  modified 
-  xn  yn 

by  1/2  and  the  test  is  repeated  until  Aq  <  Aq  max. 

(3)  At  CHRMS,  test  the  RMS: 

(a)  Compute  the  change  in  the  RMS.  If  it  is  less 
than  CONTEST  (usually  ,01),then  convergence 
has  occurred. 

(b)  If  the  change  is  greater,  exit  to  CONCOM,  test 
if  the  current  RMS  is  greater  than  the  previous 
one.  If  it  is,  the  correction  is  diverging. 

If  the  RMS  increases  on  4  consecutive  iterations, 
exit  to  DIVCOM. 

(c)  At  CHSW2: 

If  the  change  is  greater  than  CONTEST,  but  the  RMS 
is  smaller  than  the  previous  one,  the  number  of 
iterations  (location  RPT)  is  decreased  by  1  and 
tested  for  0.  If  RPT  /  0,  the  D ,C .  will  continue. 

(d)  At  CCC200: 

If  RPT  =  0  for  the  first  time,  a  comment  will  be 
printed  stating  that  CONTEST  will  be  multiplied  by 
5  and  the  D.C,  will  be  repeated  for  the  same  number 
of  iterations. 

(e)  At  CCC190: 

If  RPT  =  0  for  the  second  time,  the  run  will  be 
terminated  by  a  comment  stating  that  convergence  is 
not  determined. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS 

DESCRIPTION: 


CHECKI 

SPIRDEC 


To  initialize  subroutine  CHECK 
JMP  CHECKI 

S KNURL  =  0  DXXXXXXX 

the  last  7  bits  correspond  respectively  to  B,  i, 

fl,  U  ,  a  ,  a  ,  n 
o’  xn  yn 

If  the  bit  is  =  1,  the  element  is  to  be  corrected. 

CNF LAG  =  1,  then  do  n  only  correction  first. 

See  description. 

Program  -  CHGNXN 

14  Cells 

(1)  Sets  the  following  switches  in  the  CHECK  subroutine: 

CHSW2  =  JMP  CCC200 
CHSW1  =  JMP  CHQ 
GHSW  =  JMP  CHSW+1H 

(2)  Moves  SKNTRL  to  KNTRL.  Tests  if  n  only  correction ,  if  y 

Set  KNTRL  =  0  01  and  CHSW  =  JMP  CHN 

(3)  Sets  OLDRMS  =  F/ 1000000 

DIVFL  =  0  0 

COJNTL  =  C/HLT,  N, where  N  -  number  of  bits  =  !  in  KNIR 

(4)  uses  subroutine  CHGNXN  to  set  up  matrix  size. 
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CL.OJR 

SP1RDEC 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINES: 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


To  modify  the  output  control  words  for 
subroutine  PRNTMAT . 

JMP  CHGCWD 

CGUNTR  =  C/HLT.O;  C/HLT,N 
N  =  matrix  size 
Ls  N  -_1 

Modified  control  words  for  subroutine  PRNTMAT. 
None 


19  cells 

If  the  matrix  to  be  printed  by  subroutine  PRNTMAT 
is  a  1  x  1  or  a  7  x  7,  no  modification  of  the  output 
control  words  is  necessary.  If  the  matrix  is  a  2  x  2 
to  a  6  y.  6,  the  control  words  must  be  modified  to  pick 
up  the  correct  locations  in  buffer  MATRIXB . 

(See  subroutine  MATRIXB  for  definition  of  buffer.) 
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CHGNXN 

:T  r  JHC 


PURPOSE 
CALL  SEQUENCE- 

INPUT- 

OUTPUT: 


SUBROUTINES: 

STORAGE 

REQUIREMENTS 

DESCRIPTION- 


To  set  parameters  for  least  squares  subroutines. 

TMA  C/HLT, N  N  «  matrix  size 
IMP  CHGNXN 

:A,  =  C/HLT  ,N 


COUNTL 
COIN  TP 
CO.NTLR 
■  son 

USQP? 

GSQP3 

LSQPA 

1-SQF5 

LSQP6 

eSQSPl 

LSQSp3 

LSQUN 


fo-.-rr  -  C'HGCWD 


18  Cells 


COUNTL  =  C/HLT, N;  C/HLT, 0 
CO'JNTR  -  C/HLT, 0;  C/HLT, N 
COUNTLR  =  C/HLT,N;  C/HLT,N 

LSQPl  =  C/HLT,  TERMS  +  N;  C/JMP,  LSQMULT 
LSQP2  =  C'HlT  All  +  N2  +  1;  C/JMP,  LSQL0D2 
LSQT-3  =  C/HLT  Bll  +  N;  C/JMP,  LSQLOC 
LSQP 8  -  C/HLT,  N-l 
LSQP6  =  C/HL] ,  N  +  1,  C/HLT,  N  +  1 

1SQSP1  =  C/HLT,  N 

LSQSP3  =  C/HLT  ,N;  C  'HLT.N 

1.SQDN  =  C/TDM,  0,2;  C/.\ix0L,  N',  2 
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CMPCF 
SPIRDEC 
1  of  2 


PURPOSE:  If  range  or  angles  are  observed,  to  compute  the 

residuals  and  differential  expressions  which  are 
entered  in  the  least-squares  matrix. 


CALL  SEQUENCE: 

TJM 

CMPCFY 

for  range 

< 

or 

1 

CM 

CMPCFY 

for  angles 

JMP 

CMPCF 

TAM 

(unweighted 

residual) 

TQM 

(weighted  residual) 

INPUT: 

KNTRL 

=  o - OXXXXXXX, 

The  last  7  bits  correspond  to  the  7  elements 

B,  i,  fi,  U  ,  a  ,  a  ,  n. 
o  yn  xn 

If  =  1,  then  the  element  is  being  corrected. 

LTERMS  =*  C/HLTR,  TERMS 

and  output  of  C0EFF2  subroutine 

OUTPUT:  See  description. 

SUBROUTINES:  Program  -  CDLTB 

STORAGE 

REQUIREMENTS:  41  Cells 


DESCRIPTION: 


The  following  formulas  ar^  solved  for,  as  indicated 
by  KNTRL,  and  packed  into  the  TERMS  buffer  to  be 
entered  in  the  least  squares  matrix  by  subroutine  LSQR. 

If  P  is  observed,  the  terms  are: 


<V 

C -V)  +  (RN)  (L-U) 

1 

aP 

for  n 

<V 

a-Y)  +  <v 

(L-U) 

1 

ap 

for  a 

xn 

<v 

(i-Y)  +  %> 

(L-U) 

1 

1  ^ 

for  a 

yn 

<u  ) 

(L-V)  +  (R  ) 

(L-U) 

1 

for  U 

u 

-  u 

a 

n 

o 
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CMPCF 
SPIRDEC 
2  of  2 


DESCRIPTION: 

A,1 

(cos 

i) 

(L-V) 

-  'sir.  i'  (cos  U)  (L*W) 

r 

for 

0 

(continued) 

L 

r 

_ 

- 

ap 

5,1 

= 

(sin 

U) 

(L-W) 

r 

a 

for 

i 

L 

- 

P 

6,1 

= 

(See 

results 

of  CDLTB  subroutine, 

for 

B 

ic  is  used  to  compute  this  term) 

7,1  -  (P  -  D  )  — weighted  residual 

C  CP 

The  unweighted  residual  is  -  p  )  . 

If  A  or  o'  is  observed: 

In  the  formulas  for  n  through  i,  substitute 

i  1 

A  for  L,  and  — :  for  — — 

—  —  pa.  a 

c  A  p 

The  weighted  residual  is:  p  (A*AL)  -  ; 

c  mcctA 

The  unweighted  residual  is:  P„  (A* AL) 

d 


If  h  or  g  is  observed: 


Ir. 

for 


the  formula  for  n  through 
and  p^T  f3r 

r  •>  n 


i, 


substitute  D 


The  weighted  residual  is:  P  (D*AI)  - - 

pcah 

and  the  unweighted  residual  is  Pc  D-AL)  . 
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CNIRL 
SPIRDEC 
1  of  5 


PURPOSE:  To  control  the  differential  correction. 

CALL  SEQUENCE:  JMP  CNTRL 

JMP  (Error) 

(Normal  Return) 

INPUT:  Observations  starting  in  EBLOC 

OANDE  =  C/HLT,  OBLOC ;  C/HLT,  EBLOC 

XLO 

AXO 

AYO 

AZO  Elements  at  epoch 

HXO 

HYO 

HZO 

> 

Input  required  for  all  subroutines  used. 

Input  from  P  Card  3. 

OUTPUT:  DLTNN  =  An/n 

DLTAX  =  Aa 

xn 

DLTAY  =  Aa 

y 

DLTUO  =  Au 

o 

DLTND  =  A  Q 

DLTIN  =  A  i 

DLTB  =  Ab/b 

RMS  =  root  mean  square  of  weighted  residuals 

SUBROUTINES:  System  -  PANT 

Program  -  PAGECON,  MOVDAT,  INITIAL,  CDLTBIN,  SAV5PTS, 

ADBASH,  REVSUB ,  SAVICK,  SAVWBF ,  DIVDIF,  SAVW, 
SETW,  COEFF2 ,  RESW,  RESOUT1,  SAVELM,  COMRMS , 
LSQS 

STORAGE 

REQUIREMENTS:  106  Cells 
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CNTRL 
SP1RDEC 
2  of  5 


DESCRIPTION- 


This  subroutine  is  primarily  a  logic  routine  to 
control  use  of  the  subroutines  which  are  necessary 
for  a  differential  correction 

The  observations  are  stored  and  retrieved  in  order 
starting  at  epoch.  The  closest  observation  before 
epoch  is  the  first;  then  order  the  observations 
backward  in  tine  to  the  earliest  one  before  epoch, 
The  observations  after  epoch  are  stored  next  in 
chronological  order  Epoch  can  be  anywhere  in 
relation  to  the  observations,  i.e  before,  after, 
or  in  the  middle.  The  integration  control  is  set 
to  handle  each  situation  most  efficiently.  The 
integration  starts  at  epoch  and  goes  back  in  time 
to  the  earliest  observation,  then  everything  is 
reinitialized  and  integration  begins  at  epoch  and 
goes  forward  in  time  to  the  latest  observation. 

If  at  any  time  during  the  integration  the  time  or 
revolution  number  requested  for  the  new  epoch  is 
found  (P  Card  3) .  the  elements  are  saved  to  avoid 
repetition  of  integration.  If  the  new  epoch  falls 
outside  the  span  of  the  observations,  subroutine 
PRERES  will  continue  or  initialize  the  integration 
to  find  the  new  epoch  elements  (See  subroutine 
PRERES) 

A  The  following  are  executed  only  once: 

(1)  Set  these  locations  #  0. 


PFJ  SAV 


indicates  new  epoch  elements 
have  not  been  found 


S3.VFLEM  -  subroutine  SAVELM  has  not 
been  called 

FIRST  -  some  observations  may  be 
before  epoch 


(2)  Force  a  page  for  output 


(3)  Zero  the  least  squares  matrix  buffers 
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CNIRL- 

SPIRDEC 

3  of  5 


DESCRIPTION: 

(continued 


Set  the 

following: 

REJCNT 

=  0/0 

(REJCNT  is  number  of  rejected 
residuals) 

RCNT 

=  0/0 

(RCUT  is  number  of  accepted 
residuals) 

SUM 

=  F/0 

(SUM  is  sum  of  squares  of 
weighted  residuals) 

REV 

=  EPREV 

(EPREV  is  epoch  revolution) 

(5)  Retrieve  the  first  observation.  If  no 
observations,  take  error  exit. 

Otherwise,  determine  whether  the  first 
observation  is  before  or  after  epoch, 
i.e.  should  the  integration  be  backward 
or  forward  in  time  with  epoch  as  a  base? 

Set  necessary  switches  accordingly. 

B.  Initialize  for  the  first  integration: 

(1)  Call  subroutine  INITIAL  to  start  with  the 
epoch  elements. 

(2)  Initialize  the  CDLTB  subroutine  with  CDLTBIN. 

(3)  Use  subroutine  SAV5PT5  to  save  the  epoch 
elements  as  the  first  entry  in  the  inter¬ 
polation  buffer  (ICK) , 

C.  Begin  the  basic  integration  loop: 

(1)  Integrate  for  the  next  point  ising  ADBASH 
subroutine . 

(2)  Update  the  revolution  number  with  REVSUB 
subr  j  it  in  e . 

(3)  Save  the  clement  set  with  subroutine  SAV5PTS. 

(4)  If  the  buffer  for  interpolation  does  not 
contain  5  elements  sets,  go  to  step  C(l). 
Otherwise  go  to  D . 
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DESCRIPTION: 

(continued) 


CNTRL 
SPIRDEC 
A  of  5 


D.  Look  for  new  epoch  elements: 

(1)  If  the  new  epoch  elements  have  been  found, 
go  to  step  E.  Otherwise  determine  if 
searching  for  a  revolution  number  or  a  time. 

(2)  If  b>  revolution,  test  if  current  revolution 
number  is  equal  to  the  new  epoch  revolution 
number.  If  not,  go  to  E.  If  it  is,  save  the 
last  element  set  in  the  ICK  buffer  as  the 
element  set  corresponding  to  the  revolution 
number.  Save  other  values  necessary  for 
PRERES  subroutine  and  set  PRDSAV  =  0  indicating 
elements  have  been  found. 

(3)  If  by  time,  test  if  new  epoch  time  is  within 
the  time  span  of  the  interpolation  buffer. 

If  not,  go  to  E.  Otherwise,  interpolate  for 
the  elements  at  the  new  epoch  time.  Save 
other  necessary  values  for  PRERES  subroutine 
and  set  PRDSAV  »  0. 

E.  Search  for  elements  corresponding  to  observation 

time: 

(1)  If  the  observation  time  is  not  within  the 
time  span  of  the  interpolation  buffer,  go 
to  step  C  to  omit  the  first  time  and  add  a 
new  one . 

(2)  If  it  is  within  the  time  span,  interpolate 
for  the  elements  at  the  observation  time. 

Pass  these  elements  to  the  C0EFF2  subroutine 
which  will  compute  the  residuals  and  the 
partial  derivatives  and  enter  them  in  the 
least  squares  matrix  buffers. 

(3)  Output  the  residuals  if  requested. 

(4)  Retrieve  the  next  observation: 

If  on  the  same  side  of  epoch  as  the 
.previous  observation,  go  to  step  E; 
if  not,  go  to  step  F. 

If  end  of  observations,  go  to  step  G. 
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CNTR- 
SPTRDEC 
5  of  5 


DESCRIPTION: 

(continued) 


(a'  If  yes,  and  if  elements  have  not 
been  found,  save  the  buffers 
necessary  to  restart  integration 
by  PRERES  subroutine. 

(t ;  If  after  epoch,  reinitialize 

several  switches  and  go  to  step  B, 

G.  Test  for  new  epoch  elements: 

(1>  If  the  new  epoch  elements  have  been  found, 
go  to  step  H.  If  not,  determine  where  new 
epoch  is  in  relation  to  epoch.  Save  buffers 
to  restart  under  several  conditions.  At 
times  it  is  more  convenient  to  initialize 
for  the  integration  to  the  new  epoch  elements, 
i.e.  if  all  the  observations  are  before  epoch 
and  the  new  epoch  is  after  epoch. 

H.  Solve  for  the  delta  elements:  j 

(1)  Compote  the  root  mean  square  of  the  residuals. 

(2)  ft  there  art  enough  observations  accepted,  J 

solve  the  l^ast  squares  matrix  with  subroutine 
ISQS.  If  nor,  take  an  error  exit  to  location 

FINISH .  j 

t 

(3)  The  its  'Its  of  solving  the  least  squares  matrix 

(delta  elements)  art  in  a  b_ffer.  Unpack  the  . 

buff.tr  ard  exit. 


F.  Reinitialize  to  integrate  forward  from  epoch: 

(1)  If  already  initialized,  go  to  step  E. 

If  not,  test  if  new  epoch  elements  are 
before  epoch. 
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AK0-7101  FIGURE  6.  C.NTRL  FLOW  DIAGRAM  (2  of  2) 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


SUBROUTINES: 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


COE  EFT 
SPIRDEC 
1  of  2 


To  compote  the  coefficients  for  the  least  squares 
matrix  and  residuals  for  each  observed  quantity. 


JMP  COEFF2 

OBFLC 

= 

observed  quantitie 

SIGMA 1 

- 

SIGMAr 

- 

<?6 

STGMA3 

= 

SIGMA4 

= 

Tr 

RCSDL 

= 

( P  -  0C)  (KM) 

ARSDL 

= 

P  A.  (L  -L  )  (KM) 
c  —  —  — c 

DRSDL 

= 

Pc  D.  (L  -Lc)  (KM) 

•  • 

RRSDL 

— 

(D  -  0C)  (KM/SEC) 

RC-FLG 

-1 

L  if  accepted 

ALE  LG 

= 

*  if  rejected 

> 

L  OB' LG 

= 

RRF  LG 

- 

REJFLG 

=1 

*  if  any  rejected 

J 

A  if  all  accepted 

Modified  A  and  B  matrices  in  All  and  B1L  buffers. 

Program  -  SOBXYZ,  RHOSB,  COMXLX,  DOT PR ,  CMFlF,  RESREJ1, 
LSQR,  COMDEX-,  CORRD,  REJECT  1 ,  REJECT?. 

95  Ce  Lis 

First  computes  R's  and  L ' s  at  time  t. 

The"  tects  -  'Pr  LG  ’Or  observed  quantities  and  computes 
the  specified  residuals  and  coefficients  using  subroutines 
CMFCF  and  CORRD.  Tests  if  the  residuals  are  accepted  or 
rejected.  If  accepted, the  coefficients  are  added  to  the 
matrices  A  and  B  using  subroutine  LSQR.  If  rejected,  a 
flag  is  set,  corresponding  to  the  residual  to  indicate 
rejection. 
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C0EFF2 
SPIRDEC 
2  of  2 


DESCRIPTION 
(cont inued) 


4  switches  mjst  be  preset  before  calling  this  subroutine: 


for  1st 

pass  only 

COEFR 

JMP 

REJECT1 

COEFA 

JMF 

REJECT1 

COEFH 

JMP 

REJECT  1 

CCEFRR 

JMP 

REJECT2 

fcr  the 

2nd 

through 

COEFR 

JMP 

RESREJ 1 

COEFA 

JMP 

RESREJ1 

COEFH 

jMP 

RESREJ 1 

COEFRR 

JM? 

RESREJ 1 
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PURPOSE 

CALL  SEQUENCE 

INPUT: 


OUTPUT: 


SUBROUTINES; 

STORAGE 

REQUIREMENTS 

DESCRIPTION; 


COE  LIS 
S^IRDEC 


To  prepare  values  for  output  by  ELMOL'T 

:MP  COELTS 

XLO 
AXNO 
A'-'NO 
HXO 
r.  o 
HZO 
B 

FEV 
QO 
SI  MI 
ESQ 
P 

PTUOZ 

XNO 

LPRINT  =  L  (cleg) 

AXPRINT  =  a 

XII 

AVPRINT  = 

y  n 

HXPR1NT  -  li’ 

x 

HYPPINT  =  h 

V 

HZ  PR  IN  7  =  h 

z 

SPRINT  =  B 
PREV  =  Rev  No. 

HSUBQP  =  h  (k  ' 

q 

P APR! NT  =  (minutes) 


■None ) 


19  Cells 


Computes  HSUBQP  and  PAPRINT  -  moves 

'  sed  to  prepare  initial,  final,  and 
ou Lput  . 


othei.  values  to  output  cells, 
piediction  elements  for 
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COMDEL 

SPIRDEC 


PURPOSE: 


To  compute  AL  ,  AL  ,  AL 
x  y  z 


CALL  SEQUENCE: 

JMT  COMDEL 

INPUT: 

XLS  1 PX 

XLSUBY 

from  AZREC 

X--SUBZ 

XLX  = 

P  /  P 

XLU  = 

x  c 

P  /  P 

K 

M 
^  T 

X 

V  c 

P  /  P 

OUTPUT: 

DELTX 

z  c 

SUBROUTINES : 

DELTA 

DEulZ 

(Near) 

STORAGE 

REQUIREMENTS : 

6  Cells 

DESCRIPTION: 

DELTX  = 

XLS'BX  -  XLX 

DEL7Y  = 

XLSUBY  -  XLY 

DELTZ  = 

XLSUBZ  -  XLZ 
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PURPOSE:  To  compute  delta  q. 


CALL  SEQUENCE:  JMP  COMEQ 


INPUT:  DLTNN 

= 

An/n 

XNO 

= 

n 

0 

DLTAX 

= 

Aa 

X 

DLTAY 

= 

Aa 

y 

AXNO 

a 

xno 

AYNO 

= 

a 

yno 

QO 

qo 

EO 

= 

e 

o 

AO 

= 

a 

o 

OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


DELTAQ  =  Aq 
Philco  -  FSQRT 


13  Cells 


If  e  ^  0,  then 
o 

q  =  -  2/3  An/ri  qQ 


a 

o 


I  a  Aa  +  a 
|  xno  xn  yno 


e 

o 


If  e  =0,  then 

o 

q  =  -  2/3  An,  q  -  a“V/Aa2  +  Aa2 
^  /n  oy  xn  y 


COMDQ 

SPIRDEC 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


COMRMS 

SPIRDEC 


To  compute  root  mean  square 
IMP  COMRMS 

RCNT  =  C/'HLT,N  where  N  -  number  of  accepted 

weighted  residuals 

SUM  =  sum  of  squares  of  accepted  weighted  residuals 

ACCCNT  =  C/dLT,N 

RCNT  =  N  (floating  point) 

RMS  -  root  mean  square 

(None) 

16  Cells 


ACCCNT  =  RCNT 

RCNT  is  converted  to  F.P. 

RMS  =  VsUN/RCNT 
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COMXLX 

sfir:  fx 


PURPOSE: 

To  compete  L 

CALL  SEQUENCE 

UMP  COMXLX 

INPIT: 

RHOX  =  P 

X 

RHOY  =  0 

y 

RTOZ  =  p 

z 

RriOC  =  p 

c 

OUTPUT: 

XLX  =  L 

X 

XL!  =  L 

y 

XLZ  =  L 

z 

SUBROUTINES : 

'[None', 

STORAGE 

REQUIREMENTS : 

5  Cel  Is 

DESCRIPTION: 

XLX  -  p  /  p 
x  c 

XLY  =  p  /  0 

y  c 

XLZ  =  p  /  p 
z  c 

ill 


CORRD 
SPIRDEC 
1  of  2 


PURPOSE: 

To  compute  coefficients  of  range-rate  correction 
equations  and  the  range-rate  residuals. 

CALL  SEQUENCE: 

JMT  CORRD 

TAM  (unweighted  residual) 

TQM  (weighted  residual) 

INPUT: 

Ki'.TRI  =  0  OXXXXXXX 

LIEU  IS  =  C/HLTR  TERMS 

;nd  output  of  C0EFF2  subroutine 

iwitch  CORRDS4  must  be  preset  to  JMP  CORRDC  befor 

enti ring  this  routine. 

OUTPUT: 

Set  description. 

SUBROUTINES: 

Pr.grim  -  RDTSB,  CDLTBA 

STORAGE 

REQUIREMENTS: 

15  Cells 

DESCRIPTION: 

First  these  formulas  must  be  solved: 

EXOM  =  a  (e  cos  E  -  e^) 

EYOM  =V  1  e2  (a  e  cos  E) 

VDOT  =  rv/  r 
RSQ  =  r2 
RC'JBE  = 

RVDOT  =  rv 
XMUA32  fi3/2 

a  (e  cos  E  -  e2) 

m 

UDOTU  Vl-e2  (a  e  cos  E) 


RDOTU 


,3/2 
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DESCRIPTION: 

(continued) 


,  'RRI 
sf  ;rtfv. 


RDOrx  -  r  in  +  C>-  (RDOT_.i 


'JDOTN  =  r‘v/3  +  (U-7  )  (UDOTU) 


RDTXN  =  |  sin  +  cc)  -  a  -  (a  ;  (e  sir.  e) 

yr.  xn 


RI  T-'N  =■ 


-  c,3  (t,  +  t)  *•  ‘a  )  (e  sin  e; 

>.n  y  r. 


l/Ta5/- 

3 

r 

V^~a5,/z 

3 


HTXN  - 


.  VT7\Ta5/z 


I'f'N  = 


.VT-VTO'- 


cos  CE  +  i*-’-)  -  a  (  —  +  l'i 
xn  pa 


2 

sin  (E  +  a)  -  a  (  —  +  l'» 
yn  pa 


RPOTOR  =  i,  r 


TEMP 

TEMPI 


p  r  +  P  v  +  P  y 

xx  v  y  z  z 


=  P  v  +  P  v  +  P 


x  > 


z  z 


PW+PW+PW 
xx  y  y  z  z 


Faxing  »•  .mp-.ied  these  f.cmtlas,  the  following  quantities 
art.  s  red  in  tve  TERMS  bt.fi cr  as  described  by  the  bits 


ir  KTv'h--: 


i r. 

"Xxn 

A-  V: 
Abo 

■;A^ 

CAi 

CA;;./b 


?  Correction  coefficients 


R 


(weighted  residual) 


1U 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


CREKT 
SPIRDEC 
1  of  2 


To  apply  corrections  to  the  elements. 

JMP  CREKT 

DLTNN  =  An/n 
DLTAX  =  Aaxn 
DLTAY  =  A ayn 
DLTND  =  A 
DLTIN  =  a i 
DLTUO  =  AU 

o 

DLTB  =  AB/b 
See  description 

Philco  -  FSIN,  FCOS,  FSQRT,  FL0G2X,  F2X 

46  Cells 

The  delta  elements  are  applied  and  new  elements 
are  computed  as  follows: 


XNO 

= 

n  =  n  (1  +  An/n) 

o  o 

B 

= 

C ,  A/m  =  B  (1  +  AB/b) 
d 

UO 

= 

U  =  U  +  AU 

OO  0 

AXNO 

= 

axno  &xno  + 

Aa 

xno 

AYNO 

= 

ayno  3yno  + 

Aa 

yno 

ESQ 

= 

2  2 

e  o  -  a  + 

xno 

2 

a 

yno 

EO 

= 

e  =vrr 

0  o 

XNODEO 

= 

^0  =  ^0  ^^0 

SINO 

= 

sin 

COSO 

= 

cos 
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CREKI 
S FIRE EC 
2  of  2 


DESCRIPTION: 
(continued ) 


XT.WJL 

— 

1  -  i 

+  M 

COST 

= 

cos  i 

SINI 

= 

sin  i 

wX 

= 

(sin 

i)  (sin  Cl) 

WY 

= 

-  (sin 

i)  (cos  Cl) 

WZ 

= 

cos  i 

XLO 

= 

L  =  1 
o 

J  *4 
0 

-  fi 

AO 

= 

<x  — 

0 

[0 

2/3 

in  , 

\  cl 

P 

= 

p  =  a. 

o  (1  "  eo2) 
0  0 

RTP 

= 

F 

HXO 

= 

h 

xo 

=  wx 

H70 

= 

h 

yo 

•\F  uv 

HZO 

= 

h 

zo 

AX  9 

= 

a  : 
xo 

=  (cos  0) 

(a  )  - 
xnc 

(cos  i) 

(sin  Cl) 

(a 

yno 

rtT  0 

- 

a 

yo 

=  / 

[sin  Q) 

(a  )  + 
xnc 

(cos  i) 

(cos  0) 

(a 

yno 

AZO 

= 

a  = 

zo 

=  (sin  i) 

Sno’ 

QO 

= 

* 

a 

c 

(1  - 

eo) 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


SUBOUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


DELOUT1 

SP1RDEC 


To  print  the  delta 

JMP  DELOUT1 

DLTNN  =  An/n 
DLTAX  =  A 
DLTAY  =  A 
DLTUO  =  A“q  ) 

DLTNO  =  £$'' 

DLTIN  =  Ai 
DLTB  =  AB/B  * 

RMS  =  weighted 

DELTA  N/N 
DELTA  AXN 
DELTA  AYN 
DELTA  UO 
DELTA  NODE 
DELTA  I 
DELTA  B 
RMS 

S\stem  -  GLOP,  PANT 
Program  -  PAGECON 

55  Cell* 

Forces  a  page,  prints  headings,  values,  10  spaces,  and 
updates  PAGECON. 


elements  after  each  D.  C.  pass. 


all  floating  point 


root  mean  square 
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DFLO’  12 
S.FIRDEC 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


SUBROUTINE: 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


To  print  the  last  pass  element  corrections, 
old  and  new  rrr.s ,  and  accepted  and  rejected  count. 

JMP  DELOUT2 

DLTNN 

DLTAX 

DLTAY 

DLIJO 

DLTND 

DLTIN 

DLTB 

RMS 

ORMS 

ACCCNT 

RE.JCNT 

DELTA  N/N 
DELTA  NXN 
DELTA  AYN 
DELTA  UO 
DELTA  NODE 
DELTA  I 
DELTA  B 
OLD  RMS 
NEW  RMS 

No.  of  residuals  used  and  rejected 

System  -  PANT,  GLOP 
Program  -  PAGECON 


56  Cells 
Forces  a  page. 

Prints  the  comment  "DC  Converged  ••  The  next  corrections 
would  be...."  Then  prints  headings,  delta  elements,  old 
and  new  rms,  accepted  and  rejected  residual  count,  gives 
8  spaces  and  updates  PAGECON. 
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FURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


DELTAU 

SPIRDEC 


To  compute  Au 

JMP  DELTAU  ‘ 

XOBS 

YOBS 

ZOBS 

UX 

UY 

UZ 

DELU 


Program  -  ARC TAN 
19  Cells 

RU  =  XOBS • UX  +  YOBS-UX  +  ZOBS*UZ  =  COSDU 

SINDU  -  WX  (UY  ZOBS  -  UZ’YOBS)  +  WY  (UZ-XOBS  -  UX1 ZOBS) 
+  WZ  (UX-YOBS  -  UY-XOBS) 

DELU  =  tan  "1  (SINDU/RU)  (mod  2tt) 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 


SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


DERIV 

SPIRDEC 


To  commif-p  f-he  H  ornraH'vpe  « C  w  w  _ 

,  -  -  j. uuc  ri  exements 

with  respect  to  time,  using  the  effects  of  radiation 

pressure,  drag  and/or  bulge  perturbation  as  specified 
on  I  Girds. 

JMP  DERIV 

See  Input  for  MARTINI,  JNDRAG,  RDPRES 
3rd  SUBXYZ  subroutines. 


XDGR 

- 

• 

X 

\ 

W  +  187 

= 

d(L)/dt 

YDGR 

= 

. ' 

y 

188 

= 

d  (a^)  /’dt 

ZDGR 

z 

\ 

189 

- 

d (a  ) /d  t 

y 

DGR 

~ 

D 

190 

= 

d(a  )/dt 

Z' 

DDGR 

= 

D 

191 

= 

d(hx)/dt 

D 

= 

D 

192 

- 

d(hy)/dt 

AXGR 

AZGR 

= 

\ 

a 

\X 

a 

193 

- 

d(h  )/dt 
z 

SMLGR  =  £ 

Program  -  SUBXYZ,  MARTINI,  JNDRAG,  RDPRES 


97  Cells 

-ses  SUBXYZ  to  compute  r,  r;  computes  specified 
perturbative  effects  to  be  used;  and  computes  the 
derivatives,  storing  them  in  the  ADRASH  buffer 
t.W  +  187  to  W  +  193)  . 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 


SUBROUTINE: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


D IV  DIF 
SPIRDEC 
1  of  2 


10  interpolate  for  elements  using  a  fourth-order 
divided  difference  method. 

JMP  DIVBIF 

ICK  buffer  (See  description) 

L  =  £  (min.  since  epoch) 

TBJF  =  temporary  storage  buffer 


ICK  +  40  =  L 


41 


42 

43 

44 

45 

46 


-  a 


=  a 


=  a 


=  h 

x 

=  h 

y 

=  h 

z 


>  at  time  t 


None 

38  Cells 

G_ven  the  5  sets  of  time  and  elements  in  buffer  ICK 
*  .  subr^tine  will  interpolate  for  the  elements  at’ 
„ime  t,  where  tQ  <  t  <  t  4,  using  the  following 

formula  to  interpolate  for  each  element: 


x  =  x  +  (t^t  ) 

O  C/ 


x.-x 
1  o 

trto 


+  (t-t^Ct-tp 


/x2-xil 

«• 

lxrx  o’ 

IvhJ 

i  trto 

t9-t 
2.  o 


+  (t-to)(t-tl)(t-t2) 


x3-x2 

.  vV 

/x,„-X  X  -X 

2  1  1  O 

«• 

rt 

CO 

1 

rt 

ro 

t2"tl 

t2~tl 

fc3 

-  C1 

h-t*  /, 

t0  -  t 
3  o 
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DIvDIF 
SlIRDEC 
2  cf  2 


DESCRIPTION: 

(continued) 


+ 


(t-tQ)  (t-tp  (t-t?)  (t-t3> 


fx  _x 

4  3 

X3-5V 

' x3"x2 

Vx0 

'4  "3 

t3mt2 

4-  4. 

--  -r’l 

rt 

-  t. 

t„ 

-  t. 

V  4 

2  J 

l  3 

w 

':4‘tl') 

w 

\ 


X3"X2 


X2“X1 


t3“t2 


hmti 


.  V 


t3"tI 


X2~X1 


X-  -x 
1  o 

t.-t 
1  o 


t0-t 
2  & 


ICK  Buffer  Format 


ICK  +  0 
+  1 
+  2 
+  3 
+  4 
+  5 
+  6 
+  7 


t  32 
+  33 
+  34 
+  35 
+  36 
+  37 
+  38 
+  39 


L 

o 


a 

x 


a 

y 

a 

z 

h 

x 

h 

y 

h 

z 


> 


5  sets  of  time  and  elements 
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PJTODncir . 

^  U*-VA  K/UU  • 

CALL  SEQUENCE: 

INPUT: 


OUTPUT: 

SUBROUTINES: 

DESCRIPTION: 


m 


DOTPR 

SPIRDEC 


To  compute  dot  products  of  vector  in  the  (A),  (q; ,  (D) 
registers  with  U,  V,  W 

TMA  (L  ) 
x 

TMQ  (L  > 

y 

TMD  (L  ) 

7. 

JMP  DOTPR 

UX 

UY 

UZ 

VX 

VY 

VZ 

WX 

WY 

WZ 

ADOTU 

ADOTV 

ADOTW 


(None) 


ADOTU 

=  I 

U 

+ 

L 

u 

+ 

L 

U 

X 

X 

y 

y 

z 

z 

ADOTV 

*  L 

V 

+ 

i 

V 

+ 

L 

V 

X 

X 

y 

y 

z 

z 

ADOTW 

=  L 

W 

+ 

L 

w 

+ 

L 

W 

X 

X 

y 

y 

z 

z 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 
SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


ELMO.  I 
SPIRDEC 


To  output  initial,  final,  and  new  epoch  element?  at 
the  end  of  the  differential  correction. 

JMP  ELMOUT 

Initial  elements  in  INELT  buffer. 

New  epoch  elements  in  PREDBF  buffer 
Final  elements  in: 

XLO 

AXNO 

AYNO 

HXO 

HYO 

HZO 

B 

REV 

See  description 

Program  -  RINEL,  BCDTIM,  COELTS ,  PRERES,  PAGECON 
System  -  PANT,  GLOP 


157  Cells 
Prints  comment: 

"Initial,  Final,  and  New  Epoch  Elements"  Prints  headings, 
then  restores  initial  elements  to  output  cells  and  prints 
them.  Prepares  and  prints  final  elements.  Uses  sub¬ 
routine  PRERES  to  obtain  new  epoch  elements  and  prints 
them.  Then  prints  comment:  "End  of  DC"  and  updates 
PAGECON.  Also  prints  card  images  of  P  Cards  necessary  to 
run  CALIB . 

Punches  the  following  cards  in  P  Card  formats: 

(1)  Corrected  elements  (PI  &  P2) 

(2)  New  Epoch  elements  (PI  &  P2) 

(3)  Cards  for  calibration  program  (PI,  P2 ,  P3) 
and  the  cards  P5-P8  which  were  used  in  the 
current  SPIRDEC  run. 
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ELMOLIi 
S FIR DEC 
I  of  2 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


To  prirt  the  corrected  elements  after  each 
differential  correction. 


JMP  ELMQUT1 


PREV 
LPRINT 
AX PRINT  = 
AYPRINT  = 
HX PRINT  = 
HYPRINT  = 
HZPRINT  = 
BPRINT  = 
HSUBQP  = 
PAPRINT  = 

OBYEAR 

OBM3 

OBDAY 

FOBHR  > 
FOBHR 
FOBMIN 
FOBSEC 


rev.  no. 

I-  (deg) 
a 

xno 

a 

yno 

h 

x 

h 

y 

h 

z 

C,  A/m 
d 

Hq  (km) 
Period  (min) 


BCD  Time 


REV  ^ 

i 

i 

TIME  I 


L 

AXN 


AYN  ! 

HXO  > 

HYO 

HZO 

B 

PER.  ALT. 

PA 


These  are  not  program  locations 
bet  output  headings 
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SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


ELMOUT1 
SPIRDEC 
2  of  2 


Program  -  COELTS,  PAGECON,  BCDTIM,  PRNTMAT 
System  -  PANT,  GLOP 


53  Cells 

Uses  BCDTIM  to  prepare  time  for  output.  Uses  COELTS 
to  compute  values  and  correct  elements  to  the  output 
units.  Prints  "Corrected  Elements".  Then  prints 
headings,  values,  and  updates  PAGECON.  If  the  DC  is 
converging,  PRNTMAT  will  be  called  to  print  the 
standard  deviations  and  correlation  matrix  of  the 
delta  elements.  If  divergent,  the  last  corrected 
element  set  will  be  punched  on  P  cards  1  and  2. 
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8 


FTAFEW 

SPIRDEC 


PURPOSE: 

To  wrapup  the  binary  ephemeris  .  ■ 

CALL  SEQUENCE: 

JMP  FTAPEW 

INPUT: 

TAPBUF  =  ephemeris  buffer 

TAPCNT  =  C/HLT,  TAPBUF  +  X 

OUTPUT: 

Sentinal  of  Z's 

SUBROUTINES : 

System  -  S\S ,  SYSNO,  SYSIO 

STORAGE 

REQUIREMENTS: 

9  Cells 

DESCRIPTION: 

Writes  the  final  block  on  the  binary  ephemeris 
tape  on  logical  10  -  the  sentinel  being  Z's  in 
the  current  blo:k.  Then  the  tape  is  rewound 

with  lrckout. 
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I 

I 

|  PURPOSE: 

CALL  SEQUENCE: 

I 

I  INPUT: 

OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


GETSEN 
SP  IRT  EC 


To  retrieve  sensor  information. 

JMP  GETSEN 
(No  sensor  data) 

(Normal  return) 

BIASAD  =  C/HLT,  EBLOC ,  C/HLT,  BIBUF 

BIBUF  =  SBUF  =  Modified  buffer  of  sensor  information 
STAID  =  00000SSS 

PHIRD  =  Cp 

XLAMBA  =  X 

XOVCT  =  x  /  cos  9 
CAPZ  =  Z 

None 

10%  Cells 

Given  a  sensor  number  in  STAID,  the  routine  will 
search  through  SBUF  until  a  match  or  Z's  are  found. 

(1)  If  Z's  are  found,  the  error  exit  (+1H)  is  taken. 

(2)  If  a  match  is  found,  the  information  will  be 
unpacked  to  the  cells  listed  under  output. 
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PURPOSE: 

CALL  SEQUENCE 

INPUT: 

OUTPUT: 


SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


GEIWGI 

SPIRDEC 


:  JMP  GETWGr 

Weights) 

(Normal  Return) 

Weights  in  WBUF 

staid1  :  ooooJssfI;  C/HLT’  BIBUF 

00000SSS  Where  SSS  is  the  station  number 
SIGMA  1  =  op 
SIGMA2  =  <7p 
SIGMA3  =  cta 
SIGMA4  = 

None 


11  Cells 

£  ?  SLfs:  ls 

and  6x11  +2h-  « ~  -tsch 
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GIPAR 

SPIRDEC 


PURPOSE:  To  prepare  input  for  the  GIPAR  program. 

CALL  SEQUENCE:  JMP  GIFAR 


INPUT:  T 

= 

minutes  since  epoch 

ORGDAY 

s 

dajc  since  beginning  of  year 

ORGTM 

= 

fraction  of  epoch  day 

X 

= 

x  (in  km  or  e.r.) 

Y 

= 

y  (in  km  or  e.r.) 

Z 

= 

z  (in  km  or  e.r.) 

FTFLAG 

PFLAG 

GIPADR 

s 

left  address  is  the  next  location 
in  the  GIPAR  buffer 

OUTPUT: 


SUBROUTINES: 


An  entry  in  the  GIPAR  buffer: 


Word  0 

1 

2 

3 

4 

5 


x 

y 

z 

day  number 
fraction  of  day 
ZZZZZZZZ 


Program  -  SEPSUB 


STORAGE 

REQUIREMENTS:  21  Cells 


DESCRIPTION:  Before  the  first  entry  to  this  subroutine  GIPADR 

must  be  set  to  C/HLT,  EBLOC  +  20  and  EBLOC  +  20 
must  be  set  =  ZZZZZZZZ.  The  subroutine  will 
update  GIPADR  as  it  makes  entries  to  the  buffer. 

(1)  Convert  time  from  minutes  since  epoch  to  days 
since  the  beginning  of  the  year  and  fraction 
of  day;  then  store. 

(2)  Test  PFLAG: 

(a)  If  PFLAG  =  0,  store  x,  y,  z 

(b)  If  PFLAG  f4  0,  convert  x,  y,  z  from  km.  to  e.r. 
and  store.  PFLAG  t  0  means  subroutine  SUBOUT 
has  converted  x,  y,  z  to  km. 

(3)  Store  sentinel  of  Z's  and  update  GIPADR. 
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HEAD 

SPIRDEC 


PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


To  output  a  new  page  with  the  page  headings . 

JMP  HEAD 

Parameters  set  by  1HEAD1  subroutine. 

New  page  with  the  described  page  headings. 

System  -  PANT,  GLOP 

52  Cells 

1st  line  -  Spiral  Decay  SPDC  Program  Page  X 

2nd  line  -  Satellite  No.  =  XXX  Satellite  name  -  x-x(10  charact 

2nd  line  (cont)  -  Element  set  no.  =  xxx  Time  of  EPOCH  = 

YY  MM  DD  HH  MM  SS  .SSS 
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I HE ADI 
SPIRDEC 


PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


To  initialize  the  page  heading  routine  (HEAD) 


JMP  I  HE  ADI 
(Return) 


SATN  =  0  - 0  NNN,  Satellite  number 

ORGDA  =  days  since  19501  £ 

ORGTM  =  fraction  of  dayj 


PAGENO  =  0 — 0 

PDAY  =  MMO  —  0 

PYEAR  =  ,  A  19YY00 

SATEL  =  NNNO  —  0,  N  =  SATN 

HEAD 2  +  3  H  =  ADDR  of  BCD  MONTH 

HEAD 5  to  HEAD 5  +  5  =  (BCD  EPOCH  Time  in  cells  for 

output  call  sequence.) 


System  -  AKLOK 
Program  -  BCDTIM 


26  Cells 

Uses  AKLOK  to  get  PDAY  and  PYEAR  - 
Uses  BCDTIM  to  get  Epoch  in  BCD  ready  for  output 
Sets  up  parameters  necessary  to  output  page  heading 
as  done  by  HEAD. 


•Nil :al 

SP  \'f<DE( 


PURPOSE: 

T-o  initialize  for  the  ADBASH  s  ib  routine. 

CALL  SEQUENCE* 

IMP  INITIAL 

INPUT: 

XLO 

AXO 

AfO 

AZO 

HXO 

diO 

HZC 

OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


See  description. 
Pr  ogram  WSETUP 


9  Cells 

Moves  elements  to  W  buffer  by  calling  WSETIT 
Sets  I  (W  +  3)  =  F/0 
Sets  switch  ADP  to  1/1T15 

Sets  switches  ADDER 1  to  ADDER 7  to  a  (JMP  DERI/; 
in  the  ADBASH  subroutine. 


132 


iftlTL 

sfclRDEC 


PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 
SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


To  print  the  cards  images  of  the  parameter  cards 
(P  cards  1  thru  10)  and  process  them. 

JMP  INITL 
./MP  (ERROR) 

(NORMAL) 

Parameter  cards  (P  cards)  in  CONBUF 

See  description. 

Program  -  PAGECON,  INPUT 
System  -  PANT 


19  Cells 

Prints  card  images  of  cards  in  CONBUF;  then  calls 
INPUT  to  process  these  cards. 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 
SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


INPUT 
SPIRDEC 
1  of  10 


To  unpack  and  validate  the  input  parameter 
cards  (1-10). 

JMP  INPUT 
JMP  (error) 

(Normal  return) 

P  Cards  (1-10)  in  CONBUF 
PCOUNT  «  C/HLT,  N 

where  N  *  number  of  P  cards 
See  description. 

Program  -  XSRCH,  RR2AHL,  CARDER,  SAVCON 
System  -  FSRLOK,  FXFLT,  INITEL,  NXTELM,  FYKLOK 


370  Cells 


The  input  parameter  cards  need  not  be  in  order 
unpacked  according  to  the  number  in  Column  79 
subroutine  XSRCH.  ’ 


they  are 
using 


If  no  error  is  found  on  a  card,  a  bit  will  be  set  in 
location  CARDS: 


1/1T0  - 

P  Card 

1 

1/1T1 

P  Card 

2 

1/1T2 

P  Card 

3 

1/  1T3  - 

P  Card 

4 

1/1T4 

P  Card 

5 

1/1T6 

P  Card 

6 

1/  IT  7 

P  Card 

7 

1/1T8 

P  Card 

8 

1/1T9  - 

P  Card 

9 

1/1T10  - 

P  Card 

10 
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DESCRIPTION: 

(continued) 


All 

floating 

point 


INPUT 
SPIRDEC 
2  of  10 


The  remainder  of  any  card  in  error  is  not  validated 
beyond  the  error;  however,  an  attempt  will  be  made  to 
validate  the  remaining  cards.  Subroutine  CARDER  will 
be  used  to  output  the  card  image  of  the  erroneous  card 
with  a  comment  indicating  the  first  field  in  error; 
location  CARDSW  will  be  set  ^  1.  If  an  error  occurs 
during  XSRCH  conversion,  a  comment  will  be  printed  with 
the  card  image,  but  with  no  field  indicated. 


The  subroutine  XSRCH  will  unpack  each  of  the  cards  into 
the  locations: 


P  Card  1 

If  Col.  78  =  A 

If  Col.  78  =  B 

^  Col .  1-12 

XLO  (rad) 

X  (km) 

13  -  24 

AXNO 

Y  (km) 

l 

25  -  36 

AYNO 

Z  (km) 

37  -  49 

HXO  (er/kmnin) 

XDOT  (m/sec) 

50  -  62 

HYO  (er/kemin) 

YDOT  (m/ sec) 

^  63-75 

HZO  (er/kemin) 

ZDOT  (m/ sec) 

If 

type  B  elements  are  input,  subroutine  RR2AHL  will  be 

used  to  convert  them 

to  type  A  elements 

with  the  results 

stored  in  the  locations  for  type  A  elements.  If  no  errors 

set  bit  0  in  CARDS  = 

1 

P  Card  2 

Col.  1-5 

SATN 

(Binary) 

6-15 

SATNM,  SATNM+1 

(Binary) 

17  -  32 

YY,  YY  +  1 

(Binary) 

33  -  37 

EPREV 

(F.P.) 

38  -  40 

ELNO 

(FX.P) 

Having  unpacked  the  card,  the  following  will  be  done: 

TOY  =  0 - 0  YY,  YY  =  BCD  year 

Using  subroutine  FSKLOK; 

ORDGA  =  days  since  1950  to  epoch  (F.P,) 

ORGTM  =  fraction  of  epoch  day  (F.P.) 

Then  the  card  image  will  be  saved  using  subroutine  SAVCON, 
for  use  by  subroutine  ELMOUT. 

If  no  errors,  set  bit  1  in  CARDS  =  1 
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INPUT 
SPIRDEC 
3  of  10 


DESCRIPTION:  P  Card  3 

(continued) 

Col.  1 
2 
3 
6 

7 

8 
9 

10-25 

26-32 

33 

34 

35 

36-39 

40-42 

43-45 

46-48 

49-51 

52-58 

59 

60 

62-64 

66-68 

69-71 

72-75 

76-78 


CRAIG 

(FX.P.) 

WGTFLG 

(FX.P.) 

PW 

(FX.P.) 

BFLAG 

(FX.P.) 

DFLAG 

(FX.P.) 

RPFLAG 

(FX.P.) 

PREDFLG 

(FX.P.) 

PRTIM,  PRTIM+1 

(BCD) 

SKNTRL 

(BCD) 

RPT 

(FX.P.) 

CNF  LAG 

(FX.P.) 

DQFLAG 

(FX.P.) 

DQN 

(F.P.) 

ABSMX 

(F.P.) 

ABMX2 

(F.P.) 

XISTSG 

(F.P.) 

KAPPA 

(F.P.) 

B 

(F.P.) 

R  OF  LAG 

(FX.P.) 

RESOPT 

(FX.P.) 

CONTEST 

(F.P.) 

F10 

(F.P.) 

F10AV 

(F.P.) 

AP 

(F.P.) 

GAMMA 

(F.P.) 
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DESCRIPTION 

(continued) 


INPUT 

SPIRDEC 

4  of  10 


Having  unpacked  P  Card  3,  the  following  validations 
and  conversions  will  be  made: 

If  R OF LAG  i  0,1,2,  set  R OF LAG  =  0 

If  RPT  =  0,  then  ERROR  (No.  of  iterations  not  specified) 

If  SKNTRL  =0,  then  ERROR  (Elements  to  correct  not  specified) 
If  B  is  to  be  corrected: 

1)  If  B  =  0,  then  ERROR 
If  B  is  not  to  be  corrected: 

2)  If  DFLAG  t  0,  then  B,  F10,  F1QAV,  AP  cannot  =  0 

3)  If  DFLAG  =  0,  then  B,  F10,  F10AV,  AP  are  not  checked 
If  DQFLAG  ^  0,  then  DQN  cannot  =  0 

If  ABSMX  =  o)  then  ERROR  (Rejection  criteria  cannot  =  0) 
and/or  ABMX2  =  Oj 

If  XISTSG  =  0,  then  set  =  F/1.5 
If  KAPPA  =  0,  then  set  =  F/l 

If  RPFLAG  ^  0,  then  GAMMA  ^  0,  and  compute  RPC0N3  =  GAMMA.  B/2.2 
If  DFLAG  =  i ,  set  RHO  =  F/0 
If  PREDFLG : 

1)  =  0,  use  FSKLOK  to  convert  PRTIM,  PRTIM+1  to  days 

since  1950  and  fraction  of  day  respectively. 

2)  =1,  use  FXFLT  to  convert  PRTIM,  PRTIMHT  to  revolution 

number  and  store  in  PRTIM 

3)  =  2,  do  not  use  PRTIM,  PRTIMH 
If  CONTEST  =  F/0,  set  =  F/,01 
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DESCRIPTION: 

(continued) 


INPUT 
SPTRDEC 
5  of  10 


Convert  S KNURL  f  r  mi  BCD  to  binary  in  reverse  order: 

SKNTRL  *  OXXXXXXX  (N,  AXN,  AYN,  U,  Q,  i,  B.  BCD 
to  SKNURL  =  0 - OXXXXXXX  (B,  i  ,  0,  U,  AYN,  AXN,  N)  binary 

Convert:  DQN  from  km  to  e.r. 

ABSMX  from  km  to  e.r. 

AEMX2  from  km/sec  to  er/kemin 

If  no  errors  are  found,  set  DCFLAG  ^  0  and  bit  2  =  1  in  CARDS, 
ose  SAV'ON  to  same  card  image  for  ELMOUT. 

P  Card  4 


2-5 

LTBUF+1 

(F.P.) 

6 

LTBUF+2 

(ex  .  P . ) 

7 

LTBUF+3 

(FX.P.) 

8 

LTBUF+4 

(EX.P.) 

9 

FTFLAG 

(FX.P.) 

1C-25 

PRTIME,  PRTIME+1 

(BCD) 

C7 

PFLAG 

(FX.P.) 

41 

3TFLAG 

(FX.P.) 

42-45 

LT3UF+5 

(F.P.) 

49-51 

PKAPPA 

(F.P.) 

52-58 

PB 

(F.P.) 

66-68 

PP  10 

(F.P.) 

'£> 

1 

i 

PF10AV 

(F.P.) 

VU-Tj 

PAP 

(F.P.) 

76-78 

PGAMMA 

(F.P.) 

138 


INF'JT 

SFIRDEC 

6  of  10 


DESCRIPTION:  The  following  validations  and  conversions  will  be  made: 

(continued) 

If  PKAPPA  =  F/0,  set  =  F/l 

If  DFLAG  t  0,  then  PB,  PF10,  PF1QAV,  PAP  cannot  =  0  or  ERROR 
If  PFLAG  =  0  and  BTFLAG  =  0,  then  ERROR  (no  output  requested) 
If  BTFLAG  =  1,  then  LTBUF  +  5  cannot  ^  0,  or  ERROR 
If  FTFLAG: 

1)  =  0,  use  FSKLOK  to  convert  (BCD)  PRTIME ,  PRIIME+1 

respectively  to  days  since  1950  and  fraction  of  day. 

2)  =  1,  use  FXFLT  to  convert  PRTIME,  PRTIME+1  to 

revolution  number  and  store  in  SREVF. 

If  LTBUF+4  ^  0,  then  PGAMMA  cannot  =  0,  or  ERROR 

If  no  errors  were  found,  set  PPPFLAG  f  0  and  bit  3  in 
CARDS  =  1 

P  Card  5 


/ 

Col.  1-10 

W  + 

12 

11-20 

+ 

13 

All 

floating 

21-30 

+ 

14 

point  l 

31-40 

+ 

15 

41-50 

+ 

16 

51-60 

+ 

17 

61-70 

+ 

18 

The  buffer  W  is  used  by  the  ADBASH  subroutine.  The  val.es 
input  will  override  the  assembled  values  in  these  locations 
as  the  absolute  error  criteria. 

SAVCON  will  be  used  to  save  the  card  for  ELM0UT.  If  no 
errors,  bit  4  in  CARDS  will  be  set  =  1. 
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IN?.  I 
SfIRDEC 
7  of  10 


DESCRIPTION: 

(continued) 


P  Card  6 


10-18 

13-/7 

28-36 


JBLTF 

JBUF+1 

JBUF+2 

JBUF+3 


All 

floating 

point 


If  this  card  is  input,  it  will  override  the  assembled 
va.^es  in  JB.P.  Subroutine  SAVCON  will  be  used  to  save 
the  card  image  for  ELMOLTT  and  bit  6  in  CARDS  will  be 
set  =  1  if  no  errors. 


P  Card  9 

Col.  1-5 
6-8 
9-13 
14-18 
19-21 
22-26 
27-31 
32-34 
3  5-39 
40-44 
45-47 
48-5/ 
53-5 } 
58-50 
6 1-65 
66-70 
"'1- 73 
74-78 


PRT+0 

1 

2 

3 

4 

5 

6 
7 


10 

11 

12 

13 

14 

15 

16 
17 


All 

floating 

point 
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INPUT 
SPIRDEC 
8  of  10 


DESCRIPTION: 

(continued) 


This  card  contains  the  intervals 

desired  to  be  printed 

or  put  on  an  empheris  tape  in  a  prediction. 

One  pass  is  made 

through  the  buffer  to  count  the  number 

of  sets  of  intervals  (3  fields  = 

a  set)  . 

Then  the  buffer  is  checked  to  set  that  the  times  are  in 

order.  If  not,  an  error  exit  is 

taken . 

If  the  times  are 

in  order,  PRT  + 

18  will  contain  the 

number  of  sets  of 
will  be  set  =  1. 

intervals  (T47) 

and  bit  9  in  CARDS 

P  Card  10 

Col.  1-15 

GIT  +  1,2  ' 

16-30 

+  3,4 

31-45 

5,6 

>  All  BCD 

46-60 

7,8 

61-75 

9,10  ‘ 

Each  set  of  2  words  is  checked  for  blanks.  If  not  blanks 
use  FSKLOK  to  convert  to  floating  point  days  and  fraction 

If  all  fields  are 

blank,  go  to  error  exit. 

Check  to  see  that 

times  are  in  order.  If  not,  take  error 

exit;  otherwise, 

set  bit  10  in  CARDS  =  1. 

P  Card  7 

Col.  1-9 

CSNM  +0  \ 

10-18 

+  1 

19-27 

+  2 

28-36 

+  3 

All 

^  floating 

37-45 

+  4 

poir 

46-54 

+  5 

55-63 

+  6 

64-72 

+  7  ) 
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INPJT 
SPIRDEC 
9  of  10 


DESCRIPTION: 

P  Card  8 

(continued) 

Col. 

1-9 

CSNM  + 

8 

\ 

10-18 

+ 

9 

19-27 

+ 

10 

28-36 

+ 

11 

>  All 

37-45 

+ 

12 

floating 

46-54 

+ 

13 

point 

55-63 

+ 

14 

64-72 

+ 

15 

/ 

These  cards  contain  the  tesseral  coefficients  for  the 
MARTINI  subroutine.  SAVCON  is  used  to  save  these  card 
images  for  ELMOUT.  If  either  or  both  cards  are  input, 
TESRAL  is  set  *  1.* 

If  no  errors  on  P  Card  7  &  8  bits  7  and  8  respectively 
will  be  set  =  l  in  CARDS. 

Having  unpacked  and  validated  the  input  parameter  cards, 
further  checks  are  made: 

If  neither  P3  nor  P4  was  input,  ie,  DCFLAG  &  PPPFLAG,  =  0, 
the  job  cannot  be  run. 

If  neither  or  only  one  of  PI  and  P2  were  input,  then  the 
elements  will  be  retrieved  from  EBLOC  and  a  comment 
printed.  If  EBLOC  has  no  elements,  the  job  will  be 
terminated . 

If  PREDFLG  =  0, 

PRTIM  will  be  made  minutes  since  epoch  by  a  double 
precision  operation  of: 

PRTIM  =  (PRTIM- ORGDA)  X  1440  +  (PRTIMtl  -  ORGTM) 

If  FTFLAG  =  0, 

PRTIME  will  be  made  minutes  since  epoch  in  the  same 
manner  as  PRTIM. 
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ISPTJT 
SPIRDEC 
10  of  10 


DESCRIPTION: 

(continued) 


bsing  FYKLOK,  ORGDAY  will  be  computed  as  days  from 
beginning  of  year  to  epoch. 

If  EPREV  >  100,000  it  will  be  modulated  100,000  as 
a  5  digit  revolution  number  is  maximum. 

If  TESRAL  ^  0,  ie,  P7  and  or  P8  were  input,  location 
CARDS  will  be  tested  to  see  that  both  were  input.  If 
not,  a  comment  will  be  printed  and  the  job  terminated, 
as  both  or  neither  are  required  for  the  MARTINI  subroutine. 

Having  successfully  passed  all  preceding  checks,  location 
CARDSW  will  be  tested.  If  =  0,  no  errors  were  found  on 
any  card,  and  the  job  may  be  run.  (EXIT  +  2H) .  If 
CARDSW  ^  0,  one  or  more  cards  were  in  error,  and  the  job 
will  be  terminated  (EXIT  4-  1H)  . 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


ITAPEW 

SPIRDEC 


To  initialize  ephemeris  tape  for  output. 


JMP  ITAPEW 


SATNM 
SATNM  +  I 


Satellite  Name  BCD 


I  block  on  tape  If 
See  description. 

TAPCNT  =  C/HLT,  TAPBUF 


Program  -  TAPEW 
System  -  SYS,  SYSNO 


8  Cells 

Positions  tape  to  1st  block  writes  1  block  on  tape 
(1st  2  words  are  the  satellite  name  -  last  2  words  are 
zero),  which  fulfills  requirements  of  the  XYZLA 
program  . 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 


OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS: 


JNDRAG 
SPTRDEC 
1  of  3 


To  compute  perturbative  effects  due  to 
the  atmosphere. 


JMP  JNDRAG 

JMP  (Error),  h  <  LOLIMIT 
(Nonna  1  Return) 


DFLAG 


LOLIMIT  = 


TEMPO 


:{ 
{ 


B 

KAPPA 

ASCON2 

TriDOT 


0  do  not  compute  drag  perturbations 
1  compute  drag  perturbations 

50  km  for  D.C.  and  Prediction 
10  km  only  for  Cowell  option 

Exospheric  temperature  at  epoch  from  TEMP 
subroutine 

from  P  Card  3  or  P  Card  4 
CD  AU 

from  P  Card  3  or  P  Card  4 
F/6378.165 
0 


Output  from  S'UBXYZ  routine  at  time  t. 
XDTGR  = 

YDTGR  =  *C 
ZDTGR  =  z' 


Phil/,  - 

Program  ■* 


FASIM,  F SIN,  FSI.V,  FSQRT,  FTENX 
’.  AuH,  ANGS.JN 


337  Cells 
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JNDRAG 
SPIRDEC 
2  of  3 


DESCRIPTION: 


atmosphere  tables from 120  *  ,  (1964>  Sarnie 

Ccesa  (1962)  tables  t  0°  “• 

at  location  NICOLET  is  comprised  of2°  a  !Hle  table 
words  for  10  fas.  altitude  Increments  °C“1 

subroutines  ^however3  ^there^ '  “> 

(JNDRAG I)  and  exit  (JNSWW  \8Pecial  entrance 

r“?v~ • S“"‘ “ 

height.  6r  t0  coinPute  density  scale 


- --.waiij  . 

The  ANGSUN  sub^outine^omp0^81111^6  fr°m  8Un* 
of  the  sun,  or  and T  P  63  the  P°8ition 
o  o* 


0  = 

-1 

sin 

U 

z 

9  = 

tan 

U 

x/ 

(2)  Compute  log  p  (h^  T)  for  four  altitudes> 

(3)  or")!13'6  f°r  108  '  (h’  T)  from  «»ults 


(4) 

(5) 


Compute  P  from  log  p  (h,T) 

Compute  the  velocity  relative  to  the 

vx  =  *  +  ey 

Vy  -  y  -  9x 


atmosphere: 


V 


+  V  2+  V  2 

y  z 


146 


JNDRAG 

spirdec 

3  of  3 


DESCRIPTION: 

(continued) 


(6)  Compute  the  perturbative  accelerations 
due  to  drag. 


XDTGR  =  = 

V 

X 

(B  P 

H  i 

6376.165) 

YDTGR  =  = 

V 

y 

(B  p 

h  1 

6378.165) 

ZDTGR  =  = 

V 

z 

(B  P 

*.  1 

6378.165) 
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LPAFT 

SPJRDEC 


PURPOSE:  Given  a  quantity,  X  in  radians,  this  subroutine  will  compute 

Y  and  Z  such  that: 

1.  Y  and  Z  have  the  same  sign  as  X 

2.  Y  is  exactly  divisible  by  2  fr 

3.  -2~  <  Z  f  2^ 

4.  Y  +  Z  =  X 


CALL  SEQ:  ENCE: 

TMA  X 

JMP  LPART 

INFl'T: 

The  input  consists  only 
in  floating  point. 

of  the  quantity  X  in  the  A 

register 

OUTPUT: 

Upon  return  from  LPART, 
Q  register  contains  Z. 

the  A  register  contains  Y,  and  the 
Both  Y  and  Z  are  in  floating  point. 

SUBROUTINES: 

None 

STORAGE 

REQUIREMENTS: 

14  Words 
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LSQ 

SPIRIEC 


PURPOSE : 

To  zero  matrix  buffer  All  to  All  +  48 
and  Bll  to  Bll  +  6 

CALL  SEQUENCE: 

JMP  LSQ 

INPUT: 

All  and  Bll  buffers  must  be  located  sequentially 
in  c're. 

OUTPUT: 

F/O  in  buffer  All  to  All  +  48  and  Bll  to  Bll  +  6 

SUBROUTINES : 

'None) 

STORAGE 

REQUIREMENTS : 

4  Cells 

DESCRIPTION: 

See  purpose. 
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LSQR 

SPIRDEC 


PURPOSE : 

CALL  SEQUENCE: 
INPUT: 


OUTPUT. 

SUBROUTINES: 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


A= 


To  form  a  least  squares  matrix. 

JMP  LSQf 

C,  thru  C  ,  in 
1  n+1 

cells  TERMS  thru  TERMS  +  n 

Modified  matrix  in  All  buffer  and  Bll  buffer. 
None 

15  cells 


Add  values  in  TERMS  buffer  to  matrices  A  (All  buffer) 
and  B  (Bll  buffer) . 


-C1C2 

■  ■  •  :Cj 

C  ~ 
n 

~:c'  Vi" 

:c2  ct 

| 

'C2  C2 

'  C2 

C 

n 

:  C2  c„+i 

1 

I 

1 

1 

B= 

1 

1 

1 

I 

•.  I 

f 

1 

1 

I 

! 

'c  c. 

n  1 

Cn  C2  ' 

7  P 

•  t  •  ^  V 

n 

C 

n 

r.c  c  , . 

n  n+l 

— 

_ , 

151 


LSQS 

SPIRDEC 


PURPOSE : 

To  solve  the  matrix  equation  AX  =  B  for  X. 

CALL  SEQUENCE: 

JMP  LSQS 

INPUT: 

A  matrix  in  All  buffer 

B  matrix  in  Bll  buffer 

OUTPUT: 

Solution  X  in  LSQX  buffer 

SUBROUTINES: 

Philco  -  FMAIN,  FMAMU 

Program  -  MATRIX 

STORAGE 

REQUIREMENTS: 

15  cells 

DESCRIPTION: 

The  top  half  of  matrix  A  is  accumulated  at  each  entry 
to  LSQR.  When  entry  is  made  to  LSQS,  the  terms  of  A 
are  moved  to  the  bottom  half.  Then  the  Philco  matrix 
inversion  subroutine  (FMAIN)  and  matrix  multiplication 
subroutine  (FMAMU)  are  used  to  solve  the  equation  for 

X.  The  solution  is  left  in  buffer  LSQX  through  LSQX+(n-l) . 
After  the  A  matrix  is  inverted,  subroutine  MATRIX  is 
called  to  solve  for  the  standard  deviations  and  correlation 
matrix  of  the  delta  elements. 

(See  LSQR  for  matrix  definition.) 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


MARTINI 
SPIRDEC 
I  of  5 


To  compute  the  perturbative  effects  of  the 
earth's  bulge. 

JMP  MARTINI 


BFLA.G  j  = 

0  do  not  compute  perturbations 

1  = 

1  use  perturbations 

W  +  1 

t(min) 

AE  =  a 

e 

X  =  x 

UZ  =  U 

z 

Y  =  y 

R  =  r 

THGR  =  0 

r 

gr 

TESRAL  1  = 

=  0,  use  only  zonals 

1  = 

:  1,  use  zonals  &  tesserals 

Buffers: 

JBUF,  CSNM 

XBDGR  = 

x>. 
x  B 

YBDGR  = 

»'b 

ZBDGR  = 

*'b 

and  values 

listed  under  description 

System  -  ARCTAN 

Philco  -  FSIN,  FCOS,  FSQRT 


169  Cells 


RT  1MUZ2  =  Vl-U  2 

_  z 


if 

z 

___  Ci 

*  0 

If 

Vl-U  2  =  0 

z 

SSUBX 

sx  -  -x 

s 

X 

=  F/0 

SSUBY 

s  =  u  r  /Vi-u  2 

S 

=  -U 

y  y  z  z 

y 

z 

SSUBZ 

sz  =  -  Vl-U  2 

s 

z 

=  F/0 
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DESCRIPTION: 

(continued) 


MARTINI 
SPIRDEC 
2  of  * 


If  Vl-U  2  ft  0 

Z 

If 

Vi-u  2  =  0 

z 

ESUBX  =  E 

X 

-u  /  Vi-u  2 
y  z 

E 

X 

=  F/l 

ESUBY  =  E 

y 

u  /  Vi-u  2 

X  z 

E 

y 

=  F/0 

ESUBZ  =  E 

F/O 

E 

=  F/0 

ZSUBX 

— 

Z 

X 

=  U 

X 

Z  =  U 

X  X 

ZSUBY 

= 

Z 

y 

=  U 
y 

Z  =  U 

y  y 

ZSUBZ 

= 

z 

z 

=  U 

z 

z  =  u 

z  z 

Compute  P^  and  P 

1  for  zonals  for  n 
n 

=  2,  3,  4,  5 

1)  P 

n 

= 

1/ri 

(2n-l) 

U  P  ,  - 

z  n-1 

<B’W  Pn-2} 

where 

P 

0 

F/l,  and 

P1  -  Uz 

p' 

n 

= 

U  P 
z 

.  +  n 
n-1 

P„-l 

where 

P/1  = 

F/l 

QQ 

= 

ae/r 

QQ 

+ 

1  = 

.  2 
li/r 

If 

TESRAL 

=  0,  only  compute 

zonal  effects: 

2)  GE 

= 

8et 

=  F/0 

GS 

- 

gst 

=  F/0 

l  tesseral  effects  are 

GU 

= 

8ut 

=  F/0 
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MARTINI 
SPIRDEC 
3  cf  5 


DESCRIPTION: 

(continued) 


g;z 

=  guz 

U  5 

=  ~r  e 

r 

(n+1)  J 

("a 

n 

P 

* 

n=2 

n 

F 

n 

GEZ 

gez 

=  F/0 

n 

GSZ 

g3Z 

=  u  Vl-U  2 

*  5 

i  j 

'a  ' 

e 

2 

n=2  n 

r 

n 

r 


"here  ^  .  J $  are  stored  in  JBUF  +  0-*  JBUF  +  3  consecutively 


4>  XBBGR  .  ^  ^  +  gut)  2x  +  ^  +  get)  Ex  +  ( ^ 


YBDGR  Yb  <8uz  +  Zy  +  <get  +  get)  Ey  +  ^st+Sst5  Sy 

ZBDGR  -  zB'  =  (g.j2  +  gut)  +  (gez  +  g£t)  Ez  +  (gsz+gst)  ^ 

(and  exit) 


,f  TESERAL^  0  use  the  tesseral  harmonics.  Substitute 
fctw  foilowi-nS  for  step  (2)  and  go  to  step  (3). 


Compute  Rpm  and  R^  for  all  2  <  n  <  4 


<  n  <  4  and  m  <  n,  where 


nm 


P 

r  nm 


and  R/  -Vf-'  ^  p' 
nm  z  rnm 


and  are  solved  as  follows: 
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MARTINI 
SPIRDEC 
4  of  5 


DESCRIPTION: 

(continued) 


R22’  R33’  R44  :  R  = 
-  mm 


(2m  -1) 


m-1 


where  R^  =  F/l 

R2i»  r31»  r41  .  R  =  Py 
-  n,l  n 


Rn,  n-1  =  (2n-l)  UR.  . 
*  z  n-1 ,  n-1 


R 

nm 


1 

n-m 


-(n-fm-1)  R  _  +  (2n-l)  U  R  .. 

n-2  ,m  z  n-l,m 


This  formula  is  the  general  expression  for  R  but  is 

nm 

only  used  for  R^  to  save  program  space  and  time. 


R'  =Vl-U  2  R 


nm 


n,nrt-l 


mU  R 
z  nm 


where  R  =  F/0 

n,n+l 


The  values  of  R  are  stored  in  buffer  RBUF+0  -  RBUF+8, 
nm 

and  for  R^  in  buffer  RPBUF+0  -  RPBUF+7 . 
nm 

TEMP2  =  XE  =  tan"1  (y/x)  -  t  (RPTIM)  -  0 

gr 


Then  build  buffer  SINCOS  containing  cos  X,  sin  X  .... 
cos  4X ,  sin  4X  consecutively  in  SINCOS+O  -  SINCOS+7  , 
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MARTINI 
SPIRDEC 
5  of  5 


DESCRIPTION: 

(continued) 


Solve  the  following  for  the  tessera!  harmonics: 
where  2  <  n  <  4  and  men 


ut 


=  -  JtVTi 


4  n 


.n 


GU  =  g.  =  -  ~n  * 1-U  Z  I  (n+1) 


n=2  m*l 


R  (C  cos  m>  +S  sin  m} ) 
cm  nm  nm 


GE  =  g 


et 


*  z 


n 

Z 


r  n=2  m=l 


in 


R  (C  sin  mX 
nm  nm 


-S  cos  m>.)  (m) 
nm 


GS  =  g 


st 


"  2  L 


r  n=2  m; 


S  (if 

fI  \  r  I 


Rx  (C  cos  m>.  +  S  sin  mX) 
nm  nm  nm 


Where  C22  S£2  . C ^  are  stored  consecutively  in 

CSNW-C-CSNMf 15  .  These  values  must  be  input  on  P  cards 
7  and  8  and  will  be  a  number  or  F/0. 


Having  solved  these  equations,  go  to  step  (3). 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


MATRIX 
SPIRDEC 
1  of  2 


To  compute  a  correlation  matrix  and  the 
standard  deviations  of  the  delta  elements. 

JMP  MATRIX 

A  •*  matrix  from  LSQR  subroutine  in  buffer  All 

COUNTR  =  C/HLT,  0;  C/HLT,  N  N  =  Matrix  Size  1  £  N  <  7 

Standard  deviations  in  buffer  SIGN 

Correlation  matrix  in  buffer  MATRIXB 

(See  description) 

Philco  -  FSQRT 


53  Cells 

To  avoid  destroying  the  A  ^  matrix  in  All  buffer, 
the  half  matrix  is  moved  to  buffer  MATRIXB. 

The  variance-covariance  matrix  is  defined  by 

c2i. 

^  —  where  C.  .  are  elements  of  A 

C..C..  1J 

n  JJ 

To  compute  this,  each  column  is  divided  by  the  root 
of  the  diagonal  term,  and  then  each  row  is  divided 
by  the  root  of  the  diagonal.  The  roots  of  the 
diagonal  terms  (which  are  the  standard  deviations) 
are  stored  in  buffer  BUF  and  are  then  moved  to  the 
appropriate  locations  in  buffer  SIGN  for  output 
purposes.  The  position  in  the  buffer  is  dependent 
on  the  elements  being  corrected.  (See  PRNTMAT 
subroutine  for  buffer  positions).  MATRIX  is  called 
by  subroutine  LSQS . 
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MATRIX 
SPIRDEC 
2  of  2 


DESCRIPTION: 

(continued) 


The  correlation  matrix  is  defined: 


11 


C11  C11 


12 


C11  C22 


ii 


cn  cjj 


c 


C22  °22 


il 


V 


C22  Cjj 


11 


C.  .  C. . 

jj  n 


1  £  j  <  7 


The  values  are  stored  by  column  in  MATRIXB, 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 

SUBROUTINES: 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


MOD 2 PI 
SPIRDEC 


To  Modulate  a  number  between  0  and  2  tt  radians. 

TMA  (number) 

JMP  MOD 2 PI 
TAM  (number) 

(A)  reg  =  number  to  be  modulated  2  tt 
TWOPI  =  2  tt  radians 

(A)  =  number  (mod  2  tt) 

(None) 


3  Cells 

Adds  or  subtracts  2  tt  radians  from  the  number  in 
the  (A)  reg  until  it  is  between  0  and  2  tt  radians. 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 


OUTPUT: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


MOVBUF 
SPIRDEC 
1  of  2 


To  move  the  observation  buffer, 

TMA  OCOUNT 
JMP  MOVBUF 

Observations  in  10  word  format  and  in  6  word 
format  starting  in  location  EBLOC. 

Left  address  of  OBSREJ  =  number  of  observations  rejected 

Left  address  of  OBSLEFT  =  number  of  observations  in  10  word 

fomat 

Left  address  of  0BSPR0C  =  total  number  of  observations 

processed 

Left  address  of  OCOUNT  =  number  of  observations  processed 

to  6  word  format  in  one  group 

0A1DE  =  C/HLT,  OBLOC;  C/HLT,  EBLOC 
See  description 


13  Cells 

Because  of  core  limitations, the  observations  in  the 
10  word  format  of  OBLOC  are  converted  to  a  6  word 
format  (see  below).  Initially  up  to  984  observations 
(10  words/observation)  are  in  EBLOC. 

Subroutine  PROOBS  controls  the  formatting  of  the 
observations  into  6  word  entries.  Subroutine  MOVBUF 
is  called  by  PROOBS  and  OBVEC  (which  is  called  by  PROOBS). 

Procedure: 

(1)  Add  the  number  of  observations  that  have  been  converted 
to  a  6  word  format  (OCOUNT)  to  the  number  of  observations 
rejected  (OBSREJ).  Multiply  this  number  by  10  and  set  up 
the  move  instructions  at  location  MOVSW,  if  the  result 
is  less  than  4096;  otherwise  print"more  than  77  observa¬ 
tions  rejected"  and  exit. 

(2)  Multiply  the  number  of  observations  still  in  the  10  word 
format  (OBSLEFT)  by  10  and  add  it  to  the  total  number  of 
observations  in  a  6  word  format  (OBSPROC)  which  is 
multiplied  by  6.  This  number  is  stored  in  OCOUNT  as  the 
total  number  of  cells  to  be  moved  up. 
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DESCRIPTION: 

(continued) 


MOVBLT 
SPIRDEC 
2  of  2 


(3)  After  moving  the  buffer  up, 

set  OCOUNT  =  0 

OBSREJ  =  0 

WOBMARK  =  C/HLT,  OBLOC;  C/HLT,  EBLOC 

(4)  Store  Z's  in  the  next  available  location  at  the 
end  of  the  buffer  and  move  this  location  to  Index 
Register  0. 

_ _  EBLOC 

Up  to  984 
observations 
in  OBLOC 

format  (10  words/obs.) 


Up  to  332 
observations 
in  6  word 
format 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 


SUBROUTINES: 

STORAGE 

REQUIREMENTS: 


MOVDAI 
SFIRDEC 
1  of  2 


To  retrieve  next  processed  observation. 


JMP  M.VDAT 

JMP  (End  of  observations) 
(Return) 


Processed  observations  in  EBLOC . 

Current  address  in  EBLX  from  location  SAVOBS. 
Sensor  data  in  SBUF. 

Weights  in  WBIF. 


STAID 

OBFLG 

T 

CAPX 

CAPY 

CAPZ 

CXDOT 

CYDOT 

RANGE  A 

SIGMA 1  J 

ASUBX 

ASUBY 

ASUBZ 

DSUBX 

DSUBY 

DSUBZ  > 

XLSUBX 

XLSUBY 

XLSLEZ 

SIGMA3 

SIGMA4  / 

SIGMA2  \ 

RODOT  / 


if  range  observed 


if  angles  observed 


if  range-rate  observed 


Program  »  GETSEN,  SENL-OC ,  AZREC ,  ALREC ,  GETWC-T 


27  Cells 
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MOVDAT 
SPIRDEC 
2  of  2 


DESCRIPTION; 


1)  Unpacts  the  following  from  processed  observation 

STAID 

OBFLG 

T 

RANGE 

RODOT 

ALPHA 

DELTA 

2)  Use  subroutine  GETSEN  to  retrieve  sensor  data. 

3)  Use  subroutine  SENLOC  to  compute: 

THTA  =  0 

SINTH  =  sin  0 
COSTH  =  cos  8 

4)  If  azimuth  and  elevation  were  observed,  call 
subroutine  AZREC  to  compute  A,  D,  L. 

5)  If  right  ascension  and  declination  were  observed, 
use  subroutine  ALREC  to  compute  A,  D,  L. 

6)  Compute  the  following: 

CAPX  =  X  =  (cos  8)  (X/cos  0) 

CAPY  =  Y  =  (sin  9)  (X/cos  8) 

CXDOT  =  X  =  -Y.0 

CYDOT  =  Y  =  X  9 

7)  Retrieve  the  weights  for  this  observation  using 
subroutine  GETWGT. 

8)  Update  SAVOBS  to  next  observation  and  exit. 
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SrJklEC 
1  cr  2 


PURPOSE : 

To  retrieve  observations  and  corresponding 
weights . 

CALL  SEQUENCE: 

jMP 

NXTOB 

JMP 

(End  of  Obs) 

JMP 

(No  Sensor  Data) 

JM? 

(No  Weights; 

JM? 

(Normal  Recurm 

INPUT: 

OSTRL 

B  -  right  aaaress  =  nomter  cf  ots  +  i 

Observations  in  EBi.il 

Weigh 

ts  ir.  WE1F 

OUTPUT: 

STAID 

T 

DELTA 

ALPHA 

RANGE 

RODGT 

and  output  of  subroutine  GETSEN 
SUBROUTINES:  Program  -  GETSEN 

STORAGE 

REQUIREMENTS:  39  Cells 

DESCRIPTION:  (1)  Inpact  the  fo.j-owing  from  the  torrent  observation 

in  EBLOC: 

cta::d  -  oooocsss 

E  ?•_>'.  'A 

AL’-'Ha  >  floating  y  tirt 

RANGE 

RODGT  . 

hi.  PE  -  CCQ00C0 

(2)  Search  WBUF  to  find  weights  for  the  sensor  number 
in  STAID.  If  no  match,  exit  +3H  (No  weights  re to 
Otherwise  gc  rc  step  (Vi  , 


1 65 


NXTOB 
SPIRDEC 
2  of  2 


DESCRIPTION: 

(continued) 


(3)  Check  to  see  that  a  weight  is  in  WBUF  for 
each  observed  quantity  of  the  observation. 
If  not,  exit  +3H. 


(4)  If  all  weights  are  entered,  call  subroutine 
■GETSEN  to  retrieve  the  sensor  data.  If  no 
sensor  data,  exit  +2H.  If  sensor  data,  exit 
+4H  (normal  return) . 


Each  call  to  NXTOB  will  retrieve  the  next  observation, 
until  an  end  of  observation  return  is  taken  (+1H) . 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


SUBROUTINES: 

STORAGE 
REQUIREMENTS . 


OBVEC 
SPIRDEC 
L  cf  2 


To  compute  and/or  store  values  in  FOBLOC  for 
one  observation  at  a  time. 


JMP  OBVEC 


OBFLC-  = 
STAID  = 
T 

RANGE “ 
RODOT 

> 

ALPHA 
DELTA . 


0 - 0  XXXX  (S  cc  PROOBS  for  description j 

00000  SSS  Station  number, 
time  (minutes  since  epoch) 

floating  point 


See  POBLOC  format  in  PROOBS 
T 

THTA  =  90 

SINTH  =  sin  9 
COSTH  =  cos  9 
CAPX  =  X 
CAF£  =  Y 

and  output  cf  OBVEC Ps  OBVECQ.  AZREC  ALREC ,  RRATE 
Pr  gram  -  MOVE Vi 


L5  CeLis 
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03VEC 

spirde: 

2  of  2 


DESCRIPTION:  1)  An  entry  is  stored  in  the  new  observation  format 

in  the  address  specified  by  index  register  0: 

Word  0  TOOOOSSS 

1  time 

2  range 

3  range  rate 

4  alpha 

5  delta 

2)  Add  1  to  0BSPR0C 

Subtract  1  from  OBSLEFT 
Add  1  to  MCOUNT 

a)  if  MCOUNT  <  332,  then  exit 

b)  if  MCOUNT  =  332,  call  subroutine 
MOVBUF  to  move  the  buffer  up. 

(See  MOVBUF  for  description.) 

Then  exit. 
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PAGE CON 
CT’IR  DEC 


PURPOSE : 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


To  keep  a  count  (TL5)  of  the  lines  output/page  and  force 
a  new  page  with  headings  when  count  >  55. 

IMA  C/HLT,N 
JMP  PAGECON 
JAZ  (  ) 

N  =  No.  of  lines  of  output,  set  N  *  55  to  force  page 
(  .*  )  =  0  means  new  page  was  output 

See  Call  Sequence 

(1)  Updated  line  count 
and/or 

(2)  New  page  with  headings  when  lin^  count  >  55. 
Program-HEAD 


7  Cells 

Updates  line  count  (LINECNT)  and  jumps  to  subroutine  HEAD 
(to  output  headings)  when  LINECNT  >55.  If  the  previous 
JMP  to  PAGECON  put  out  a  new  page  a  second  one  cannot  be 
forced . 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


PCONTRL 
SPIRDEC 
1  of  8 


To  control  the  prediction  option. 

JMP  PCONTRL 
JMP  (Error) 

(Normal) 

P  Card  4 

An  element  set  from  P  Cards  1  and  2  or  6  Card 
element  set  or  from  the  differential  correction. 

Prediction  ephemeris  as  requested  -  printed 
and/or  written  on  a  binary  tape. 

(See  description) 

System  -  FYKLOK 

Program  -  ITAPEW,  SUBOUT1,  SEPSUB,  INITIAL,  SAV5PTS , 

TEMP,  SAVW,  SETW,  SUBXYZ,  RESW,  TAPEW,  FTAPEW , 
SUBOUT,  AD  BASH,  REVSOB,  DIVDIF,  CALU,  GIPAR 

298  Cells 

This  subroutine  contains  the  logic  necessary  to  control  a 
prediction.  There  are  several  options;  some  can  be  combined 
but  not  all: 

(1)  Predict  by  revolution  number  or  time. 

(2)  Output  hard  copy  and/or  binary  tape. 

(3)  Prediction  within  one  or  more  time  intervals  at 
specified  time  increments. 

(4)  Predict  backward  or  forward  in  time  from  epoch. 

(5)  Request  points  (maximum  of  5)  to  be  left  in  core 
for  a  GIPAR  run. 
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DESCRIPTION: 

(continued) 


FCONTRL 
SPIRDEC 
2  cl  8 


Because  several  of  these  options  can  be  combined,  the 
logic  is  rather  complicated. 

A.  The  subroutine  is  initialized  as  follows: 

(1)  Set  several  switches. 

(2)  Initialize  the  buffer  to  be  used  by  GIPAR 
(starts  in  EBLOC  +  20) . 

(3)  Set  switches  controlling  a  prediction  by 
time  or  revolution  number. 

(4)  Initialize  the  binary  tape  if  requested. 

(5)  Initialize  the  printed  output  if  requested. 

(6)  Move  some  values  from  P  Card  4  to  locations 
used.  If  the  prediction  follows  a  D  .C . , 
same  results,  such  as  B,  override  the  input 
on  P  Card  4. 

(7)  Test  mode: 

(a)  If  prediction  only: 

(1)  Move  more  values  from  P  Card  4 

(2)  Compute  yB/2.2 

(3)  Set  REV  =  epoch  revolution  number 

(b)  If  D ,C .  and  prediction: 

(1)  Move  new  epoch  elements  from  buffer 
to  locations  used . 

(2)  If  time  prediction,  change  final  time 
from  minutes  since  epoch  to  minutes 
since  new  epoch. 

(3)  Convert  new  epoch  time  from  minutes 
since  epoch  to  a  base  epoch. 
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PGONTRL 
SPIRDEC 
3  of  8 


DESCRIPTION: 

(4) 

Initialize  page  heading  routine 

(continued) 

(IHEAD1)  to  print  new  epoch  in 
page  headings. 

(5) 

Call  subroutines  THGRC ,  BEGIN. 

(8) 

Set  REV  *  new  epoch  revolution 
number . 

(8)  Initialize  interpolation  buffer: 

(a)  Call  subroutine  INITIAL  to  initialize 
ADBASH . 

(b)  Save  initial  elements  as  first  set  in 
the  interpolation  buffer. 

(c)  Set  T  =  0  (printed  output  time) 

BT  =  0  (binary  tape  output  time) 

(d)  Call  subroutine  TEMP  to  compute 
temperature  at  epoch. 

(9)  Test  if  GIPAR  points  requested: 

(a)  If  requested,  prediction  interval  points 
will  be  ignored  -  if  not,  go  to  step  A(10) . 

(b)  Convert  GIPAR  time  from  Cays  and  fractions 
to  minutes  since  epoch. 

(c)  Delete  GIPAR  points  not  between  epoch 
(or  new  epoch)  and  the  final  time. 

(10)  Test  if  predicting  forward  or  backward  in  time 

and  set  switches  and  locations  accordingly. 

(a)  If  predicting  backward  and  if  requesting 
prediction  intervals,  changes  signs  of 
interval  buffer  to  negative. 
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PC  ONTRL 
SPIRDEC 
4  of  8 


DESCRIPTION 

(continued) 


(11)  Test  if  printed  output  requested: 

(a)  If  not,  go  to  step  B. 

(b)  If  yes,  test  if  GIPAR  points  requested 

(1)  If  not,  go  to  A(ll) (c) . 

(2)  If  requested,  set  a  switch  for  GIPAR 

logic  and  move  the  first  time  requested 
to  T. 

(c)  Test  if  prediction  intervals  requested 

(1)  If  not,  go  to  step  B. 

(2)  Otherwise,  move  start  time,  end  time, 
and  time  increment  to  T,  LTBUF+1  and* 
ENDT  respectively;  also  move  start 
time  to  BT. 


B.  Basic  integration  loop: 

(1)  Call  AD BASH  to  integrate  to  next  point. 

(2)  Test  if  drag  perturbations  are  greater  than  an 
epsilon;  if  so,  switch  to  Cowell  option  for  decay. 

(3)  Update,  revolution  number  if  necessary. 

(4)  Save  the  new  point  in  the  interpolation  buffer. 

(5)  If  5  points  are  in  the  buffer  for  interpolation, 
go  to  step  B(6);  if  not,  go  to  step  B(l). 

(6)  Test  if  current  revolution  number  equals  final 
revolution  number  if  a  revolution  prediction  was 
requested: 

(a)  if  not  equal,  go  to  step  C 

(b)  if  equal,  go  to  step  D(5). 
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PCONTRL 
SPJRDEC 
5  of  8 


DESCRIPTION 

(continued) 


C  Binary  tape  loop: 

(1)  If  no  tape  output  was  requested  go  to  step  D. 

(2)  Test  if  binary  tape  time  is  within  the  time  span 
of  the  interpolation  buffer. 

(a)  If  not,  test  if  printed  output  was  requested: 

( 1)  Go  to  step  B ( 1) ,  if  not  requested. 

(2)  Go  to  step  D,  if  requested 

(b)  If  in  the  range  of  the  interpolation  buffer 


(1)  Interpolate  for  elements  at  the  time 

(2)  Use  SUBXYZ  subroutine  to  convert  to 
r,  r 

(3)  Call  subroutine  TAPEW  to  output  the 
point 

(4)  Update  the  time 

(5)  Test  if  time  or  revolution  prediction 

(a)  If  time,  test  if  1  st  time 
>  final  time;  if  not  >  go  to 
step  C  (2) (b)  (5) (b) ;  if  equal  , 
wrap  up  the  binary  tape  and 
turn  off  tape  option;  if  >, 
set  t  equal  final  time,  and 
loop  to  C (2 )  once  more 


(b) 


For  revolution  or  time 
prediction  intervals; 
to  step  C (2) ;  if  yes , 
end  time;  if  <,  go  to 
if  >,  ■  t  time  =  end 

to  step  C (2) , 


,  test  for 
if  not ,  go 
test  if  time 
step  C  (2)  ; 
time  and  go 
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DESCRIPTION: 

(continued) 


PC  OMR  L 
SPIRDEC 
6  of  8 


D.  Printed  output  loop: 

(1)  If  time  prediction,  test  if  time  >  final  time: 

(a)  If  >,  go  to  D (5) . 

(b)  If  <,  go  to  step  D '2) . 

(2)  Test  if  time  is  in  span  of  interpolation  buffer 

(a)  If  not,  go  to  step  B. 

(b)  If  yes,  interpolate  for  elements  at  the 
time;  convert  to  r,  r,  and  output  the  point. 

(3)  If  GIPAR  points  were  requested,  call  subroutine 
GIPAR  to  convert  and  store  the  values 

(a)  If  revolution  prediction,  go  to  D (3)  (b) (2)  (a) . 

(b)  If  time  prediction  -  test  if  more  times 

(1)  If  not,  turn  off  prediction  interval 
option  and  go  to  D(5)(b)(for  revolution 
prediction)  or  to  D(5)(a)  (for  time 
prediction) . 

(2)  If  yes,  compare  next  time  to  final  time: 

(a)  if  £  put  next  time  in  T  and  go 
to  D(l) 

(b)  if  >  go  to  D  (3)  (b)(1). 

(4)  If  GIPAR  points  not  requested,  update  time,  and  test 
for  prediction  intervals: 

(a)  If  no  intervals,  go  to  D  .2) . 

(b)  If  intervals: 

(1)  If  time  £  end  time,  go  to  D(2) 
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DESCRIPTION: 

(continued) 


(2)  If  time  >  end  time,  set  time  =  end 
time  and  go  to  D(2)  the  first  time; 
second  time,  check  for  end  of  intervals, 
if  end,  go  to  D(5)(b)  (for  revolution 
prediction),  to  D(5)(a)  (for  time 
prediction),  if  not  end,  move  next 
interval  times  from  buffer,  reset  switch 
and  go  to  D(2) . 

(5)  Switch  at  D(l)(a) 

(a)  First  time,  set  time  to  final  time,  set 
switch  to  second  time,  go  to  step  C, 

(b)  Second  time,  test  if  tape  output: 

wrap  up  tape,  if  yes 
then  exit  +2H  (normal  exit) 


E,  Cowell  option. 

When  a  satellite  decays  to  about  90  km.  altitude,  the 

drag  coefficient  is  such  that  integration  should 

continue  in  the  Cowell  mode. 

(1)  A  test  is  made  for  backward  integration,  which 
is  not  allowed  for  decay. 

(2)  Since  the  interpolation  buffer  now  contains  N 
elements,  it  is  necessary  to  integrate  backw.nds 
for  sufficients  points  to  replace  the  N  ,M  elements 
with  r,  r,  for  the  same  time  range, 

(3)  Set  the  integration  to  go  forward  integrating  tor 

r  r  elements  until  the  interpolation  buffet  cunt ..it  is 
5  points . 

(4)  The  first  time  the  buf  fer  is  full,  save  the  inlet - 
potation  buffer,  time  for  restart  in  the  decay- 
corridor  , 

(5)  Test  for  tape  output,  If  requested,  test  if  the 
output,  time  is  in  the  range  of  the  interpolation 
buffer.  If  so,  interpolate  for  the  elements, 
output  the  point,  update  the  time.  Continue 
looping  until  the  time  range  is  excluded,  thin 
go  to  E (6)  , 
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DESCRIPTION: 

(continued) 


PC  ON7R1 
SIIRLEC 
8  of  8 


(6)  Test  for  printed  output.  Follow  the  same 
procedure  as  step  E(5).  Go  to  step  E(3) 
when  the  time  range  is  excluded . 

The  loop  of  steps (3) through (6)will  continue 
until  the  CDERIV  subroutine  (called  by  ADBASH) 
exits  to  step  E(~).  This  is  done  when  the 
vehicle  drops  below  10  km.  At  this  point,  the 
decay  corridor  is  produced. 

(7)  Print  comment  that  -vehicle  decayed. 

(a)  If  first  entrance,  double  the  value  of  B  and 
print  a  comment  to  this  effect.  If  tape  output 
was  requested,  wrap  up  the  epheir.eris  tape  and 
turn  off  the  tape  output  option.  Restore  the 
values  saved  at  step  E(4;  and  go  to  step  E(5) 
to  produce  an  empheris  for  a  drag  coefficient 
of  twice  the  original  value. 

(b)  If  second  entrance,  multiply  the  original  B 
by  1/2,  print  a  comment  and  follow  the  same 
procedure  as  E(7)(a)  only  with  a  different 
value  of  B. 

(c)  If  third  entrance,  exit  (+2H)  the  decay 
corridor  and  prediction  are  completed. 


177 


\START  / 
PCONTRL 


A 


FIGURE  8.  PCONTROL  FLOW  DIAGRAM  (1  of  t 

178 


G*n  Adm 
Jon  60 


ARO‘7101 


TART  J 
PCONTRL 


SET  ALL 
SWITCHES 
TO  1 


INITIALIZE 
BUFFERS 
FOR  GIPAR 
POINTS 


SET  A  SWITCH 
DEPENDENT  ON 
TIME  C»  REV. 
PREDICTION 


s binary 

OUTPUT  TAPE 
REQUESTED 


r  HAS  X 
A  D.C . 
BEEN 

PERFORMER 


UNPACK 

MORE 

PARAMETERS 
FROM  P  CARD  4 


COMPUTE  a 
INITIALIZE 
REV.  NO. 


MOVE  NEW 
EPOCH 
ELEMENTS 


PREDICTION 
BY  REV. 
VOR  TIME  > 


SET  FIN# 
TIME  RE  LAI 
TO  STAR! 
TIME 


COMPUTE 
CLASSICAL 
ELEMENTS  AT  / 
NEW  EPOCH  / 


SET  REV.  NO. 
TO  REV.  NO. 
OF  NEW 
EPOCH _ 


INITIALIZE 

INTEGRATION 

BUFFER 


IN 

INTERPOLATION 
V  BUFFER _ / 


GIPAR  \ 
POINTS 
REQUESTED 


CONVERT  GIPAR 
TIMES  TO 

DELETE 

GIPAR  POINTS 

MINUTES  SINCE 
EPOCH 

NUI  WllrllN 
PREDICTION 

TTME  SPAN 

A/  X 

PREDH 
^  BY  TIM 


REVOLU 


BACKWARD 


''PREDI 
FORWA 
OR  BAC 


SET  AD BASH 
TO 

INTEGRATE 

BACKWARD 


'PRINTED  \ 
OUTPUT 
REQUESTED 


Gen  Adm 

Jon  6C 


AR0-7101 


FIGURE  8, 


POONTROL  FLOW  DIAGRAM  (I  of 

178 


REDICT 

ORWARD 

BACKW 


9 


^ G I  PAR  X 
POINTS 
REQUESTED 


SET 

GIPAR 

SWITCHES 


AD BASH 
INTEGRATE 
FROM  T 
TO  T  +  ^T 


/  SWITCH N 
TO  COWELL 
EQUATIONS 


-  1  1  * 

SET  NO  GIPAR 

/^REDICTIOlX, 

v  YES 

MOVE  FIRST 
INTERVAL  SET, 

SWITCH 

— intervals 

UPDATE 

- 

X^EQUESTETy' 

ADDRESSES 

. 

/  5  \ 

yr  POINTS  X. 

[N  INTERPOLATION 

S  BUFFER 


TREDICT 
BY  REV. 
OR  TIME 


PCTAPE 

SET 

SWITCHES 
TO  1 


BINARY‘S 

TAPE 

REQUESTE- 


XIS  \ 

CURRENT  X 
REV.  EQUAL 
.FINAL  REV*'' 


SAVE 
CURRENT 
TIME _ 


DIVDIF 
INTERPOLATE 
FOR  OUTPUT 
POINT 


SAVW 

SAVE  THE 
INTEGRATION 
BUFFER 


SETW  \ 

/  1 

SET  THE  \ 

/  C‘ 

INTEGRATION  V-i 

n  x, 

BUFFER  TO  / 

\ 

THE  OUTPUT  / 

\  X. 

/PREDIC1 

BY 

TIME  OR 
X  REV .. 


/  IS\ 
CURRENTS. 
'  TIME 
^GREATER  THAN, 
— X  FINAL 
NO  X  TIMZ^ 


FTAPEW 

WRAP  UP 
THE  BINARY 
OUTPUT 


TURN  OFF 
THE  BINARY 
TAPE  OUTPUT 
OPTION 


SET  T  TO 
FINAL  TIME 


''PREDICTION  > 
INTERVALS 
.REQUESTED  ^ 


/CURRENT  X. 
TIME  PAST  X^ 
THE  END  OF  THI 
XlNTERVAL^ 


REV, 

sw  \ 

C  SWITCH 


=  1 

"time 


LAST  I 
TIM 


Gen  Adm 
Jan  60 


ARD-7101 


A 


FIGURE  8.  PCONTRL  FLOW  DIAGRAM  (2  of  4) 

179 


1 


DIAGRAM  (3  of  4) 
180 


PHLAH 

SPIRDEC 


PURPOSE: 

CALL  SEQUENCE: 


To  compute  CP,  X,  h 
JMP  PHLAH 


INPUT: 

UX 

PI 

=  U 

z 

=  it  radians 

XIMFSQ 

-  (1-f)2 

TWOPI 

=  2  tt  radians 

X 

=  X 

Y 

=  y 

RPTIM 

=  rad/solar  min. 

W  +  1 

t  (min.  since  epoch) 

THGR 

=  9 

gr 

OUTPUT: 

PHI 

-  ®  (deg) 

XLAMD 

-  \  (deg)  -180°  (W)  < 

H 

=  h  (km) 

LFLAG 

=  E  or  W  for  X  (BCD) 

SUBROUTINE: 

Program  - 

MOD  2PI ,  CALH 

System  - 

ARC TAN 

Philco  - 

FSQRT 

STORAGE 

REQUIREMENTS: 

18  Cells 

DESCRIPTION: 

PHI  =  cp 

=  tan"1  /Vl  -  Uz2 

THETA  = 

tan_1(y/x) 

(1  -  f )• 


XLAMB  =  -  t  (.0043752691)  -  THGR  +  THETA 


Uses  CALH  subroutine  to  compute  h. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT. 


OUTPUT: 


SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


PRERES 
SPIRDEC 
1  of  3 


To  retrieve  or  compute  tne  new  epoch  elements. 


JMP  PRERES 

SOLDUZ  = 
SREV 

PRDSAV  = 

EPREV 
SAVELEM  = 


last  U  from  DC 
z 

last  rev.  no.  from  DC 
0  means  new  epoch  elements  found  in  DC 
1  continue  integration  until  new  epoch  time 
epoch  rev.  no. 

0;  continue  DC  integration  to  compute  elements 
1;  initialize  integration  to  compute  elements 


and  input  required  for  subroutines. 


xn 


yn 


LPRINT  =  L  (deg) 

AX PRINT  =  a 
AYPRINT  = 

HX  PRINT  =  h 
HYPRINT  =  h 

3 

HZPRINT  =  h" 

2 

BPRINT  =  B 
HSUBQP  =  h 
PAPRINT  = 
and  time  in  BCD  format. 


(km) 


Pa  (min) 


Program  -  INITIAL,  RESKBF ,  RESICK,  ADBASH ,  REVSUB, 
SAV5PTS,  DIVDIF,  SAVICK,  SAVWBF ,  BCDIIM 


51  Cells 

There  are  2  main  paths  to  take  in  this  subroutine: 

(1)  to  retrieve  the  new  epoch  elements  which  were 
computed  in  the  DC,  i.e.,  the  new  epoch  elements 
were  in  the  span  of  the  observations ,  (2)  to  continue 
the  integration  started  by  the  DC  until  the  new  epoch 
time  is  reached . 
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PRERES 
SPIRDEC 
2  of  3 


DESCRIPTION: 

(continued) 


If  the  new  epoch  elements  were  found  in  the  D.C, 
(PRDSAV  =  0),  the  values  that  were  saved  at  the 
time  are  restored  and  the  rest  of  the  values 
computed: 


Restore:  BOO  +  0  =  a 


=  a  — 

xno 

- > 

AX PRINT 

—  a  _ 

AY PRINT 

yno 

- 

=  b  _ 

BPRINT 

=  L  (rad) — 

— > 

LPRINT 

_  , 

HX PRINT 

—  [  l  ~ 

X 

> 

_  h  _ 

hyprint 

y 

-  h 

hzprint 

Z 

_ 

Compute: 

HSUBQP  =  h^  (km)  =  |a  (1-e)  -  lj  6378.165 


km 


e  .r 


where  BOO  +  3  =  a 
BOO  +  4  =  e^ 


PAPRINT  =  P  (min)  =  2n 
a 

1  (l-.5e2) 

2 

(3  sin  i- 1)  P3JA02  +  1 

2 

2 

1 

P 

i 

n 

o 

where  BOO  +  5  =  sin  i 


BOO  +  6  —  p 
BOO  +  7  =  n 

o 

Use  subroutine  BCDTIM  to  get  time  in  BCD  format, 
where  BOO  +  8  =  t  (min  since  epoch)  . 
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PRERES 
SPIRDEC 
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DESCRIPTION: 

(continued) 


II.  If  the  new  epoch  elements  have  not  been  found, 
then  test  to  see  if  the  D.C.  integration  muse 
be  initialized  (SAVELEM  *  0).  If  so,  the  Adams- 
Bashforth  subroutine  must  be  initialized,  then 
integration  will  begin  at  epoch.  If  not,  the 
Adams -Bashf or th  buffer  will  be  restored  to  the 
time  of  the  last  observation,  then  integration 
will  continue  until  the  elements  are  found  for 
the  new  epoch  time . 

Having  found  the  new  epoch  elements,  the 
procedure  will  be  the  same  as  that  listed 
under  I. 

This  routine  is  called  only  by  the  ELMOUT 
subroutine  when  the  elements  are  to  be  printed. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


PR  INTO 
SPIRDEC 


To  print  the  weights  and  biases  being  used 
in  the  differential  correction. 

JMP  PRINTO 

Weights  in  WBUF 
Biases  in  BIBUF 

WANDBI  =  C/HLT,  WBUF;  C/HLT,  BIBUF 

Weights  and  biases  printed  on  hard  copy. 
(See  description) 

Program  -  PAGECON 
System  -  PANT,  GLOP 


92  Cells 

This  subroutine  is  called  by  option,  if  Column  3 
on  P  Card  3=1.  The  weights  and  biases  are 
printed  as  follows: 

(1)  Print  heading 

(2)  Retrieve  an  entry  from  WBUF  and  BIBUF. 

Convert  to  output  format. 

(3)  Output  the  Sensor  Number  and  the  weights 
and  biases. 

(4)  When  a  new  page  is  necessary,  go  to  step  (1) . 
Otherwise,  go  to  step  (2). 

The  subroutine  will  exit  when  a  word  (00000ZZZ)  is 
found  in  WBUF  . 


186 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


PROOBS 
SFTRIEC 
1  of  2 


To  set  up  the  processed  observation. 
JMP  PROOBS 


Observation  biases  in  BIBUF  buffer,  stored  by 
station  number. 

Observations  -  start  location  is  EBLOC*  terminated 
by  Z's. 

file  in  SBLOC . 

?!  C/,HLT»  N+1»  N  =  number  of  observations. 
C/HLT,  (Address),  Address  =  1st  location  to  be 
used  for  a  processed  observation 
PREDFLG  =  new  epoch  option 


Sensor 

OSTROB 

SAVOBS 


Processed  observation  buffer,  modified  sensor  buffer. 
(See  description) 


Program  -  BIAS,  SORTOB,  SETSBUF,  NXTOB,  OBVEC,  PAGECON 
MOVBUF  N , 

System  -  GLOP 


55  Cells 

(1)  Use  BIAS  subroutine  to  make  a  pass  through  the 
observations  and  apply  biases. 

(2)  Call  subroutine  SORTOB  to  sort  the  observations 
in  the  order  to  be  processed. 

(3)  Use  subroutine  SETSBUF  to  make  a  pass  through 
the  observations  and  move  sensor  information 
required  co  SBUF. 

(4)  Sec  0C07NJ,  MCO,NT,  OBSPROC,  OBSRF;  =  0/0 

Set  OBSLEFT  =  C/HLT,  N;  N  =  number  of  observations 

(5)  Retrieve  an  observation  using  NXTOB: 

(a)  if 'no  sigmas"return,  skip  observation, 
print  comment,  subtract  1  from  OBSLEFT 
and  add  1  to  OBSREJ. 

(b)  if  no  sensor  return,  follow  same  procedure 
as  (a) . 
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PROOBS 
SPIRDEC 
2  of  2 


DESCRIPTION:  (6)  Set  the  bits  in  OBFLG  corresponding  to  the 

(continued)  observed  quantities: 

1/1T47  -  P 
1/1T46  -  A  &  h 
1/1T45  -  a  &  6 
1/1T44  -  P 

(7)  Use  subroutine  OBVEC  to  compute  values  and 
store  them. 

Continue  steps  (5V0)  until  an  "end  of  observations 
return"  from  NXTOB;  then 

(8)  If  OCOUNT  t  0,  use  subroutine  MOVBUF  to  move 
remaining  observation  buffer  up. 

(9)  If  PREDFLG  =  2,  set  PRTIM  =  time  of  the  last 
observation. 

See  the  listed  subroutines  for  output  format  and 
block  formats. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


PRNTMAT 

SPIRDEC 


To  print  the  standard  deviations  and  the 
correlation  matrix  of  the  delta  elements. 

JMP  PRNTMAT 


SCOUNTR 

= 

C/HLT 

SIGN  +  0 

rs 

o  n 

1 

= 

a  a 

xn 

2 

o  a 

yn 

3 

= 

<7  U 

o 

4 

= 

o  n 

5 

= 

°  i 

6 

= 

B 

0,  C/HLT,  N 


Where  N  matrix  size 


MATRIXB  -  MATRIXB  +  27 
Contains  the  correlation  matrix. 

It  is  of  variable  size  depending  on  the 
elements  being  corrected. 


See  description. 


System  -  PANT,  GLOP 


74  Cells 

Prints  headings  and  values  for  the  standard 
deviations.  Then  print  the  headings  for  the 
matrix  and  N  lines  of  the  matrix  as  specified 
by  location  SCOUNTR. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


RDPRES 

SPIRDEC 


To  compute  the  perturbative  acceleration  due  to 
direct  solar  radiation  pressure,  r ^ . 

JMP  RDPRES 


RPFLAG 


no  perturbations 
compute  perturbations 


RPC0N3  =  -v-B/2.2 
XRDGR  =  xx 

r 

YRDGR  =  y' 

]  r 

ZRDGR  =  z' 

r 

Philco  -  FSIN,  FCOS 
Program  -  ANSSUN 


34  Cells 

If  RPFLAG  =  0,  the  subroutine  sets  XRDGR,  YRDGR,  and 
ZRDGR  =  F/0  and  exits.  Otherwise  it  computes  the  following: 

Call  subroutine  ANGSUN  to  compute  L^,  then: 


cosy  = 

i0  •  £ 

r 

If  cos 

♦  >  0, 

then  satellite 

Xs  = 

r 

(  v  • 

W.l)  Lx0 

n 

>-■ 

/ 

• 

(  y  • 

B/2-2)  Ly0 

z'  = 

r 

(  v  • 

E/2.2)  Lz0 

If  cos 

+  <  o, 

then  compute: 

sin  (  i|r  +  T)  )  =  sin  y  cos  71  +  sin  T|  cos  \|r 

=  (  cos^y  -1)  (l/r^-l)  +  cos  \|r / r 

If  sin  (  \li  +  T  )  >  0,  the  satellite  is  illuminated  so  compute 

r  .  Otherwise  rN  is  set  =  0. 
r  r 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


RDISB 

SPIRDEC 


To  compute  range  rate. 


JMP  RDTSB 


XDOT  =  x 
YDOT  =  y 
ZDOT  =  z 
CXDOT  =  X 
CYDOT  =  Y 
CZDOT  =  Z 
XLX  =  L 

x 

XLY  =  L 

y 

XLZ  =  L 

z 

RHODT  =  <5 

c 

RODTX  =  p 

x 

RODTY  =  p 

.y 

RODTZ  =  p 

z 

^None; 


8  Cells 

FODTX  = 

XDOT  +  CXDOT  = 

P  = 

X 

+  X 

X 

rctty  = 

i DOT  +  CYLOI  = 

•  TJ  . 

II 

y 

+  Y 

RODTZ  - 

ZTOT 

P  = 

z 

z 

RHODT  = 

0  L  +  p  L  + 

t>  L 

- 

p 

xx  y  y 

z  z 

c 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINE: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


READOBS 

SPIRDEC 


To  read  observations  from  logical  tape  0. 

JMP  READOBS 

OANDE  =  C/HLT,  OBLOC;  C/HLT,  EBLOC 
and  observations  on  logical  0. 

Observations  in  EBLOC 

SAVOBS  -  left  address  is  the  location  of  the  Z's 
terminating  the  observations 

System  -  SYS,  SYSNO,  SYSIO 


13  Cells 

This  subroutine  is  called  only  when  OUTOPT  =  1. 

OUTOPT  is  the  output  option  from  the  SPSJOB  Card . 

The  option  should  be  set  to  1  with  more  than  492 
but  not  more  than  984  observations  in  the  input. 

In  this  case,  the  Executive  routine  will  write 
all  the  observations  on  logical  0.  SPIRDEC  will 
read  them  in  only  once. 

Subroutine  READOBS  does  the  following: 

(1)  Rewind  logical  0 

EBLOC — >Index  Register  3 

(2)  Set  MCOUNT  =  0 

Read  1  block  into  the  address  at  location  T100; 
it  is  EBLOC  initially. 

(3)  Test  every  10th  word  for  Z's  in  the  block 
just  read: 

(a)  if  no  Z's  are  found,  add  10  to  the  address 
at  TI00,  add  10  to  Index  Register  3,  add  1 
to  MCOUNT: 

if  MCOUNT  =  12,  go  to  step  (2) 

if  MCOUNT  <  12,  continue  searching  for  Z's. 

(b)  if  Z's  are  found,  save  the  location  of  Z's 
in  the  left  address  of  SAVOBS  and  exit. 


192 


REJECI1 

SPIRDEC 


PURPOSE 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 

SUBROUTINES: 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


To  accept  or  reject  a  residual  for  ran^e 
for  1st  pass.  ° 


or  angles 


TMA  (residual) 
JMP  REJECT 1 
JAZ  (rejected) 


ABSMX  = 
SIGMAI  = 


rejection  criteria  (e.r.)  for  range  and  angles 
weight  corresponding  to  the  observed  quantity 


RCNT  =  C/HLT,  N;  C/HLT,  0 
SUM  =  F.P. 

REJCNT  =  C/HLT,  N;  C/HLT,  0 


Program:  RESREJ1 


b  Cells 


ifABS^ab^1Utvo^1Ue  °f  the  unwei8hted  residual 
-  ABSMX,  then  RCNT  is  increased  by  1,  and  the  weivht-prl 
residual  squared  is  added  to  SUM. 

,the"  R™T  “d  S»  «•  affected. 
RETFLT  (rejected  count)  is  increased  by  1,  the 

to'indicat^re jec tion^iSkS *  ^  **  (A)  is  c^red 
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REJECT2 

SPIRDEC 


PURPOSE:  To  accept  or  reject  a  range  rate  residual  for  1st  pass. 

CALL  SEQUENCE:  TMA  (residual) 

JMP  REJECT2 
JAZ  (rejected) 

INPUT:  ABMX2  =  rejection  criteria  (e.r./kemin) 

SIGMAI  =  a. 

o 

OUTPUT:  RCNT  =  C/HLT ,  N;  C/HLT,  0 

SUM  =  Floating  point 

REJCNT  =  C/HLT,  N;  C/HLT,  0 

SUBROUTINES:  Program:  RESREJ1 

STORAGE 

REQUIREMENTS:  6  Cells 

DESCRIPTION:  If  | unweighted  residual)  <  ABMX2,  then  RCNT  is  increased 

by  1  and  the  weighted  residual  (squared)  is  added  to  SUM. 

If  | unweighted  residual)  >  ABMX2,  then  RCNT  and  SUM  are 
not  affected,  but  the  REJCNT  (rejection  count)  is  increased 
by  1,  the  REJFLG  is  set  to  astericks,  and  the  (A)  register 
is  cleared  to  indicate  rejection. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


RESICK 

SPIRDEC 


To  restore  elements  to  continue  integration. 
JMP  RESICK 

BOO  +  16  -  BOO  +  55 

Contains  the  last  5  sets  of  elements 

ICK  +  0  -  ICK  +  39 

Contains  the  same  element  sets 

None 


4  cells 

To  continue  integration  for  a  prediction  or  the 
new  epoch  elements,  the  ICK  buffer  must  be  restored 
for  the  DIVDIF  subroutine. 
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RES OUT 1 
SPIRDEC 
1  of  2 


PURPOSE:  To  output  observation  residuals. 

CALL  SEQUENCE:  JMP  RES0UT1 

INPUT:  OBFLG  =  flag  for  observed  quantities 

RHOC  =  Pc 
RGSDL  =  Ap 
ARSDL  =  AA 
DRSDL  =  Ah 
RRSDL  =  Ap 

OUTPUT:  P,  a  &  6  or  A  &  h,  P  residuals 

These  residuals  flagged  with  asterisk  if  rejected 
Also,  VMAG,  DELU,  U,  BETA,  Time,  Station  Number 

SUBROUTINES:  Program  -  DELTAU,  RES0UT2 

System  -  GLOP,  PAGECON 

STORAGE 

REQUIREMENTS:  80  Cells 

DESCRIPTION:  If  only  range  rate  was  observed,  VMAG  is  not  computed 

otherwise 

VMAG  ="^(RGSDL)2  +  (ARSDL)2  +  (DRSDL)2 

If  angles  were  not  observed,  U,  At,  and  BETA  are  not 
computed,  otherwise: 

XCZS  =  p  L  -  X  =  x 

C  X  o 


YOBS  -  o  L  -  Y  =  y 
c  y  J  o 


ZOBS 

OBSR 
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RESOUT 1 
SPIRDEC 
2  of  2 


DESCRIPTION: 

(continued) 


RESOUTI  uses  subroutine  DELTAU  to  compute  Au. 
Then  computes  U: 


U 


SINU 

COSU 


+ 


DELU 


J 


Computes  At  (DELU) 


A  t  =  (-  Au/  V  P)  (r2  /  k  ) 

e 

BETA  =  sin-1  ^  (xo  Wx  +  yQ  Wy  +  ^  vy 
where  -90°  <  BETA  <  90° 

Tests  RESOPT: 

if  RESOPT  =  0,  then  print  angles  in  deg. 
if  RESOPT  =  1,  then  print  angles  in  km. 
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praposE- 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 


SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


RESOUT2  & 

RESOUT4 

SPIRUEC 


To  convert  time  to  BCD  and  output  page  heading  if 
necessary . 

JMP  RES0UT2 

T  =  t  (Min.) 

STAID  =  0  OXXX,  X  =  Sta.  No. 

OBYEAR  N 
OBMO 

OBDAY  \  BCD  time 

FOBHR 

FOBMIN 

FOBSEC  j 

STAID  -  Left  Justified 

Program  -  BCDTIM,  PAC-ECON 
System  -  PANT 


43  Cells 

RESOUT2  converts  time  to  BCD,  tests  if  new  page  headings 
are  needed.  If  so,  it  calls  RES0UT4  to  print  page 
headings  . 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


RESREJ 1 
SPIRDEC 


-a.  w 


_  _ 

passes  except  the  first 


TMQ 

TMA 

JMP 

JAZ 


(weighted  residual) 
(residual) 

RESREJ 1 

(residual  rejected) 


MAX  =  weighted  rejection  criteria 

ASTK  =  w /******** 


RCNT 

SUM 

REJCNT  = 
REJFLG  = 

(None) 


C/HLT,  N 
Floating  point 
C/HLT,  N 
*  or  blank 


7  Cells 

If  |  weighted  residual  |  >  MAX  then  rf  tpht  ,• 

REJFLG  =  *,  and  (A)  reg.  =  0  18  increased  bX  1, 

If  |  weighted  residual j  <  MAX,  then  RCNT  is  increased  bv  1 
(weighted  residual)  is  added  to  SUM,  and  REJFLG  =  blank^  ’ 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


SUBROUTINES: 

STORAGE 

REQUIREMENTS 

DESCRIPTION: 


RESW 

SPIRDEC 


To  restore  values  to  the  Adams-Bashf or th  buffer. 


JMP  RESW 


QQ  +  1  =  L 

2  =  a 

x 

3  =  a 

y 

A  =  a 

z 

5  =  h 

x 

6  =  h 

y 

7  =  h 

z 


W  +  5  =  L 

6  =  a 

x 

7  =  a 

y 

8  =  a 

z 

9  =  h 

x 

10  =  h 

y 

11  =  h 

z 


Program  -  SETW 


2  Cells 

Restores  values  saved  in  QQ  +  1  to  QQ  +  /  to  W  +  5  to  W  +  11 
to  continue  the  Adams-Bashfor th  integration. 
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RESWBF 

SPIRDEC 


PURPOSE : 

To  restore  W  buffer 

CALL  SEQUENCE: 

JMP  RESWBF 

INPUT: 

PREDBF  buffer  (11  cells) 

OUTPUT: 

W  t  1  to  W  +  11  restored  from  PREDBF  buffer  (11  cells) 

SUBROUTINES : 

Program  -  SAVWBF 

STORAGE 
REQUIREMENTS : 

3  Cells 

DESCRIPTION: 

Restores  W  buffer  to  last  point  in  ADBASH  subroutine 
m  order  to  restart  integration. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT : 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


REVSUB 

SP1RDEC 


To  update  the  revolution  number. 

JMP  REVSUB 

UZ  =  U  at  t 

z 

OLDUZ  =  U  at  t  -  It 
z 

OLDUZ  =  U  at  t 

7. 

REV  =  Rev  No  at  t 

None 


5  Cells 

A  switch  CNTSW  must  be  preset  before  calling  the 
subroutine: 

If  integrating  forward  from  epoch: 

CNTSW  =  C/JAN,  CNTB;  C/TMA,  F/I 

If  integrating  backward: 

CNTSW  =  C/JAP,  CNTB;  C/TMA,  F/-1 

If  the  signs  of  UZ  and  OLDUZ  are  the  same, 
the  rev.  no.  is  not  modified. 

If  the  sign«  are  different,  then  a  node  has  been 
crossed.  If  it  is  the  ascending  node,  the  rev.  no, 
modi f  red  by  1 
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RHOSB 

SPIRDEC 


PURPOSE:  To  compute  range 

CALL  SEQUENCE:  JMP  RHOSB 

INPUT:  X  =  x 

Y  =  y 
Z  =  z 
CAPX  =  X 
CAPY  =  Y 
CAPZ  =  Z 

OUTPUT:  RHOC  =  p 

RHOX  =  o£ 

RHOY  =  p 

RHOZ  =  oz 

SUBROUTINES:  Philco  -  FSQRT 

STORAGE 

REQUIREMENTS:  LL  Cells 

DESCRIPTION:  RHOX  =  X-H3APX  (ie.p  =  x  +  X) 

x  4 

RHOY  =  Y+CAPY 
RHOZ  =  Z+CAPZ 

RHOC  =  pc  =  '/(RHOX)2  +  (RHOY)2  +  (RHOZ)2 
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i 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 


SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


RINEL 

SPIRDEC 


To  restore  initial  elements  to  output  cells. 
jMP  RINEL 


INELT  buffer 


LPRINT  = 
AX PRINT  = 
AY  PR  I  NT  = 
HX PRINT  = 
HYPRINT  = 
HZ PR I NT  = 
B PRINT  = 
HSUBQP  = 
PAPRINT  = 


L  (deg) 
a 

xno 

a 

yno 

h 

x 

h 

y 

h 

z 


H  (km) 
period  (min.) 


Program  -  SINEL 


2  Cells 

Moves  value  from  INELT  buffer  to  output  cells. 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 


SUBROUTINES: 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


RR2AHL 
SPIRDEC 
1  of  2 


To  convert  r,  r  elements  to  a,  h,  L 

JMP  RR2AHL 
(return) 

X  =  x  (km) 

Y  =  y  (km) 

Z  =  z  (km) 

XDOT  =  x  (km/ sec) 

YDOT  =  y  (km/sec) 

ZDOT  =  z  (km/ sec) 

HXO  =  h 

xo 

HYO  =  h 

yo 

HZO  =  h 

zo 

AXNO  =  a 

xno 

AYNO  *  a 

yno 

XLO  =  L 

o 

+  cells  on  next  page 

Philco  -  FSQRT 
System  -  ARCTAN 


94  Cells 

Converts  from  r,  r  to  a,  h,  L  using  the  following 
formulation. 

Called  by  subroutine  INPUT  only. 
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RR2AHL 

SF"Rr£C 

2  ci  2 


RR2AHL  Formulation 


HXO  =  h 

xo 

=  y 

z  -  z 

y 

HYO  =  h 

-  z 

X  -  X 

z 

yo 

HZO  =  h 

=  X 

y  -  y 

X 

zo 

P  =  h  2  + 

h 

2  +  h 

2 

xo 

yo  zo 

=  P 


RTP  =  \  p 

WX  =  W  =  HXO  /  ",  p 

R  =  ^ x2  +  y2  +  z2  =  r 
UX  =  U  -  x/r. 


— »y 


x  — -,V  ,  z 


RDOT  = 

r  =  (xx 

+ 

yv 

+  ZZ) 

It 

ES1NEV 

=  r  A  P 

= 

e  sin 

V 

ECOSV  = 

p/r  -1 

= 

e 

cos  V 

ESQ  =  (e  sin  v) 

2 

+ 

(e  cos 

■> 

v)fc 

EO  =  e 

=  ‘  e2 

o 

o 

VX  =  v 

X 

=  (rx  - 

xr) 

2*  , 

'  .xf 

)  X 

n 

cd 

II 

< 

p/l-e2 

LLO  =  1 

=  tan 

1 

-  Vx 

=  e 


0  <  1  <  360 

-  o  — 


►y  > z 


XNO  =  N  =  K  /  a3/2 

AX  =  a  =  U  e  cos  v  -  V  e  sin  v  , 
xx  x 

AXNO  =  a  =  (a  W  -  a  W  )  /  \ 1-W  2 

xno  y  x  x  y  z 

AYNO  =  a  VT-Wz2  •  Wz  a  W  +  a  .  W  )  /  1-W  2 

7.  V  V  H  tr'  '  _ 


v  -  E  =  tan 


x  x  y  y ' 

i  .  2 

i  e  sin  v  ■.  1  +  1-e  +  ecosv) 


|  2  2 
j  * 'e  sin  v)  +  cos  v  +  1)  (1  +  1-e  ) 


x  — i-A  ( z 


QQ+  2  -v-M-v-E+re  sinv  /1-e2  /p 


XLO  =  L  =  1  -  <’v  -M ) 

o  o 
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SAVCON 

SPIRDEC 


PURPOSE : 

CALL  SEQUENCE: 


INPUT: 


OUTPUT: 

SUBROUTINES: 


To  save  card  images  of  the  parameter  cards 

TMD  C/HLT,  SCONBUF+X 
JMP  SAVCON 


Index  register  3  =  1st 
card  image  in  C0N3UF 


location  of  a  parameter 


See  description. 


None 


STORAGE 
REQUIREMENTS : 


DESCRIPTION: 


The  parameter  card  image  is  moved  to  th 
CONBUF,  to  be  used  by  subroutine  E» 
output  purposes . 


buffer, 

for 


SCONBUF  +  0  to  +  9  P  Card  2 

+  10  to  +  19  P  card  3 

+  20  to  +  29  P  card  5 

+  30  to  +  39  ?  card  6 

+  40  to  +  49  P  card  7 

+  50  to  +  59  P  Card  8 
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SPTKDt: 


PURPOSE 

CALL  SEQ 
INPUT: 

QUIP JT 

subrout: 

STORAGE 

REQUIREM 

DS3CR1PT 


'  '  r  •  i.  -  :  -CrSio'V  *  0  CO  Z  l  Cl,  e 

i  v.r.e.'c  :io~  co  've  new  epoch  e  l ament  3 
JEN1,  r  Jyy 

5  aeco  o.  c..^v.i  in  the  ICK  t 
REV  *-  ■  ;',e  _  in  number 

GLLTi  -  i  .  -he  piTvio.s  cite 

W  —  aL'B.V'S  o  . 

See  desc.t  ■>  o~ 

N’ES.  Program  -  S.-n-WSE'.  SAVICK 

ENTS;  5  Tells 

TON.  This  subroutine  only  called  by  the  D.C<  control 

subroutine  (CMRL)  _f  the  new  epoch  t^me  requested 
is  outside  the  scan  of  the  observations  and  on  the 
.pc  site  6ice  :i  the  old  epoch  from  the  last 
obssi' v-aticn .  A  third  requirement  fox  ~3ing  this 
rot  tine  is  that  the  old  epoch  must  be  within  the 
span  cf  the  observations, 

'h~  autre u: » re  v;  .. 

'  ye  SAVWo?  to  save  W  +  1— >W  +  11 

se  oA  ICK  to  save  ICR  — ^ICE  d  3'- 

rive  :.r  SR£V 

'  1  S  OLD  HZ 

;  r  Sr. Vr. i  -  ■  d  .  .  .  .  f-  t  ••  jrf 

W?  d  'i  te<.! 

”a  'ues  s-  <.<  i  »  t  - 1’  •  R ' 


tier 
pc  in  t. 
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PURPOSE: 

CALL  SEQUENCE: 

INPUT: 


OUTPUT: 


SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


SAVICK 

SPIRDEC 


To  save  values  to  compute  the  new  epoch  elements. 

TMD  C/HLT,  ICK+X;  C/HLT,  B00f8 
JMP  SAVICK 

AXN  =  a 

xn 

AYN  =  a 

yn 

B  =  CD  A/m 
A  =  a 
ESQ  =  e2 
SINI  =  sin  i 
P  =  p 
XN  =  n 

ICK  buffer  from  DIVDIF  subroutine 
BOO  +  0  =  a 

xn 

+  1  =  a 

yn 

+  2  =  B 

+  3  =  a 

+  4  =  e2 

+  5  =  sin  i 

+  6  =  p 

+  7  =  n 

+  8  thru  +15  =  new  epoch  elements 

+  16  thru  +55  =  last  5  sets  of  integrated  elements 

Program-  SAW,  SEW,  SUBXYZ,  RESW 

14  Cells 

Saves  the  new  epoch  elements  in  BOO  +  8  thru  BOO  +  15. 
Values  saved  in  BOO  +  0  thru  BOO  +  7  are  necessary  to 
prepare  the  new  epoch  elements  for  output.  In  BOO  +  lb 
thru  BOO  +  55  the  last  5  sets  of  elements  are  saved  to 
continue  integration  if  a  prediction  is  desired. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


SAVW 

SPIRDEC 


To  save  the  elements  in  the  Adams-Bashf or th  buffer. 


JMP  SAW 


W  +  5  =  L 

6  =  a 

x 

7  =  a 

y 

8  =  a 

z 

9  =  h 

x 

10  =  h 

y 

11  =  h 

z 


QQ  +  1  -  L 

2  =  a 

x 

3  =  a 

y 

4  =  a 

z 

5  =  h 

x 

6  =  h 

y 

7  =  h 

z 


Program  -  SETW 


2  Cel  Is 


Moves  values  from  W  +  5  thru  W  +  11  to  QQ  +  1  thru  QQ  +  7 
so  that  other  values  may  be  placed  in  W  +  5  thru  W  +  11 
to  be  able  to  use  the  SUBXYZ  subroutine. 
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SAVWBF 

SPIRDEC 


PURPOSE : 

To  save  W  buffer 

CALL  SEQUENCE: 

JMP  SAVWBF 

INPUT: 

W  buffer 

OUTPUT: 

W  +  1  to  W  +  11 - *PREDBF  +0  to  PREDBF  +  10 

SUBROUTINES: 

(None) 

STORAGE 
REQUIREMENTS : 

5  Cells 

DESCRIPTION: 

Saves  11  cells  of  the  W  buffer  in  order  to  restart 
ADBASH  subroutine. 
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SAV5PTS 

SPIRDEC 


PURPOSE : 


To  save  elements  for  interpolation 


CALL  SEQUENCE: 


JMP  SAV5PTS 


INPUT: 


W  +  1  =  t  (min) 

5  -  L 


6  =  a 


8  =  a 

9  =  h 

> 

10  =  h 

11  =  h 


OUTPUT: 


See  description. 

ICK  +  0  =  t 

+  1  =  L 

+  2  =  a 


+  39 


5  sets 


SUBROUTINES: 


None 


STORAGE 
REQUIREMENTS : 


8  Cells 


DESCRIPTION: 


Buffer  ICK  contains  5  sets  of  elements  at  t^  -  t<.  . 
When  t |  is  no  longer  needed,  elements  at  t„  -  t,.  are 
moved  to  replaced  t.  -  t. .  And  the  elements  at  t, 

]-  Q 

in  the  W  buffer  are  moved  to  replace  the  elements  at 
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SENLOC 

SPIRDEC 


PURPOSE: 

CALL  SEQUENCE: 


To  compute  6,  sin  9,  cos  6 


TMA  T 
JMP  SENLOC 


INPUT: 

T 

= 

t  -  minutes 

XLAMBA 

— 

\ 

E 

THGR 

= 

e 

gr 

RPTIM 

= 

.0043752691 

OUTPUT: 

THTA 

= 

9  at  time  t 

SI  NTH 

= 

sin  9 

COSTH 

= 

cos  9 

SUBROUTINES : 

Philco  - 

FCOS,  FSIN 

STORAGE 

REQUIREMENTS: 

6  Cells 

DESCRIPTION: 

THTA  = 

e 

=  t  (.00437 

gr 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINE: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


SETSBUF 
SPIRDEC 
1  of  2 


To  set  up  a  modified  SBLOC  in  SBUF. 
JMP  SETSBUF 

BIASAD  =  C/HLT,  EBLOC ;  C /HLT ,  BIBUF 
Observations  starting  in  EBLOC 
Sensors  in  SBLOC 

Modified  SBLOC  in  buffer  SBUF 

See  buffer  layout  under  description. 

System  -  SGET 


19  Cells 

Starting  with  the  first  observation,  extract  the 
sensor  number: 

(1)  if  the  sensor  information  is  in  SBUF,  go  to 
the  next  observation 

(2)  if  the  sensor  is  not  in  SBUF,  call  subroutine 
SGET  to  unpack  the  entry  in  SBLOC,  and  store 
only  5  quantities  for  the  sensor  in  SBUF,  then 
go  to  the  next  observation 

(3)  if  SBUF  is  full  (30  sensors) ,  exit  without  error 
indication.  The  retrieval  routine  GETSEN  has  an 
error  exit. 

The  sensor  information  must  be  moved,  since  SBLOC  is 
used  for  a  working  buffer  and  for  temporary  storage. 
SBUF  is  the  same  buffer  as  LIBUF;  BIBUF  (bias  buffer) 
is  no  longer  needed  at  this  point. 
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SETSBiF 
SPIRDEC 
2  of  2 


DESCRIPTION: 

SBUF  Format 

(continued) 

SBUF  +  0 

00000SSS 

Sensor  Number 

1 

m  (radians) 

PHIRD 

2 

X  (radians) 

XLAMBA 

3 

x  /  cos  Q 

XOVCT 

4 

Z 

CAPZ 

1 

1 

1 

l 

1 

145 

OOOOOSSS 

146 

CD 

147 

X 

148 

x  /  cos  6 

149 

Z 

150 

OOOOOZZZ 

5  words  /  entry  -  to  a  maximum  of  30  sensors  - 
by  00000ZZZ  after  the  last  entry. 


terminated 
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SETW 

SPIRDEC 


PURPOSE:  To  move  interpolated  elements  at  time  t  to  the 

Adams-Bashf or th  buffer. 

CALL  SEQUENCE:  JMP  SETW 

INPUT:  ICK  +  40  -  L 

4  L 

42 

43 

44 

43 

44 

OUTPUT:  W  +  5  =  L 

6  =  a 

x 

7  =  a 

y 

8  -  a 

z 

9  =  h 

x 

10  -  h 

y 

11  =  h 

z 

SUBROUTINES:  None. 

STORAGE 

REQUIREMENTS:  4  Cells 

DESCRIPTION:  ICK  +  40  thru  ICK  +47  contains  elements  at  time  t 

computed  by  the  DIVDIF  subroutine.  To  compute  the 
x,  y,  z  position  at  this  time,  these  values  must  be 
ia  W  +  3  thru  W  +  11  to  use  the  SIIRXYZ  subroutine. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS 

DESCRIPTION: 


SINEL 

SPIRDEC 


To  save  initial  elements  for  output. 


JMP  SINEL 


LPRINT  = 
AX PRINT  = 
AY PRINT  = 
HXPRINT  = 
HYPRINT  = 
HZPRINT  = 
B PRINT  = 
HSUBQP  = 
PAPRINT  = 


Lq  (deg) 
a 

xno 

a 

yno 

h 

x  o 
h 

yo 

h 

z  o 

CpA/m 

perigee  altitude  (km) 
period  (minutes) 


Input  values  in  buffer 
INELT  -  INELT+8 


None 


4  Cells 

Moves  input  values  to  the  INELT  buffer  in  the 
order  indicated. 
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SORTOB 

SPIRDiX 


PURPOSE: 

To  sort  the 

observations  by  time. 

CALL  SEQUENCE: 

JMP  SORTOB 

INPUT: 

MCOUNT 

C/HLT,  0;  C/HLT,  N 

N  -  nunbe-  of  negaci Te  observation  times 

OANDE 

C/HLT,  O3L0C;  C'HLT,  BLOC 

OUTPUT: 

See  description 

SUBROUTINES: 

None 

STORAGE 

REQUIREMENTS: 

33  Cells 

DESCRIPTION: 

(1)  Set  the  right  aduress  of  TEMPI  and  TEMP 2  to  the 
number  of  observations  minus  one. 

Set  switch  SORTOB3  to  :nrt  the  conservations  in 
true  time  order. 

(2)  EBLOC — 5>Incax  Register  4 
EBLOC  +  10  — >Iudex  Register  5 
C/HLT,  OBLOC;  C  /:  ILT ,  EBLOC  — >  QQ 

(3)  PREDBF  — > Index  Register  6 

(4)  Sort  observations  in  time  order  -  make  as  many 
passes  through  the  observations  as  the  number  in 
TEMPI,  using  TE'D'2  to  tell  when  one  pass  is  complete. 
TEMP2  vik.1  match  TEMPI  at  tie  start  of  each  pass. 

(5)  Test  MCOUNT: 

(a)  if  =  0,  exit 

(b)  if  ^  fl,  then  there  are  some  observations 
before  epoch,  (i.e.  some  negative  times). 

‘These  must  he  resorted  in  reverse  order. 

Example:  buffer  order  would  be 

-2,  . -100,  2  _ 100 

since  observations  will  be 
retrieved  in  this  order. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


SUBOl'T 

SPIRDEC 


To  output  the  prediction  ephemeris. 

JMP  SUBOUT 
NEWPAGE  = 

PFLAG 
X 
Y 
Z 

XDOT 
YDOT 
ZDOT 

X  ~ 

HEDLIN 
HEDLINI 

> 

VALLIN 
VALLINI 

According  to  value  of  PFLAG.  (See  description.) 

System  -  GLOP,  PANT 

Program  -  PAGECON,  PHLAH,  BCDTIM 

65  cells 

If  NEWPAGE  =  0,  then  outputs  headings  at  the  top  of 
a  new  page,  using  locations  HEDLIN  and  HEDLINI. 

If  PFLAG  =  "  •  3,  subroutine  PHLAH  is  called  to 

compute  0,  h. 

If  PFLAG  =  1  or  3 ,  then  x,  y,  z,  x,  y,  z,  are 
converted  to  output  units  of  km  and  km/sec. 


page  indicator 
print  option 
x 

y 

z 

X 

y 

z 

t  (min) 

from  SUBOUTI 


The  values  requested  by  PFLAG  are  then  printed 
and  PAGECON  is  called  to  update  the  line  count. 


PURPOSE : 

CALL  SEQUENCE: 

INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 
REQUIREMENTS : 

DESCRIPTION: 


SUBOUTI 

SPIRDEC 


To  initialize  subroutine  SUBOUT. 

JMP  SUBOUTI 

PFLAG  =  0,  1,  2,  or  3  T47  (print  option  :rom  P  card  4) 

See  description. 

Program  -  PAGECON 


14  cells 

Sets  up  parameters  for  SUBOUT  (HEDLIN,  HEDLINI, 
VALLIN,  VALLINI)  according  to  IT  LAG. 

Then  use  PAGECON  to  force  a  page  and  sets 
NEWPAGE  =  0. 

If  PFLAG  =  0  or  2: 

HEDLIN  =  C/TMA,  LINE4;  C/TIJ,  L:  NE43 

HEDLINI  =  0  -  0 

VALLIN  =  C/TMA,  P0C7D;  C/TIJ,  HEDL23 
VALLINI  ■-  0  -  0 

If  LAG  =  1 

HEDLIN  =  C/TMA,  LINE4 ;  C/TIJ,  LINE4A 

HEDLINI  =  0 - 0 

VALLIN  =  C/TMA,  POC7D ;  C/TIJ,  HEDL21 
VALLINI  =  0 - 0 

And  sets  switch  SUBOUT2  to  jump  to  SUB0UT7 . 

If  PFLAG  =  3 

HEDLIN  =  C/TMA,  LINE4 ;  "/TIJ,  LINE4B 
HEDLINI  =  C/HLT,  0;  C/TIJ,  LINE4A 
VALLIN  =  C/TMA,  P0C7D;  C/TIJ,  HEDL23 
VALLINI  =  C/HLT,  0;  C/TIJ,  HEDL21 
And  sets  switch  SUB0UT2  to  jump  to  SUB0UT8. 
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S'TBXYZ 
S  PIRDEC 
1  of  3 


PURPOSE:  To  compute  _r,  jr  from  A,  h,  L 

CALL  SEQUENCE:  JMP  SUBXYZ 
JMP  (ERROR) 

JMP  NORMAL 

INPUT:  W  +  5  =  L 

W  +  6  =  a 

x 

W  +  7  =  a 

y 

W  +  8  =  a 

z 

W  +  9  =  h 

x 

W  +  10=  h 

y 

W  +  11*  h 

z 

OUTPUT:  X  =  x 

Y  =  y 
Z  *  z 
XDOT  =  x 
YD OT  =  y 
ZDOT  =  z 

plus  necessary  intermedia r.e  quantities  listed  on  following 
pages . 

SUBROUTINES:  Program  -  ARCTAN 


STORAGE 

REQUIREMENTS.  100  cells 
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SUBXYZ 

SPIRDEC 


Kepler's  Equation: 
(1) 


COSEO  =  cos  (E  +  cjo  ) 

SINEO  =  sin  (E  +  a'  ) 

ECOSE  =  SINEO  (a  )  +  COSEO  (a  ) 
yn  xn 


Test: 


QQ  =  1  -  ECOSE 

ESINE  =  SINEO  (a  )  -  COSEO 
xn 


EPSI  = 


ESINE  +  U  -  EOI  i 
1  -•  ECOSE  i 


) 


if 


<  €  , 


iterate  again  to  Eq .  (1)  after 
replacing  to  a  maximum  of  50  times,  1 
then  take  error  exit  from  sub- 


EOl  by 
EOI  -  EPSI 


routine 

i 

if  >  t  , 

continue 

R  =  A  (1 

-  ECOSE)  =  r 

RDOT  =  ( 

A '  a 

V  p  /  r) 

ESINE  = 

r 

RVDOT  =  ( 

\  a  A  P  /  r) 

,  1-e2  = 

:  rv 

COSU  =  ESINE 

(1-  1-e2 

a 

-a  +  COSEO  a/r 

1 

yn 

xn 

SINU  =  SINEO 

-  a  -  a 

ESINE 

(1-  / 1-e2  a/ r 

yn  xn 

AR  =  a/r 

UX  =  COSU 

N  +  SINU  ‘ 

M 

X 

X 

UY  =  COSU 

'  N 

!  +  SINU  ’ 

M 

y 

y 

UZ  =  COSU 

'  N 

+  SINU  ‘ 

M 

z 

z 

VX  =  COSU 

'  M 

-  SINU  • 

N 

X 

X 

VY  =  COSU 

'  M 

-  SINU  • 

N 

y 

y 

VZ  =  COSU 

'  M 

-  SINU  • 

N 

z 


z 
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Sl'BXYZ 
SPIRDEC 
3  of  3 


X  =  r  (  Ux  )  =  x 

Y  =  r  (  Uy  )  =  y 

Z  =  r  (  U  )  =  z 

z 

XDOT  =  r  TJ  +  r  v  V  =  x 
x  x 

YDOT  =•  r  U  +  r  v  V  =  y 

y  y 

ZDOT  =  r  lT  +  r  v  V  =  z 
z  z 


223 


PURPOSE: 

CALL  SEQUENCE: 
INPUT  : 


OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


TAPEW 

SPIRDEC 


To  write  a  binary  ephemeris  tape  on  logical  1 

JMF  TAPEW 

W+l  *  t  (min) 

X  (e.r.) 

Y  (e.r.) 

Z  (e.r.) 

XDOT  (e.r./kemin) 

YDOT  (e.r./kemin) 

ZDOT  (e.r./kemin) 

TAPCNT  =  C/HLT,  TAPBUF  +  X 

One  block  on  the  binary  ephemeris  tape  for  XYZLA 
subroutine  when  buffer  is  full. 

System  -  SYS,  SYSNO,  SYS 10 


16  Cells 

Adds  t,  x,  y,  z,  x,  y,  z  to  tape  output  buffer  (TAPBUF). 
When  the  buffer  is  full,  one  block  (128  words)  is 
written  on  tape.  TAPCNT  is  updated. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINES : 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


TEMP 

SPIRDEC 


To  compute  temperature  at  epoch. 

JMP  TEMP 

FIOAV  =  F10 

FLO  =  F10 

AP  '  Ap 

ORGDAY  =  days  from  beginning  of  year  to  epoch =  D 
TEMPO  =  Tq  (°K) 

None 


19 

cells 

To 

=  974°  + 

4.2° 

<F10  ■ 

-  150)  + 

0.004° 

<F10 

-  150)2 

T' 

Ao 

=  T  + 
io 

1.9° 

(F10  ■ 

f’lo)  - 

r 

f 

\ 

!0.39°+0 

f  *1 

T 

lo 

=  Ti  + 

F~10 

sin  <4n 

(D-60); 

.15° 

sin<2^  (D- 150 

1 

365  •  ' 

l  365  J 
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THGRC 

SPIRDEC 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 


OUTPUT: 


SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


To  compute  0  at  epoch 
JMP  THGRC 

-  0  OY ,  Y  =  BCD  year 

ORGTM  =  fraction  of  epoch  day 
THGRO  =  0  from  TLC  subroutine 

ORGDAY  =  days  since  beginning  of  year 

TO  -  days  and  fraction  since  beginning  of  year 
THGR  =  0 

ORGTM  =  minutes  in  epoch  day 
Program  -  TLC,  MOD2PI 

13  Cells 

THGR  =  9gr+  (Days  since  beginning  of  year)  (0.9856473354) 
(fraction  of  day)  (360.9856473) 

Converts  ORGTM  from  fraction  of  day  to  minutes 


+ 
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WF  TC-H  r 
spifdei 

1  of  2 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


To  read  the  weight  tape  and  set  up  the 
weight  and  bias  buffers. 

JMP  WEIGHT 

Weights  and  biases  on  logical  tape  7 
WAND B I  =  C/HLT,  WBUF;  C/HLT,  BIBUF 
NXTWGTA  =  C/HLT,  ABIT  +120;  C/HLT,  0 

See  description 

System  -  PANT,  SYS,  SYSNO,  SYSIO 
Program  -  XSRCH 


82  Cells 

This  subroutine  is  called  if  WGTFLG  =  1 
(Column  2  on  P  Card  3) ,  which  indicates  that 
a  weight  tape  is  mounted  and  should  be  read. 

The  following  procedure  is  executed: 

(1)  Read  a  block  from  logical  tape  7  initially, 
then  only  read  a  new  block  when  all  cards 
in  the  block  have  been  processed: 

(a)  if  the  first  word  of  a  card  is  ENDSIGMA, 
the  tape  is  complete,  rewind  logical  7 
and  exit . 

(b)  if  not  ENDSIGMA,  unpack  one  card  using 
XSRCH 

(1)  Set  WCTBUF  =  OOOOOSSS 

where  SSS  =  Sensor  Number 

(1/.  Save  the  address  of  the  next 
card  in  NXTWGTA. 
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DESCRIPTION: 

(continued) 


WEIGHT 
SPIRDEC 
2  of  2 


(2)  Search  WBUF  to  find  a  word  of  G0000ZZZ  or 
the  same  Sensor  Number. 

(a)  if  the  Sensor  Number  is  found,  replace 
the  assembled  values  with  the  new  values 
after  conversion  as  fellows: 


WBJF 

(Weight  Buffer) 

BIBJF 

(Bias  buffer) 

Word  0 

ooooosss 

Word  0 

OOOOOSSS 

1 

l/tfp  (e.r.) 

-Bp  'e.r.) 

2 

i/<?A  (rad) 

O 

-Ba  (rad) 

3 

U\  (rad) 

3 

~Bh  (rad) 

4 

l/sp  (e.r/K) 

4 

“Bp  (e.r/K) 

5 

-B  (min) 

rr-> 

Go  to  Step  (1) . 

(b)  if  the  sensor  number  is  not  found,  but 
00000ZZZ  is  found  check  to  see  if  there 
is  room  for  another  entry 

(1)  if  there  is  no  room  (i.e,  30  sensors 
have  been  entered) ,  print  a  comment 
and  go  to  step  (l) ,  skipping  the  entry 

(L)  if  there  is  room, follow  same  procedure 

of  storing  and  computing  values  at  (2)  (a) , 
also  add  a  word  of  OOOOOZZZ  at  the  end 
of  each  buffer. 
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W  ILE  JR 
SFIRDEC 


PURPOSE: 

CALL  SEQUENCE: 
INPUT: 

OUTPUT: 

SUBROUTINES: 

STORAGE 

REQUIREMENTS: 

DESCRIPTION: 


To  print  the  input  observations. 

JMP  WILBUR 

Observations  starting  in  location  EBLOC , 
terminated  by  a  word  of  Z’^. 

Observations  printed  in  card  format  on  bard  copy. 

Program  -  PAGECON,  SEPSUB 
System  -  PANT,  GLOP,  DKLOK 


164  Cells 

This  subroutine  is  called  by  option,  if  Column  l  on 
P  Card  3=1.  The  observations  are  in  OBLOC  format 
and  must  be  converted  to  card  format  for  output. 

The  observations  will  be  printed  in  the  same  order 
as  they  are  in  the  input  deck,  as  they  have  not  yet 
been  sorted.  In  addition,  the  column  numbers  will 
be  printed  at  the  top  of  each  page.  A  space  will 
appear  between  each  field  of  the  observation  card 
for  clarity. 
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PURPOSE: 

CALL  SEQUENCE: 
INPUT. 


OUTPUT: 


SUBROUTINES : 

SIORAGE 

REQUIREMENTS: 

DESCRIPTION: 


WSEI.P 

SPIP.DEC 


To  set  up  W  buffer  for  Adams -Bashf orth  integration. 
JMP  WSETuP 


XLO 

= 

L 

o 

AXO 

= 

a 

xo 

AYO 

= 

a 

yo 

AZO 

— 

a 

zo 

HXO 

= 

h 

xo 

HYO 

Vi 

yo 

HZO 

= 

h 

zo 

W  + 

1 

=  F/O 

w  + 

5 

=  XLO 

w  + 

6 

=  AXO 

w  + 

/ 

=  AYO 

w  + 

8 

=  AZO 

w  + 

9 

=  HXO 

W  +10 

=  HYO 

W  +11 

=  HZO 

(None) 

LO  Cel Ls 

Moves  lament  see  at  epoch,  i.e.  (W  +  L  =  T>  to 
W  buffer  for  subroutine  AD BASH . 
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APPENDIX 

ADAMS -BASHF OR TH  INIEGRAIxON" 


The  Adams-Bashf or th  (A-B)  method  together  with  the  Runge-Kutta  (R  -K) 
method  of  integration  are  employed  in  this  subroutine.  R-K  is  used  as  a  starting 
procedure  so  that  an  adequate  number  of  step-wise  solutions  may  be  obtained  to 
build  a  difference  table  needed  by  the  A-B  method.  The  interval  of  integration  is 
automatically  varied  to  keep  the  discrepancy  between  the  integrated  values  and  an 
absolute  error  check  within  prescribed  limits. 

1-1.  R-K  PROCEDURE  FOR  STARTING 


To  start  the  Adams-Bashf or th  procedure,  a  series  of  R-K  integration 
steps  (at  equal  intervals  of  X)  are  generated  until  enough  points  are  available 
to  set  up  the  difference  table  needed  by  A-B. 


The  system  of  equations  to  be  solved  is  given  by. 


Yi =  fi  V  V 


V  \ 

‘  7' 


Y.  (X  )  =  Y. 

1  c  1C 


l  1,  2 ,  ..*?7 


Let  Y  be  the  value  of  \  at  X-X  and  f  the  derivative  of  \  at  X=X  ^nd 
let  h  be  the  step  size  of  the  i  ndependent  variaBle  X  note  that  the  subscript.  ’’  i,! 
has  been  omitted  for  simpLK  it>  .  Ihe  R-K  method  uses  the  classical  iourth -order 
formulas: 


K,  =  hf  (X  ,  Y  ) 

1  r  n 

K  =  hf  (X  +  ■;  h,  Y  +  ~  K1  ) 
1  n  2  n  2  1 
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K  =  hf  (X  +  ±  h,  Y  +  \  K  ) 
3  n  L  n  2  2 


K  =  hf  (X  +  h,  Y  +  JO 
4  n  n  3' 

Y  .  =  Y  +  \  (K,  +  2K„  +  2K  +  K.) 
n+ 1  n  6  1  ^  3  4 


a.  Interval  Control  and  Modifications 


The  method  of  interval  control,  while  in  the  R-K  section,  consists  of 
generating,  from  the  point  X  ,  Y  and  its  derivative  f  ,  two  R-K  points  Y  , 

Y  , „  and  their  derivatives  f  f  With  this  information,  a  Simpson's  Rule 

integration  is  performed  on  the  integral: 

Ay  =  fn+2  f .  dx 

Jj-n  J 


By  Simpson's  Rule: 


(f  +  4  f  , .  +  f  ,0) 
n  n+ 1  n+2 


The  following  evaluation  is  made  for  each  of  the  seven  equations: 


e.  =  absolute  error  control  supplied  by  the  user  for  each  equation. 


If,  for  any  of  £h^  seven  equations,  the  value  ^  equals  or  _xceeds  0, 

then  h  is  reduced  by  10  '  ,  and  the  integration  is  restarted  from  X  ,  Y  .  This 

procedure  is  repeated  until  all  ft  are  less  than  0:  then  the  values  of^  T*,  Y 

„  ,  _  ‘  ,  .  ,  .  „  rn-2  ’  n+2 

are  returned  to  the  user  as  valid  points. 


b.  Difference  Table  Construction 


After  the  appropriate  number*  of  R-K  integration  points  have  been 
generate^,  and  seven  values  of  f  (j  =  0,  1,  ..,,6^  at  equal  steps  are  available, 
then  the  difference  table  requiredJby  A-B  is  constructed. 


*The  program  generates  4  R-K  points  per  1  A-B  point,  therefore,  25  R-K  points 
must  be  available  to  start  A-B. 


232 


233 


rhe  predictor  is  of  the  form 


6  r 

'  =  Y  +  h  r-  B  VJ 

n+1  n  r  Vn 


Bo=  1 


Bl  =  1/2 


B2  =  5/12 


B3  =  3/8 


B.  =  251/720 
4 


B^  =  95/288 


B6  =  19087/60480 


The  corrector  is  of  the  form 

6  .  r 
Y  ,  =  Y  +  h  v-  s  n  r 

p+1.  n  4_  r  Vn+1 

r=0 


b;  =  -i/2 

B^  =  -1/12 
-  -1/24 

— / 

B.'  =  -19/720 
4 

Bg  =  -3/160 
B^  =  -863/60480 
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The  diif^ren  es  VJ  r ,  .sad  in  the  correacor  are  competed  by 

Vn  1  0 
evaluating  the  derivative  cf  the  predictor  f^+^  ^  ;  then 

Vr  _  r- 1  r-1 

trt-L  ‘Vnfl  '  V"  ’  6 

-3.  A-B  INTERVAL  CONTROL  AND  MODIFICATION 

Aftar  the  corrector  has  been  computed,  the  quantity  6  Le,  evaluated. 
^  =  Jpredictor  -  corrector)  -  e 

e.  =  absolute  error  control 

i 

If,  for  any  equation,  £>0;  then  the  interval,  h,  will  be  decreased 
by  1/2.  New  differences,^  (r=i,  2,  ...6),  based  on  half  the  old 
interval  are  then  computed  by: 


Vn1 

vn2 
v„3 
v„4 
1  v„5 
v„6 


0  0 

0  0 

0  0 

0  0 


vnl 


.125  .0625 

.0390625  .02734375 

.0205078125 

.25  .125 

.078125  .0546875 

.041015625 

.125 

.09375  .0703125 

.0546875 

0 

.0625  .0625 

.0546875 

0 

0  .03125 

.0390625 

0 

0  0 

.015625 

r  0 
n 

IT 

and  Vn5  repeated  attempts  are  made 

to  integrate 

until,  for 

all  equations,  ^  <•  0, 

vnl 
vn2 
vn3 

\v: 

I  v 3 

I  Vn 

V„6 


When  all  6  <  0,  a  new  set  of  differences  (\7r,  i  will  be 
—  Vnfl 

computed  by  evaluating  the  derivative  of  th  e  corrector  ^7?  )? 

then  evaluating:  n+^ 


V 


,r-l 


r-1 


n+i 


v;+;  •  Vn 


r=l ,  2, 


.  .  .  5  6 


The  A-3  routine  will  check  if  a  larger  interval  h  may  be  -sed  and  still 
keep  the  errcr  within  prescribed  limits.  The  doubling  criteria  are. 

1  .  vcie  m_st  have  previously  beer,  seven  integrations 
by  AB  withci-t  a  change  in  h. 

2.  All  6  <•  -  .955*? 


dif f er enc es 


If  these  crit 

r  -iSbC 


vr; 


la  are  satisfied.,  then  h  is  doubled  and  a  new  set  cf 
v  t*ice  the  old  interval  are  computed  by: 


O  1  ! 
Vn+1  i 

1 

- 

-» 

- 1 

0 

0 

0 

0 

Vr.+1 

i 

Vn+1 

0 

4 

-4 

1 

n 

0 

i  i 

i 

! 

1 

1  i 

|  Vn+1 

i 

3 

n+i 

0 

0 

8 

-12 

6 

-1 

’  3 

-'n-*-l 

o  4  i 

^n+1 

1 

) 

1 

0 

0 

0 

16 

-32 

24 

1 

!Vn *1 

! 

v 5  i 

Vn+1  ! 

- 

0 

0 

0 

32 

-80 

! 

Vn+1 

O  6  1 

Vn+i  1 

r 

0 

0 

0 

0 

6^- 

j 

! 

! 

|V„‘l 

V,+:  -  Vr.U 

_  _ _  x 

The  \  , .  are  returned  to  the  _ser  at  this  point.  Ihe  -ew  \j  and  h 
n+ 1  Vn+1 

(if  any;  tfill  be  tried  on  the  next  step. 


/  4  .SPECIAL  CONSIDERATION  FOR  y 

1  ;n 

The  first  equation,  \ .  (located  at  W+5)  ,  is  carried  modulo  T77  for  all 

computations  within  ADBASH .  i’hi£‘  feature  allows  absolute  (only',  error  control  to 

be  used  for  this  variable.  Hie  value  Y  ,  .  ,  however,  will  not  be  modulo  2”  b_t 

will  be  Y.  (as  entered  by  the  user)  plus  AY,  by  integration. 

In  "  J-9n  y 

The  "modulo  217"  feature  restricts  the  user  to  input  Y  in  radians,  and  also 

restricts  the  derivative  rev  tire  to  accept  Y,  in  radians  (mod  Vt+s  . 

!,r. 
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.FIGURE  I-I.  ADAMS- BASHFORTH  FLOW  DIAGRAM  (2  of  6) 
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.FIGURE  I-I.  ADAMS- BASH  FORTH  FLOW  DIAGRAM  (3  of  6) 
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AD  18 


REINITIALIZE 
DIFFERENCE 
TABLE  AREA 


1 

1 

t 

X 

o 

-  X 

FIGURE  1-1.  ADAMS- BASHFORTH  FLOW  DIAGRAM  (4  of  6) 
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PREDICTOR 
-  CORRECTOR 


^'errorN, 

s  WITHIN 
PRESCRIBED 
JLIMITS  FOR 

\  / 


UPDATE  FLAG 
FOR  CONTINUOUS 
A-B  STEPS 


ADDER  7_ 


DERIV 


FIGURE  1-1 


ADAMS- BASHFORTH  FLOW  DIAGRAM  (5  of  6) 
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FIGURE  1-1.  ADAMS- BASHFORTH  FLOW  DIAGRAM  (6  of  6) 
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Description  of  W  Buffer  (Work  Region) 
Word 


W 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 


+10 
+11 
+  12 
+13 

+14 

+  15 

+16 

+  17 

+  18 


Not  Used 
T 

AT  to  be  attempted  on  next  entry  to  ADBASH 
T  at  step  start 

Integration  direction;  +  or  -  any  non  zero  number 


L 

a 


hA 

hy 


z  -I 


Y' 


Abrolute  Error  Control 


For  the  variables  listed  above  each  column,  the  following  locations  contain: 


L 

a 

X 

a 

y 

a 

z 

h 

X 

h 

h 

z 

W  +  36 

W  +  60 

W  +  84 

W  +  108 

W  +  132 

W  + 

"i56 

W  +  180 

+  37 

+  61 

+  85 

+  10’* 

+  133 

+ 

157 

+  181 

+  38 

+  62 

+  86 

+  no 

+  134 

+ 

158 

+  182 

+  39 

+  63 

+  87 

+  in 

+  135 

+ 

159 

+  183 

+  40 

+  64 

+  88 

+  112 

+  136 

+ 

160 

+  184 

+  41 

+  65 

+  89 

+  113 

+  13  7 

+ 

161 

+  185 

+  42 

+  66 

+  90 

+  114 

+  138 

+ 

162 

+  186 
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Derivatives 


W  +  187 
+  188 
+  189 
+  190 
+  191 
+  192 
+  193 
+  194 
+  195 
+  196 
+  197 
+  198 
+  199 
+  200 
+  201 
+  202 


•  i 


a 

.x 

a 


L  total  revolutions  (radians) 

L  partial  revolution  (radians) 
^  1 


Y's  at  step  start 


locations 

For  the 
contain: 

variables'  listed  above  each 

column , 

the  following 

sets  of 

© 

L 

a  a  a 

x  _ 2 _  z 

h 

X 

h 

y 

h 

2 

W  +  31 

W  +  55  W  +  79  W  +  103 

W  +  127 

W  +  151  W 

+  175 

© 

+  32 

+56  +80  +  104 

+  128 

+  152 

+  176 

@ 

+  33 

+  57  +  81i  +  105 

+  129 

+  153 

+  177 

© 

+  34 

+58  +82  +  106 

+  130 

+  154 

-t-  178 

© 

+  35 

+59  +83  +  107 

+  131 

+  155 

+  179 

1  .  After  a  Runge  Kutta  step: 

© 

Y  at  the  half  step 

(?) 

V  at  the  step  start 

© 

A  Y  <  f  the  second  half  step 

© 

A  1  of  the  first  half  step 

@ 

Y  at  the  step  start 
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2.  After  an  Adams  B^.shforth  step: 

0  Predictor 

© 

© 

© 

0  j Predictor  -  Corrector) 

The  predictor  and  corrector  constants,  B'  ,  are  stored  as  follows: 


W  +  26 

B 

o 

+  27 

B  / 
o 

+  50 

Bi 

+  51 

Bi 

+  74 

B2 

+  75 

B2 

+  98 

B3 

+  99 

B3 

+  122 

B4 

+  123 

B4 

+  146 

B5 

+  147 

B5 

+  170 

B6 

+  171 

K 

All  locations  not  mentioned  are  used  for  temporary  working  storage. 
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I -6. HOW  TO  USE  AD BASH 


The  remainder  of  this  .-Appendix  gives  specific  instructions  in  the  use 
of  the  Fver/-a.r.. 


ADBASH  has  been  designed  to  integrate  a  system  of  first  order  differential 
equations  over  some  range  of  X  from  XQto  X  fiua]..  Given  a  set  of  initial  conditions, 
ADBASH  will,  on  each  successive  jump  to  the  subroutine,  give  a  step  by  step  solution 
to  the  system  of  equations.  The  subroutine  will  automatically  control  the  step  size 
so  as  to  keep  the  error  within  the  limits  specified  by  the  user. 

To  obtain  values  of  Y^,  Y^,  . ..,  Y^,  at  discrete  values  of  X,  the  user 
should  code  his  program  to  perform  some  type  of  interpolation,  using  points  computed 
by  ADBASH  that  surround  the  point  of  interest.  This  type  of  coding  will  allow  ADBASH 
to  work  efficiently  and  use  an  optimum  step  size. 

Any  error  exits  from  ADBASH  must  be  coded  into  ADBASH  by  the  user.  An 
error  exit  is  provided  at  a  +  1H . 

T  1  ■  DERIVATIVE  SUBROUTINE  (SUPPLIED  BY  USER) 

At  each  entry,  the  derivative  subroutine  must  compute  the  derivatives  of 
the  system,  using  the  current  values  of  X  and  the  Y.,  and  store  the  results  in  the 
W  buffer.  1 

The  first  instruction  of  the  derivative  subroutine  must  be:  TJM  exit. 

The  return  instruction  of  t'  e  derivative  subroutine  must  be:  exit  JMP  0. 

Input  Output 

W  +  1  X 


+  5 

Yj  (mod  2tt; 

w  + 

187 

y; 

+  6 

Y2 

+ 

188 

Y2 

+  7 

V 

*3 

+ 

189 

Y3 

+  8 

Y4  ’ 

+ 

190 

Y4 

+  9 

V5 

+ 

191 

Y5 

+  10 

Y6 

+ 

192 

Y6 

+  11 

Y7 

+ 

193 

To  minimize  time  the  derivative  subroutine  should  perform  only  the  function 
stated  above 
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:  3.  INITIALIZATION 


Before  the  first  of  any  series  of  calls  to  AD BASH,  the  following  must 
be  defined: 


W  (work)  Region 


W  +  1 
+  4 
+  5 
+  6 
+  7 
+  8 
+  9 
+10 
+11 
+  12 
+13 
+14 
+15 
+16 
+  17 
+  18 


X 

o 

direction  of  integration,  any  +  or  -  any  non  zero  number 


Absolute 

error 

control 


Control  Word 


ADP  -  set  this  to  1/1T15 
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Derivative  Subroutine  Calling  Setup 


The  name  of  the  user  supplied  derivative  subroutine  must  be  inserted 
into  the  address  portion  of  seven  locations  within  ADBASH ,  These 
locations  are* 

ADDER 1,  ADDER2 ,  ADDER3 ,  ADDER4 , 

ADDER5 ,  ADDER6,  ADDER 7 

This  may  be  accomplished  as  follows: 

TIJ  NAME 

TJM  ADDER 1 

TJM  ADDER2 


IJM  ADDER 7 

Optional  Initialization 


RKINDT  (initial  step  size  in  fit.  pt.) 

This  is  defined  as  0.35  within  ADBASH  but  may  be 
changed  by  the  user. 

ADFERRK  (number  of  R-K  steps  taken  to  generate  one  A-B  point  » 
This  is  defined  as  4RK/1AB  within  ADBASH  but  may  be 
changed  to  any  value  below: 


1/1  115  =  2RK/1AB 
1/1  T17  =  8RK/1AB 
1/1  118  ='  16RK/1AB 
1/1  T19  =  32RK/1AB 
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1-9 


SAMPLE  FLOW  CHART 


The  following  diagram  illustrates  a  typical  application 
using  ADBASH: 


FIGURE  1-2  .  ADBASH  APPLICATION 
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I  ’ 0 .  TO  USE  RUNGE-KUTTA  INTEGRATION  ONLY 


Daring  a  series  of  jumps  to  ADBASH,  the  user  may  restrict  the  procedure 
to  use  Runge-Kutta  integration  only. 

To  use  R-K  only: 

1.  Set  ADP  =  1/1T15  before  the  first  of  a  series  of  jumps 
to  ADBASK . 

2.  Before  the  second  and  each  successive  call  to  ADBASH, 
set  ADP  as  follows: 

Option  I:  ADP  =  D/2.1B15;  This  will  permit 

ADBASH  to  decrease  or  increase  the 
step  size  as  prescribed  by  the 
error  Cunciols. 

Option  II:  ADP  =  D/2.2B15;  This  will  permit 
ADBASH  only  to  decrease  the  step 
size  as  prescribed  by  the  error 
controls . 
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